title: Learning Generative Deception Strategies in Combinatorial Masking Games

abstract: Deception is a crucial tool in the cyberdefence repertoire, enabling defenders to leverage their informational advantage to reduce the likelihood of successful attacks. One way deception can be employed is through obscuring, or masking, some of the information about how systems are configured, increasing attacker's uncertainty about their targets. We present a novel game-theoretic model of the resulting defender-attacker interaction, where the defender chooses a subset of attributes to mask, while the attacker responds by choosing an exploit to execute. The strategies of both players have combinatorial structure with complex informational dependencies, and therefore even representing these strategies is not trivial. First, we show that the problem of computing an equilibrium of the resulting zero-sum defender-attacker game can be represented as a linear program with a combinatorial number of system configuration variables and constraints, and develop a constraint generation approach for solving this problem. Next, we present a novel highly scalable approach for approximately solving such games by representing the strategies of both players as neural networks. The key idea is to represent the defender's mixed strategy using a deep neural network generator, and then using alternating gradient-descent-ascent algorithm, analogous to the training of Generative Adversarial Networks. Our experiments, as well as a case study, demonstrate the efficacy of the proposed approach.


title: FPGA-based low-cost synchronized fiber network for experimental setups in space

abstract: Custom experiment setups in physics often require control electronics to execute actions and measurements on a small time scale. When further constraints limit the experiment's environment, for example when the experiment is inside a sounding rocket, conventional network systems will not suffice those constraints because of weight, heat or budget limitations. This paper proposes a network architecture with a time resolution of less than 1 ns over a pair of plastic fibers while using low-cost commercial hardware. The plastic fibers in comparison to copper fibers have a low weight and additionally can isolate parts of the setup galvanically. Data rates of 40 Mbit/s enable the network to transfer large amounts of measurements and configuration data over the network. Proof-of-concept implementations of network endpoints and switches on small FPGAs are analyzed in terms of synchronicity, data rate and resource usage. Using commercial parts the resolution of 1 ns is reached with a standard deviation of less than 100 ps. Compared to a copper wire implementation the weight is reduced by about one order of magnitude. With its low weight at a low cost, the network is useful in space or laboratory setups which require high time resolution.


title: Obstacle-aware Waypoint Generation for Long-range Guidance of Deep-Reinforcement-Learning-based Navigation Approaches

abstract: Navigation of mobile robots within crowded environments is an essential task in various use cases, such as delivery, health care, or logistics. Deep Reinforcement Learning (DRL) emerged as an alternative method to replace overly conservative approaches and promises more efficient and flexible navigation. However, Deep Reinforcement Learning is limited to local navigation due to its myopic nature. Previous research works proposed various ways to combine Deep Reinforcement Learning with conventional methods but a common problem is the complexity of highly dynamic environments due to the unpredictability of humans and other objects within the environment. In this paper, we propose a hierarchical waypoint generator, which considers moving obstacles and thus generates safer and more robust waypoints for Deep-Reinforcement-Learning-based local planners. Therefore, we utilize Delaunay Triangulation to encode obstacles and incorporate an extended hybrid A-Star approach to efficiently search for an optimal solution in the time-state space. We compared our waypoint generator against two baseline approaches and outperform them in terms of safety, efficiency, and robustness.


title: The Transfer Matrix Method and The Theory of Finite Periodic Systems. From Heterostructures to Superlattices

abstract: Long-period systems and superlattices, with additional periodicity, have new effects on the energy spectrum and wave functions. Most approaches adjust theories for infinite systems, which is acceptable for large but not small number of unit cells $n$. In the past 30 years, a theory based entirely on transfer matrices was developed, where the finiteness of $n$ is an essential condition. The theory of finite periodic systems (TFPS) is also valid for any number of propagating modes, and arbitrary potential profiles (or refractive indices). We review this theory, the transfer matrix definition, symmetry properties, group representations, and relations with the scattering amplitudes. We summarize the derivation of multichannel matrix polynomials (which reduce to Chebyshev polynomials in the one-propagating mode limit), the analytical formulas for resonant states, energy eigenvalues, eigenfunctions, parity symmetries, and discrete dispersion relations, for superlattices with different confinement characteristics. After showing the inconsistencies and limitations of hybrid approaches that combine the transfer-matrix method with Floquet's theorem, we review some applications of the TFPS to multichannel negative resistance, ballistic transistors, channel coupling, spintronics, superluminal, and optical antimatter effects. We review two high-resolution experiments using superlattices: tunneling time in photonic band-gap and optical response of blue-emitting diodes, and show extremely accurate theoretical predictions.


title: Turn-to-Diarize: Online Speaker Diarization Constrained by Transformer Transducer Speaker Turn Detection

abstract: In this paper, we present a novel speaker diarization system for streaming on-device applications. In this system, we use a transformer transducer to detect the speaker turns, represent each speaker turn by a speaker embedding, then cluster these embeddings with constraints from the detected speaker turns. Compared with conventional clustering-based diarization systems, our system largely reduces the computational cost of clustering due to the sparsity of speaker turns. Unlike other supervised speaker diarization systems which require annotations of time-stamped speaker labels for training, our system only requires including speaker turn tokens during the transcribing process, which largely reduces the human efforts involved in data collection.


title: Bayesian Inference of a Social Graph with Trace Feasibility Guarantees

abstract: Network inference is the process of deciding what is the true unknown graph underlying a set of interactions between nodes. There is a vast literature on the subject, but most known methods have an important drawback: the inferred graph is not guaranteed to explain every interaction from the input trace. We consider this an important issue since such inferred graph cannot be used as input for applications that require a reliable estimate of the true graph. On the other hand, a graph having trace feasibility guarantees can help us better understand the true (hidden) interactions that may have taken place between nodes of interest. The inference of such graph is the goal of this paper. Firstly, given an activity log from a social network, we introduce a set of constraints that take into consideration all the hidden paths that are possible between the nodes of the trace, given their timestamps of interaction. Then, we develop a nontrivial modification of the Expectation-Maximization algorithm by Newman [1], that we call Constrained-EM, which incorporates the constraints and a set of auxiliary variables into the inference process to guide it towards the feasibility of the trace. Experimental results on real-world data from Twitter confirm that Constrained-EM generates a posterior distribution of graphs that explains all the events observed in the trace while presenting the desired properties of a scale-free, small-world graph. Our method also outperforms established methods in terms of feasibility and quality of the inferred graph.


title: The equivariant coarse Baum-Connes conjecture for metric spaces with proper group actions

abstract: The equivariant coarse Baum-Connes conjecture interpolates between the Baum-Connes conjecture for a discrete group and the coarse Baum-Connes conjecture for a proper metric space. In this paper, we study this conjecture under certain assumptions. More precisely, assume that a countable discrete group $\Gamma$ acts properly and isometrically on a discrete metric space $X$ with bounded geometry, not necessarily cocompact. We show that if the quotient space $X/\Gamma$ admits a coarse embedding into Hilbert space and $\Gamma$ is amenable, and that the $\Gamma$-orbits in $X$ are uniformly equivariantly coarsely equivalent to each other, then the equivariant coarse Baum-Connes conjecture holds for $(X, \Gamma)$. Along the way, we prove a $K$-theoretic amenability statement for the $\Gamma$-space $X$ under the same assumptions as above, namely, the canonical quotient map from the maximal equivariant Roe algebra of $X$ to the reduced equivariant Roe algebra of $X$ induces an isomorphism on $K$-theory.


title: A Learned Stereo Depth System for Robotic Manipulation in Homes

abstract: We present a passive stereo depth system that produces dense and accurate point clouds optimized for human environments, including dark, textureless, thin, reflective and specular surfaces and objects, at 2560x2048 resolution, with 384 disparities, in 30 ms. The system consists of an algorithm combining learned stereo matching with engineered filtering, a training and data-mixing methodology, and a sensor hardware design. Our architecture is 15x faster than approaches that perform similarly on the Middlebury and Flying Things Stereo Benchmarks. To effectively supervise the training of this model, we combine real data labelled using off-the-shelf depth sensors, as well as a number of different rendered, simulated labeled datasets. We demonstrate the efficacy of our system by presenting a large number of qualitative results in the form of depth maps and point-clouds, experiments validating the metric accuracy of our system and comparisons to other sensors on challenging objects and scenes. We also show the competitiveness of our algorithm compared to state-of-the-art learned models using the Middlebury and FlyingThings datasets.


title: A high-order shock capturing discontinuous Galerkin-finite-difference hybrid method for GRMHD

abstract: We present a discontinuous Galerkin-finite-difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple. An unlimited discontinuous Galerkin candidate solution is computed for the next time step. If the candidate solution is inadmissible, the time step is retaken using robust finite-difference methods. Because of its a posteriori nature, the hybrid scheme inherits the best properties of both methods. It is high-order with exponential convergence in smooth regions, while robustly handling discontinuities. We give a detailed description of how we transfer the solution between the discontinuous Galerkin and finite-difference solvers, and the troubled-cell indicators necessary to robustly handle slow-moving discontinuities and simulate magnetized neutron stars. We demonstrate the efficacy of the proposed method using a suite of standard and very challenging 1d, 2d, and 3d relativistic magnetohydrodynamics test problems. The hybrid scheme is designed from the ground up to efficiently simulate astrophysical problems such as the inspiral, coalescence, and merger of two neutron stars.


title: Optical Characterization of Disordered Yb-doped Silica Glass Anderson Localizing Optical Fiber

abstract: We investigate and report the optical and laser characteristics of a ytterbium-doped transverse Anderson localizing optical fiber to develop a fundamental understanding of the light propagation, generation, and amplification processes in this novel fiber. Ultimately, the goal based on the measurements and calculations conducted herein is to design and build a random fiber laser with a highly directional beam. The measurements are based on certain observations of the laser pump propagation and amplified spontaneous emission generation in this fiber. Judicious approximations are used in the propagation equations to obtain the relevant desired parameters in simple theoretical fits to the experimental observations, without resorting to speculations based on the intended construction from the fiber preform.


title: Treatment Effects in Market Equilibrium

abstract: Policy-relevant treatment effect estimation in a marketplace setting requires taking into account both the direct benefit of the treatment and any spillovers induced by changes to the market equilibrium. The standard way to address these challenges is to evaluate interventions via cluster-randomized experiments, where each cluster corresponds to an isolated market. This approach, however, cannot be used when we only have access to a single market (or a small number of markets). Here, we show how to identify and estimate policy-relevant treatment effects using a unit-level randomized trial run within a single large market. A standard Bernoulli-randomized trial allows consistent estimation of direct effects, and of treatment heterogeneity measures that can be used for welfare-improving targeting. Estimating spillovers - as well as providing confidence intervals for the direct effect - requires estimates of price elasticities, which we provide using an augmented experimental design. Our results rely on all spillovers being mediated via the (observed) prices of a finite number of traded goods, and the market power of any single unit decaying as the market gets large. We illustrate our results using a simulation calibrated to a conditional cash transfer experiment in the Philippines.


title: Decentralized Control of Two Agents with Nested Accessible Information

abstract: In this paper, we investigate a decentralized stochastic control problem with two agents, where a part of the memory of the second agent is also available to the first agent at each instance of time. We derive a structural form for optimal control strategies which allows us to restrict their domain to a set which does not grow in size with time. We also present a dynamic programming (DP) decomposition which can utilize our results to derive optimal strategies for arbitrarily long time horizons. Since obtaining optimal control strategies by solving this DP decomposition is computationally intensive, we present potential resolutions in the form of simplified strategies by imposing additional conditions on our model, and an approximation technique which can be used to implement our results with a bounded loss of optimality.


title: Deep Learning with Kernel Flow Regularization for Time Series Forecasting

abstract: Long Short-Term Memory (LSTM) neural networks have been widely used for time series forecasting problems. However, LSTMs are prone to overfitting and performance reduction during test phases. Several different regularization techniques have been shown in literature to prevent overfitting problems in neural networks. In this paper, first, we introduce application of kernel flow methods for time series forecasting in general. Afterward, we examine the effectiveness of applying kernel flow regularization on LSTM layers to avoid overfitting problems. We describe a regularization method by applying kernel flow loss function on LSTM layers. In experimental results, we show that kernel flow outperforms baseline models on time series forecasting benchmarks. We also compare the effect of dropout and kernel flow regularization techniques on LSTMs. The experimental results illustrate that kernel flow achieves similar regularization effect to dropout. It also shows that the best results is obtained using both kernel flow and dropout regularizations with early stopping on LSTM layers on some time series datasets (e.g. power-load demand forecasts).


title: A comprehensive revisit of select Galileo/NIMS observations of Europa

abstract: The Galileo Near Infrared Mapping Spectrometer (NIMS) collected spectra of Europa in the 0.7-5.2 $\mu$m wavelength region, which have been critical to improving our understanding of the surface composition of this moon. However, most of the work done to get constraints on abundances of species like water ice, hydrated sulfuric acid, hydrated salts and oxides have used proxy methods, such as absorption strength of spectral features or fitting a linear mixture of laboratory generated spectra. Such techniques neglect the effect of parameters degenerate with the abundances, such as the average grain-size of particles, or the porosity of the regolith. In this work we revisit three Galileo NIMS spectra, collected from observations of the trailing hemisphere of Europa, and use a Bayesian inference framework, with the Hapke reflectance model, to reassess Europa's surface composition. Our framework has several quantitative improvements relative to prior analyses: (1) simultaneous inclusion of amorphous and crystalline water ice, sulfuric-acid-octahydrate (SAO), CO$_2$, and SO$_2$; (2) physical parameters like regolith porosity and radiation-induced band-center shift; and (3) tools to quantify confidence in the presence of each species included in the model, constrain their parameters, and explore solution degeneracies. We find that SAO strongly dominates the composition in the spectra considered in this study, while both forms of water ice are detected at varying confidence levels. We find no evidence of either CO$_2$ or SO$_2$ in any of the spectra; we further show through a theoretical analysis that it is highly unlikely that these species are detectable in any 1-2.5 $\mu$m Galileo NIMS data.


title: Effect of Photon Counting Shot Noise on Total Internal Reflection Microscopy

abstract: Total internal reflection microscopy (TIRM) measures changes in the distance between a colloidal particle and a transparent substrate by measuring the intensity of light scattered by the particle when it is illuminated by the evanescent field that is created from light totally internally reflected at the substrate interface. From these measurements, the height-dependent effective potential $\varphi(z)$ between the colloidal particle and the substrate can be measured. The spatial resolution with which TIRM can resolve the height $z$ and effective potential $\varphi(z)$ is limited by the intrinsic shot noise of the photon counting process used to measure the scattered light intensity. We develop a model to determine the spatial resolution with which TIRM can measure $\varphi(z)$ and verify its validity with simulations and experiments. We further establish the critical role of photon-counting statistics and the intensity integration time $\tau$ in TIRM measurements, which is a trade-off between narrowing the width of the photon counting distribution and capturing the instantaneous position of the probe particle.


title: The $\Pi^1_2$ Consequences of a Theory

abstract: We develop the abstract framework for a proof-theoretic analysis of theories with scope beyond ordinal numbers, resulting in an analog of Ordinal Analysis aimed at the study of theorems of complexity $\Pi^1_2$. This is done by replacing the use of ordinal numbers by particularly uniform, wellfoundedness preserving functors in the category of linear orders. Generalizing the notion of a proof-theoretic ordinal, we define the functorial $\Pi^1_2$ norm of a theory and prove its existence and uniqueness for $\Pi^1_2$-sound theories. From this, we further abstract a definition of the $\Sigma^1_2$- and $\Pi^1_2$-soundness ordinals of a theory; these quantify, respectively, the maximum strength of true $\Sigma^1_2$ theorems and minimum strength of false $\Pi^1_2$ theorems of a given theory. We study these ordinals, developing a proof-theoretic classification theory for recursively enumerable extensions of $\mathsf{ACA}_0$ Using techniques from infinitary and categorical proof theory, generalized recursion theory, constructibility, and forcing, we prove that an admissible ordinal is the $\Pi^1_2$-soundness ordinal of some recursively enumerable extension of $\mathsf{ACA}_0$ if and only if it is not parameter-free $\Sigma^1_1$-reflecting. We show that the $\Sigma^1_2$-soundness ordinal of $\mathsf{ACA}_0$ is $\omega_1^{ck}$ and characterize the $\Sigma^1_2$-soundness ordinals of recursively enumerable, $\Sigma^1_2$-sound extensions of $\Pi^1_1{-}\mathsf{CA}_0$.


title: Formation of polar terrestrial circumbinary planets

abstract: All circumbinary planets currently detected are in orbits that are almost coplanar to the binary orbit. While misaligned circumbinary planets are more difficult to detect, observations of polar aligned circumbinary gas and debris disks around eccentric binaries suggest that polar planet formation may be possible. A polar aligned planet has a stable orbit that is inclined by 90 degrees to the orbital plane of the binary with an angular momentum vector that is aligned to the binary eccentricity vector. With n- body simulations we model polar terrestrial planet formation using hydrodynamic gas disk simulations to motivate the initial particle distribution. Terrestrial planet formation around an eccentric binary is more likely in a polar alignment than in a coplanar alignment. Similar planetary systems form in a polar alignment around an eccentric binary and a coplanar alignment around a circular binary. The polar planetary systems are stable even with the effects of general relativity. Planetary orbits around an eccentric binary exhibit tilt and eccentricity oscillations at all inclinations, however, the oscillations are larger in the coplanar case than the polar case. We suggest that polar aligned terrestrial planets will be found in the future.


title: Modeling Dynamic Attributes for Next Basket Recommendation

abstract: Traditional approaches to next item and next basket recommendation typically extract users' interests based on their past interactions and associated static contextual information (e.g. a user id or item category). However, extracted interests can be inaccurate and become obsolete. Dynamic attributes, such as user income changes, item price changes (etc.), change over time. Such dynamics can intrinsically reflect the evolution of users' interests. We argue that modeling such dynamic attributes can boost recommendation performance. However, properly integrating them into user interest models is challenging since attribute dynamics can be diverse such as time-interval aware, periodic patterns (etc.), and they represent users' behaviors from different perspectives, which can happen asynchronously with interactions. Besides dynamic attributes, items in each basket contain complex interdependencies which might be beneficial but nontrivial to effectively capture. To address these challenges, we propose a novel Attentive network to model Dynamic attributes (named AnDa). AnDa separately encodes dynamic attributes and basket item sequences. We design a periodic aware encoder to allow the model to capture various temporal patterns from dynamic attributes. To effectively learn useful item relationships, intra-basket attention module is proposed. Experimental results on three real-world datasets demonstrate that our method consistently outperforms the state-of-the-art.


title: A note on mean field games of controls with state constraints: existence of mild solutions

abstract: We show the existence of "mild solutions" for a first-order mean field game of controls under the state constraint that trajectories be confined in a closed and bounded set in euclidean space. This extends the results of Cannarsa and Capuani to the case of a mean field game of controls. Our controls are velocities and we find that the existence of an equilibrium is complicated by the requirement that they should have enough regularity. We solve this by imposing a small Lipschitz constant on the dependence of the Lagrangian on the joint measure of states and controls, and showing that regular paths can be approximated within the same class of functions despite the constraint.


title: Sparse multi-reference alignment: sample complexity and computational hardness

abstract: Motivated by the problem of determining the atomic structure of macromolecules using single-particle cryo-electron microscopy (cryo-EM), we study the sample and computational complexities of the sparse multi-reference alignment (MRA) model: the problem of estimating a sparse signal from its noisy, circularly shifted copies. Based on its tight connection to the crystallographic phase retrieval problem, we establish that if the number of observations is proportional to the square of the variance of the noise, then the sparse MRA problem is statistically feasible for sufficiently sparse signals. To investigate its computational hardness, we consider three types of computational frameworks: projection-based algorithms, bispectrum inversion, and convex relaxations. We show that a state-of-the-art projection-based algorithm achieves the optimal estimation rate, but its computational complexity is exponential in the sparsity level. The bispectrum framework provides a statistical-computational trade-off: it requires more observations (so its estimation rate is suboptimal), but its computational load is provably polynomial in the signal's length. The convex relaxation approach provides polynomial time algorithms (with a large exponent) that recover sufficiently sparse signals at the optimal estimation rate. We conclude the paper by discussing potential statistical and algorithmic implications for cryo-EM.


title: The time-identity tradeoff

abstract: Distinguishability plays a major role in quantum and statistical physics. When particles are identical their wave function must be either symmetric or antisymmetric under permutations and the number of microscopic states, which determines entropy, is counted up to permutations. When the particles are distinguishable, wavefunctions have no symmetry and each permutation is a different microstate. This binary and discontinuous classification raises a few questions: one may wonder what happens if particles are almost identical, or when the property that distinguishes between them is irrelevant to the physical interactions in a given system. Here I sketch a general answer to these questions. For any pair of non-identical particles there is a timescale, $\tau_d$, required for a measurement to resolve the differences between them. Below $\tau_d$, particles seem identical, above it - different, and the uncertainty principle provides a lower bound for $\tau_d$. Thermal systems admit a conjugate temperature scale, $T_d$. Above this temperature the system appears to equilibrate before it resolves the differences between particles, below this temperature the system identifies these differences before equilibration. As the physical differences between particles decline towards zero, $\tau_d \to \infty$ and $T_d \to 0$.


title: Relaxation approach for learning neural network regularizers for a class of identification problems

abstract: The present paper deals with the data-driven design of regularizers in the form of artificial neural networks, for solving certain inverse problems formulated as optimal control problems. These regularizers aim at improving accuracy, wellposedness or compensating uncertainties for a given class of optimal control problems (inner-problems). Parameterized as neural networks, their weights are chosen in order to reduce a misfit between data and observations of the state solution of the inner-optimal control problems. Learning these weights constitutes the outer-problem. Based on necessary first-order optimality conditions for the inner-problems, a relaxation approach is proposed in order to implement efficient solving of these inner-problems, namely the forward operator of the outer-problem. Optimality conditions are derived for the latter, and are implemented in numerical illustrations dealing with the inverse conductivity problem. The numerical tests show the feasibility of the relaxation approach, first for rediscovering standard $L^2$-regularizers, and next for designing regularizers that compensate unknown noise on the observed state of the inner-problem.


title: Reflected spectroscopy of small exoplanets I: determining the atmospheric composition of sub-Neptunes planets

abstract: Direct imaging of widely separated exoplanets from space will obtain their reflected light spectra and measure their atmospheric properties, and small and temperate planets will be the focus of the next generation telescopes. In this work, we used our Bayesian retrieval algorithm ExoReL$^{\Re}$ to determine the constraints on the atmospheric properties of sub-Neptune planets from observations taken with a HabEx-like telescope. Small and temperate planets may have a non-H$_2$-dominated atmosphere, and therefore, we introduced the compositional analysis technique in our framework to explore the bulk atmospheric chemistry composition without any prior knowledge about it. We have developed a novel set of prior functions for the compositional analysis free parameters. We compared the performances of the framework with the flat prior and the novel prior and we reported a better performance when using the novel priors set. We found that the retrieval algorithm can not only identify the dominant gas of the atmosphere but also to constraint other less abundant gases with high statistical confidence without any prior information on the composition. The results presented here demonstrates that reflected light spectroscopy can characterize small exoplanets with diverse atmospheric composition. The Bayesian framework should be applied to design the instrument and the observation plan of exoplanet direct imaging experiments in the future.


title: The fourth virial coefficient for hard spheres in even dimension

abstract: The fourth virial coefficient is calculated exactly for a fluid of hard spheres in even dimensions. For this purpose the complete star cluster integral is expressed as the sum of two three-folded integrals only involving spherical angular coordinates. These integrals are solved anallytically for any even dimension d and working with existing expressions for the other terms of the fourth cluster integral we obtain an expression for the fourth virial coefficient B4(d) for even d that sums a finite number of simple terms, with the number of terms increasing with d.


title: Deep Reinforcement Learning-Based Long-Range Autonomous Valet Parking for Smart Cities

abstract: In this paper, to reduce the congestion rate at the city center and increase the quality of experience (QoE) of each user, the framework of long-range autonomous valet parking (LAVP) is presented, where an Autonomous Vehicle (AV) is deployed in the city, which can pick up, drop off users at their required spots, and then drive to the car park out of city center autonomously. In this framework, we aim to minimize the overall distance of the AV, while guarantee all users are served, i.e., picking up, and dropping off users at their required spots through optimizing the path planning of the AV and number of serving time slots. To this end, we first propose a learning based algorithm, which is named as Double-Layer Ant Colony Optimization (DL-ACO) algorithm to solve the above problem in an iterative way. Then, to make the real-time decision, while consider the dynamic environment (i.e., the AV may pick up and drop off users from different locations), we further present a deep reinforcement learning (DRL) based algorithm, which is known as deep Q network (DQN). The experimental results show that the DL-ACO and DQN-based algorithms both achieve the considerable performance.


title: Constraining Time Dependent Dark Matter Signals from the Sun

abstract: Dark matter (DM) particles captured by the Sun can produce high energy electrons outside the Sun through annihilating into meta-stable mediators. The corresponding cosmic-ray electron signals observed by the space-based experiments will be time dependent due to the orbital motion of the space-based detectors. The shape of this time dependence is predictable given the orbital information of the detectors. Since the high-energy CR electron (with energy E>100 GeV) fluxes are expected to be constant in time, non-observation of such time variation can be used to place upper limits on the DM annihilation cross section. We analyze the time dependence of dark matter cosmic-ray signals in three space-based experiments: AMS-02, DAMPE and CALET. Under the assumption that no time dependent signal is observed, we derive the 95% C.L. exclusion limits on the signal strength from the current data. We map our limits onto the parameter space of the dark photon model and find that the constraints are comparable with that derived from the supernova SN1987A.


title: Simulation of plasma emission in magnetized plasmas

abstract: The recent Parker Solar Probe (PSP) observations of type III radio bursts show that the effects of finite background magnetic field can be an important factor in the interpretation of data. In the present paper, the effects of background magnetic field on the plasma emission process, which is believed to be the main emission mechanism for solar coronal and interplanetary type III radio bursts, are investigated by means of the particle-in-cell simulation method. The effects of ambient magnetic field are systematically surveyed by varying the ratio of plasma frequency to electron gyro-frequency. The present study shows that for a sufficiently strong ambient magnetic field, the wave-particle interaction processes lead to a highly field-aligned longitudinal mode excitation and anisotropic electron velocity distribution function, accompanied by a significantly enhanced plasma emission at the second harmonic plasma frequency. For such a case, the polarization of the harmonic emission is almost entirely in the sense of extraordinary mode. On the other hand, for moderate strengths of the ambient magnetic field, the interpretation of the simulation result is less than clear. The underlying nonlinear mode coupling processes indicate that to properly understand and interpret the simulation results require sophisticated analyses involving interactions among magnetized plasma normal modes including the two transverse modes of the magneto-active plasma, namely, extraordinary and ordinary modes, as well as electron-cyclotron-whistler, plasma oscillation, and upper-hybrid modes. At present, a nonlinear theory suitable for quantitatively analyzing such complex mode-coupling processes in magnetized plasmas is incomplete, which calls for further theoretical research, but the present simulation results could provide a guide for future theoretical efforts.


title: A Small Variation of the Circular Hodograph Theorem and the Best Elliptical Trajectory of the Planets

abstract: A small variation of the circular shape of the hodograph theorem states that for every elliptical solution of the two-body problem, it is possible to find an appropriate inertial frame such that the speed of the bodies is constant. We use this result and data from the NASA JPL Horizon Web Interface to find the best fitting ellipse for the trajectory of Mercury, Venus, Earth, Mars, and Jupiter. The process requires us to find procedures to obtain the plane and ellipse that best fit a collection of points in space. We show that if we aim for the plane that minimizes the sum of the square distances from the given points to the unknown plane, we obtain three planes that appear to divide the set of points equally into octants, one of these being our desired plane of best fit. We provide a detailed proof of the hodograph theorem.


title: WebFlow: Scalable and Decentralized Routing for Payment Channel Networks with High Resource Utilization

abstract: Payment channel networks (PCNs) have been designed and utilized to address the scalability challenge and throughput limitation of blockchains. Routing is a core problem of PCNs. An ideal PCN routing method needs to achieve 1) high scalability that can maintain low per-node memory and communication cost for large PCNs, 2) high resource utilization of payment channels, and 3) the privacy of users. However, none of the existing PCN systems consider all these requirements. In this work, we propose WebFlow, a distributed routing solution for PCNs, which only requires each user to maintain localized information and can be used for massive-scale networks with high resource utilization. We make use of two distributed data structures: multi-hop Delaunay triangulation (MDT) originally proposed for wireless networks and our innovation called distributed Voronoi diagram. We propose new protocols to generate a virtual Euclidean space in order to apply MDT to PCNs and use the distributed Voronoi diagram to enhance routing privacy. We conduct extensive simulations and prototype implementation to further evaluate WebFlow. The results using real and synthetic PCN topologies and transaction traces show that WebFlow can achieve extremely low per-node overhead and a high success rate compared to existing methods.


title: SLO beyond the Hardware Isolation Limits

abstract: Performance isolation is a keystone for SLO guarantees with shared resources in cloud and datacenter environments. To meet SLO requirements, the state of the art relies on hardware QoS support (e.g., Intel RDT) to allocate shared resources such as last-level caches and memory bandwidth for co-located latency-critical applications. As a result, the number of latency-critical applications that can be deployed on a physical machine is bounded by the hardware allocation capability. Unfortunately, such hardware capability is very limited. For example, Intel Xeon E5 v3 processors support at most four partitions for last-level caches, i.e., at most four applications can have dedicated resource allocation. This paper discusses the feasibility and unexplored challenges of providing SLO guarantees beyond the limits of hardware capability. We present CoCo to show the feasibility and the benefits. CoCo schedules applications to time-share interference-free partitions as a transparent software layer. Our evaluation shows that CoCo outperforms non-partitioned and round-robin approaches by up to 9x and 1.2x.


title: Performance Improvement of Dimmable OFDM-Visible Light Communication using Subcarrier Index Modulation and Reed Solomon Encoding

abstract: In this paper, we propose a new subcarrier index modulation scheme for orthogonal frequency division multiplexing (OFDM), which incorporates the Reed-Solomon (RS) encoding in a visible light communication (VLC) system. In this scheme, the incoming bits are first encoded using the RS encoder, then a set of symbols in the resulting RS codeword are punctured, and the remaining symbols are modulated and mapped onto the OFDM subcarriers. This system is referred to as RS-OFDM-IM. Unlike the traditional subcarrier index modulation (SIM) schemes, the proposed scheme operates based on conveying extra information by inactivating the selected subcarriers, which facilitates simultaneous clipping noise reduction and spectral efficiency enhancement in OFDM-VLC. The bit error rate (BER) and throughput of the proposed technique is theoretically and numerically analyzed. The simulation results show the superiority of the proposed technique as compared to the coded DCO-OFDM without SIM and the classical SIM in OFDM-VLC in a system with practical clipping conditions.


title: Exact Learning of Qualitative Constraint Networks from Membership Queries

abstract: A Qualitative Constraint Network (QCN) is a constraint graph for representing problems under qualitative temporal and spatial relations, among others. More formally, a QCN includes a set of entities, and a list of qualitative constraints defining the possible scenarios between these entities. These latter constraints are expressed as disjunctions of binary relations capturing the (incomplete) knowledge between the involved entities. QCNs are very effective in representing a wide variety of real-world applications, including scheduling and planning, configuration and Geographic Information Systems (GIS). It is however challenging to elicit, from the user, the QCN representing a given problem. To overcome this difficulty in practice, we propose a new algorithm for learning, through membership queries, a QCN from a non expert. In this paper, membership queries are asked in order to elicit temporal or spatial relationships between pairs of temporal or spatial entities. In order to improve the time performance of our learning algorithm in practice, constraint propagation, through transitive closure, as well as ordering heuristics, are enforced. The goal here is to reduce the number of membership queries needed to reach the target QCN. In order to assess the practical effect of constraint propagation and ordering heuristics, we conducted several experiments on randomly generated temporal and spatial constraint network instances. The results of the experiments are very encouraging and promising.


title: Convergence of Langevin-Simulated Annealing algorithms with multiplicative noise

abstract: We study the convergence of Langevin-Simulated Annealing type algorithms with multiplicative noise, i.e. for $V : \mathbb{R}^d \to \mathbb{R}$ a potential function to minimize, we consider the stochastic equation $dY_t = - \sigma \sigma^\top \nabla V(Y_t) dt + a(t)\sigma(Y_t)dW_t + a(t)^2\Upsilon(Y_t)dt$, where $(W_t)$ is a Brownian motion, where $\sigma : \mathbb{R}^d \to \mathcal{M}_d(\mathbb{R})$ is an adaptive (multiplicative) noise, where $a : \mathbb{R}^+ \to \mathbb{R}^+$ is a function decreasing to $0$ and where $\Upsilon$ is a correction term. This setting can be applied to optimization problems arising in Machine Learning. The case where $\sigma$ is a constant matrix has been extensively studied however little attention has been paid to the general case. We prove the convergence for the $L^1$-Wasserstein distance of $Y_t$ and of the associated Euler-scheme $\bar{Y}_t$ to some measure $\nu^\star$ which is supported by $\text{argmin}(V)$ and give rates of convergence to the instantaneous Gibbs measure $\nu_{a(t)}$ of density $\propto \exp(-2V(x)/a(t)^2)$. To do so, we first consider the case where $a$ is a piecewise constant function. We find again the classical schedule $a(t) = A\log^{-1/2}(t)$. We then prove the convergence for the general case by giving bounds for the Wasserstein distance to the stepwise constant case using ergodicity properties.


title: The Ice Chamber for Astrophysics-Astrochemistry (ICA): A New Experimental Facility for Ion Impact Studies of Astrophysical Ice Analogues

abstract: The Ice Chamber for Astrophysics-Astrochemistry (ICA) is a new laboratory end-station located at the Institute for Nuclear Research (Atomki) in Debrecen, Hungary. The ICA has been specifically designed for the study of the physico-chemical properties of astrophysical ice analogues and their chemical evolution when subjected to ionising radiation and thermal processing. The ICA is an ultra-high vacuum compatible chamber containing a series of IR-transparent substrates mounted in a copper holder connected to a closed-cycle cryostat capable of being cooled down to 20 K, itself mounted on a 360{\deg} rotation stage and a z-linear manipulator. Ices are deposited onto the substrates via background deposition of dosed gases. Ice structure and chemical composition are monitored by means of FTIR absorbance spectroscopy in transmission mode, although use of reflectance mode is possible by using metallic substrates. Pre-prepared ices may be processed in a variety of ways. A 2 MV Tandetron accelerator is capable of delivering a wide variety of high-energy ions into the ICA, which simulates ice processing by cosmic rays, the solar wind, or magnetospheric ions. The ICA is also equipped with an electron gun which may be used for electron impact radiolysis of ices. Thermal processing of both deposited and processed ices may be monitored by means of both FTIR spectroscopy and quadrupole mass spectrometry. In this paper, we provide a detailed description of the ICA set-up, as well as an overview of preliminary results obtained and future plans.


title: Matricial Archimedean order unit spaces and quantum correlations

abstract: We introduce a matricial analogue of an Archimedean order unit space, which we call a $k$-AOU space. We develop the category of $k$-AOU spaces and $k$-positive maps and exhibit functors from this category to the category of operator systems and completely positive maps. We also demonstrate the existence of injective envelopes and C*-envelopes in the category of $k$-AOU spaces. Finally, we show that finite-dimensional quantum correlations can be characterized in terms of states on finite-dimensional $k$-AOU spaces. Combined with previous work, this yields a reformulation of Tsirelson's conjecture in terms of operator systems and $k$-AOU spaces.


title: A Multi-Agent Deep Reinforcement Learning Coordination Framework for Connected and Automated Vehicles at Merging Roadways

abstract: The steady increase in the number of vehicles operating on the highways continues to exacerbate congestion, accidents, energy consumption, and greenhouse gas emissions. Emerging mobility systems, e.g., connected and automated vehicles (CAVs), have the potential to directly address these issues and improve transportation network efficiency and safety. In this paper, we consider a highway merging scenario and propose a framework for coordinating CAVs such that stop-and-go driving is eliminated. We use a decentralized form of the actor-critic approach to deep reinforcement learning$-$multi-agent deep deterministic policy gradient. We demonstrate the coordination of CAVs through numerical simulations and show that a smooth traffic flow is achieved by eliminating stop-and-go driving. Videos and plots of the simulation results can be found at this supplemental $\href{https://sites.google.com/view/ud-ids-lab/MADRL}{\text{site}}$.


title: Modeling calcium dynamics in neurons with endoplasmic reticulum: existence, uniqueness and an implicit-explicit finite element scheme

abstract: Like many other biological processes, calcium dynamics in neurons containing an endoplasmic reticulum are governed by diffusion-reaction equations on interface-separated domains. Interface conditions are typically described by systems of ordinary differential equations that provide fluxes across the interfaces. Using the calcium model as an example of this class of ODE-flux boundary interface problems, we prove the existence, uniqueness and boundedness of the solution by applying comparison theorem, fundamental solution of the parabolic operator and a strategy used in Picard's existence theorem. Then we propose and analyze an efficient implicit-explicit finite element scheme which is implicit for the parabolic operator and explicit for the nonlinear terms. We show that the stability does not depend on the spatial mesh size. Also the optimal convergence rate in $H^1$ norm is obtained. Numerical experiments illustrate the theoretical results.


title: Electron Irradiation and Thermal Chemistry Studies of Interstellar and Planetary Ice Analogues at the ICA Astrochemistry Facility

abstract: The modelling of molecular excitation and dissociation processes relevant to astrochemistry requires the validation of theories by comparison with data generated from laboratory experimentation. The newly commissioned Ice Chamber for Astrophysics-Astrochemistry (ICA) allows for the study of astrophysical ice analogues and their evolution when subjected to energetic processing, thus simulating the processes and alterations interstellar icy grain mantles and icy outer Solar System bodies undergo. ICA is an ultra-high vacuum compatible chamber containing a series of IR-transparent substrates upon which the ice analogues may be deposited at temperatures of down to 20 K. Processing of the ices may be performed in one of three ways: (i) ion impacts with projectiles delivered by a 2 MV Tandetron-type accelerator, (ii) electron irradiation from a gun fitted directly to the chamber, and (iii) thermal processing across a temperature range of 20-300 K. The physico-chemical evolution of the ices is studied in situ using FTIR absorbance spectroscopy and quadrupole mass spectrometry. In this paper, we present an overview of the ICA facility with a focus on characterising the electron beams used for electron impact studies, as well as reporting the preliminary results obtained during electron irradiation and thermal processing of selected ices.


title: Removing the deadwood from DFT/MRCI wave functions: The p-DFT/MRCI method

abstract: The combined density functional theory and multireference configuration interaction (DFT/MRCI) method is a powerful tool for the calculation of excited electronic states of large molecules. There exists, however, a large amount of superfluous configurations in a typical DFT/MRCI wave function. We show that this deadwood may be effectively removed using a simple configuration pruning algorithm based on second-order Epstein-Nesbet perturbation theory. The resulting method, which we denote p-DFT/MRCI, is shown to result in orders of magnitude saving in computational timings, while retaining the accuracy of the original DFT/MRCI method.


title: Theory of overparametrization in quantum neural networks

abstract: The prospect of achieving quantum advantage with Quantum Neural Networks (QNNs) is exciting. Understanding how QNN properties (e.g., the number of parameters $M$) affect the loss landscape is crucial to the design of scalable QNN architectures. Here, we rigorously analyze the overparametrization phenomenon in QNNs with periodic structure. We define overparametrization as the regime where the QNN has more than a critical number of parameters $M_c$ that allows it to explore all relevant directions in state space. Our main results show that the dimension of the Lie algebra obtained from the generators of the QNN is an upper bound for $M_c$, and for the maximal rank that the quantum Fisher information and Hessian matrices can reach. Underparametrized QNNs have spurious local minima in the loss landscape that start disappearing when $M\geq M_c$. Thus, the overparametrization onset corresponds to a computational phase transition where the QNN trainability is greatly improved by a more favorable landscape. We then connect the notion of overparametrization to the QNN capacity, so that when a QNN is overparametrized, its capacity achieves its maximum possible value. We run numerical simulations for eigensolver, compilation, and autoencoding applications to showcase the overparametrization computational phase transition. We note that our results also apply to variational quantum algorithms and quantum optimal control.


title: Security Review of Ethereum Beacon Clients

abstract: The beacon chain is the backbone of the Ethereum's evolution towards a proof-of-stake-based scalable network. Beacon clients are the applications implementing the services required to operate the beacon chain, namely validators, beacon nodes, and slashers. Security defects in beacon clients could lead to loss of funds, consensus rules violation, network congestion, and other inconveniences. We reported more than 35 issues to the beacon client developers, including various security improvements, specification inconsistencies, missing security checks, exposure to known vulnerabilities. None of our findings appears to be high-severity. We covered the four main beacon clients, namely Lighthouse (Rust), Nimbus (Nim), Prysm (Go), and Teku (Java). We looked for bugs in the logic and implementation of the new security-critical components (BLS signatures, slashing, networking protocols, and API) over a 3-month project that followed a preliminary analysis of BLS signatures code. We focused on Lighthouse and Prysm, the most popular clients, and thus the highest-value targets. Furthermore, we identify protocol-level issues, including replay attacks and incomplete forward secrecy. In addition, we reviewed the network fingerprints of beacon clients, discussing the information obtainable from passive and active searches, and we analyzed the supply chain risk related to third-party dependencies, providing indicators and recommendations to reduce the risk of backdoors and unpatchable vulnerabilities. Our results suggest that despite intense scrutiny by security auditors and independent researchers, the complexity and constant evolution of a platform like Ethereum requires regular expert review and thorough SSDLC practices.


title: Improved optimization strategies for deep Multi-Task Networks

abstract: In Multi-Task Learning (MTL), it is a common practice to train multi-task networks by optimizing an objective function, which is a weighted average of the task-specific objective functions. Although the computational advantages of this strategy are clear, the complexity of the resulting loss landscape has not been studied in the literature. Arguably, its optimization may be more difficult than a separate optimization of the constituting task-specific objectives. In this work, we investigate the benefits of such an alternative, by alternating independent gradient descent steps on the different task-specific objective functions and we formulate a novel way to combine this approach with state-of-the-art optimizers. As the separation of task-specific objectives comes at the cost of increased computational time, we propose a random task grouping as a trade-off between better optimization and computational efficiency. Experimental results over three well-known visual MTL datasets show better overall absolute performance on losses and standard metrics compared to an averaged objective function and other state-of-the-art MTL methods. In particular, our method shows the most benefits when dealing with tasks of different nature and it enables a wider exploration of the shared parameter space. We also show that our random grouping strategy allows to trade-off between these benefits and computational efficiency.


title: Safe Policy Learning through Extrapolation: Application to Pre-trial Risk Assessment

abstract: Algorithmic recommendations and decisions have become ubiquitous in today's society. Many of these and other data-driven policies, especially in the realm of public policy, are based on known, deterministic rules to ensure their transparency and interpretability. For example, algorithmic pre-trial risk assessments, which serve as our motivating application, provide relatively simple, deterministic classification scores and recommendations to help judges make release decisions. How can we use the data based on existing deterministic policies to learn new and better policies? Unfortunately, prior methods for policy learning are not applicable because they require existing policies to be stochastic rather than deterministic. We develop a robust optimization approach that partially identifies the expected utility of a policy, and then finds an optimal policy by minimizing the worst-case regret. The resulting policy is conservative but has a statistical safety guarantee, allowing the policy-maker to limit the probability of producing a worse outcome than the existing policy. We extend this approach to common and important settings where humans make decisions with the aid of algorithmic recommendations. Lastly, we apply the proposed methodology to a unique field experiment on pre-trial risk assessment instruments. We derive new classification and recommendation rules that retain the transparency and interpretability of the existing instrument while potentially leading to better overall outcomes at a lower cost.


title: Simple and Effective Zero-shot Cross-lingual Phoneme Recognition

abstract: Recent progress in self-training, self-supervised pretraining and unsupervised learning enabled well performing speech recognition systems without any labeled data. However, in many cases there is labeled data available for related languages which is not utilized by these methods. This paper extends previous work on zero-shot cross-lingual transfer learning by fine-tuning a multilingually pretrained wav2vec 2.0 model to transcribe unseen languages. This is done by mapping phonemes of the training languages to the target language using articulatory features. Experiments show that this simple method significantly outperforms prior work which introduced task-specific architectures and used only part of a monolingually pretrained model.


title: On the viability of the evolution of the universe with Geometric Inflation

abstract: We perform a general analysis of the cosmological viability of Geometric Inflation. We show that the evolution of the universe, from inflation to the present day, can be seen from the addition of an infinite tower of curvature invariants into the Hilbert-Einstein action. The main epochs of the Universe can be reproduced: Inflation, Big Bang Nucleosynthesis, and Late-time acceleration driven by the cosmological constant. The slow-roll condition is a robust prediction of the theory. Inflation possesses a graceful exit with enough number of $e-$folds between the limit imposed by the Planck density and the exit of the exponential expansion to solve the horizon problem and the absence of topological defects. We also provide some scenarios where the energy scale of the theory can be calibrated.


title: Paint4Poem: A Dataset for Artistic Visualization of Classical Chinese Poems

abstract: In this work we propose a new task: artistic visualization of classical Chinese poems, where the goal is to generatepaintings of a certain artistic style for classical Chinese poems. For this purpose, we construct a new dataset called Paint4Poem. Thefirst part of Paint4Poem consists of 301 high-quality poem-painting pairs collected manually from an influential modern Chinese artistFeng Zikai. As its small scale poses challenges for effectively training poem-to-painting generation models, we introduce the secondpart of Paint4Poem, which consists of 3,648 caption-painting pairs collected manually from Feng Zikai's paintings and 89,204 poem-painting pairs collected automatically from the web. We expect the former to help learning the artist painting style as it containshis most paintings, and the latter to help learning the semantic relevance between poems and paintings. Further, we analyze Paint4Poem regarding poem diversity, painting style, and the semantic relevance between poems and paintings. We create abenchmark for Paint4Poem: we train two representative text-to-image generation models: AttnGAN and MirrorGAN, and evaluate theirperformance regarding painting pictorial quality, painting stylistic relevance, and semantic relevance between poems and paintings.The results indicate that the models are able to generate paintings that have good pictorial quality and mimic Feng Zikai's style, but thereflection of poem semantics is limited. The dataset also poses many interesting research directions on this task, including transferlearning, few-shot learning, text-to-image generation for low-resource data etc. The dataset is publicly available.(https://github.com/paint4poem/paint4poem)


title: Optimal Decision Making in High-Throughput Virtual Screening Pipelines

abstract: The need for efficient computational screening of molecular candidates that possess desired properties frequently arises in various scientific and engineering problems, including drug discovery and materials design. However, the large size of the search space containing the candidates and the substantial computational cost of high-fidelity property prediction models makes screening practically challenging. In this work, we propose a general framework for constructing and optimizing a virtual screening (HTVS) pipeline that consists of multi-fidelity models. The central idea is to optimally allocate the computational resources to models with varying costs and accuracy to optimize the return-on-computational-investment (ROCI). Based on both simulated as well as real data, we demonstrate that the proposed optimal HTVS framework can significantly accelerate screening virtually without any degradation in terms of accuracy. Furthermore, it enables an adaptive operational strategy for HTVS, where one can trade accuracy for efficiency.


title: From the AKNS system to the matrix Schroedinger equation with vanishing potentials: Direct and inverse problems

abstract: We relate the scattering theory of the focusing AKNS system with vanishing boundary conditions to that of the matrix Schroedinger equation. The corresponding Miura transformation which allows this connection, converts the focusing matrix nonlinear Schroedinger (NLS) equation into a new nonlocal integrable equation. We apply the matrix triplet method to derive the multisoliton solutions of the nonlocal integrable equation, thus proposing a new method to solve the matrix NLS equation.


title: From data to reduced-order models via generalized balanced truncation

abstract: This paper proposes a data-driven model reduction approach on the basis of noisy data. Firstly, the concept of data reduction is introduced. In particular, we show that the set of reduced-order models obtained by applying a Petrov-Galerkin projection to all systems explaining the data characterized in a large-dimensional quadratic matrix inequality (QMI) can again be characterized in a lower-dimensional QMI. Next, we develop a data-driven generalized balanced truncation method that relies on two steps. First, we provide necessary and sufficient conditions such that systems explaining the data have common generalized Gramians. Second, these common generalized Gramians are used to construct projection matrices that allow to characterize a class of reduced-order models via generalized balanced truncation in terms of a lower-dimensional QMI by applying the data reductionconcept. Additionally, we present alternative procedures to compute a priori and a posteriori upper bounds with respect to the true system generating the data. Finally, the proposed techniques are illustrated by means of application to an example of a system of a cart with a double-pendulum.


title: Revealing pion and kaon structure via generalised parton distributions

abstract: Clear windows onto emergent hadron mass (EHM) and modulations thereof by Higgs boson interactions are provided by observable measures of pion and kaon structure, many of which are accessible via generalised parton distributions (GPDs). Beginning with algebraic GPD Ansaetze, constrained entirely by hadron-scale $\pi$ and $K$ valence-parton distribution functions (DFs), in whose forms both EHM and Higgs boson influences are manifest, numerous illustrations are provided. They include the properties of electromagnetic form factors, impact parameter space GPDs, gravitational form factors and associated pressure profiles, and the character and consequences of all-orders evolution. The analyses predict that mass-squared gravitational form factors are stiffer than electromagnetic form factors; reveal that $K$ pressure profiles are tighter than $\pi$ profiles, with both mesons sustaining near-core pressures at magnitudes similar to that expected at the core of neutron stars; deliver parameter-free predictions for $\pi$ and $K$ valence, glue, and sea GPDs at the resolving scale $\zeta=2$GeV; and predict that at this scale the fraction of meson mass-squared carried by glue and sea combined matches that lodged with the valence degrees-of-freedom, with a similar statement holding for mass-squared radii.


title: An effective field theory approach for the $S = 3/2$ bilayer honeycomb antiferromagnet

abstract: The spin-3/2 Heisenberg antiferromagnet on the bilayer honeycomb lattice is a minimal model to describe the magnetic behavior of Bi$_3$Mn$_4$O$_{12}$(NO$_3$). We study this model with frustrating inter-layer second-neighbor couplings, taking into account quantum and thermal fluctuations. We use a path integral formulation in terms of coherent states to describe the low energy physics of the model. We show that for a particular point in the parameter space, close to the experimental estimated couplings, a continuum classical degeneracy is lifted by both quantum and thermal fluctuations, and a collinear state is then selected by an order by disorder mechanism. Our results provide a global perspective in the understanding of the experimental observations.


title: Entropy scaling law and the quantum marginal problem: simplification and generalization

abstract: Recently, we introduced a solution to the quantum marginal problem relevant to two-dimensional quantum many-body systems [I. H. Kim, Phys. Rev. X, 11, 021039]. One of the conditions was that the marginals are internally translationally invariant. We show that this condition can be replaced by a weaker condition, namely the local consistency of the marginals. This extends the applicability of the solution to any quantum many-body states in two dimensions that satisfy the entropy scaling law, with or without symmetry. We also significantly simplify the proof by advocating the usage of the maximum-entropy principle.


title: Photometry and model of near-Earth asteroid 2021 DW1 from one apparition

abstract: On 4 March 2021 at 9 UTC a 30-m in diameter near-Earth asteroid 2021 DW1 passed the Earth at a distance of 570000 km, reaching the maximum brightness of V=14.6 mag. We observed it photometrically from 2 March, when it was visible at V=16.5 mag, until 7 March (V=18.2 mag). During that time 2021 DW1 swept a 170 degrees long arc in the northern sky, spanning solar phase angles in the range from 36 to 86 degrees. This made it an excellent target for physical characterisation, including spin axis and shape derivation. Convex inversion of the asteroid lightcurves gives a sidereal period of rotation P=0.013760 +/- 0.000001 h, and two solutions for the spin axis ecliptic coordinates: (A) lambda_1=57 +/- 10, beta_1=29 +/- 10, and (B) lambda_2=67 +/- 10, beta_2=-40 +/- 10. The magnitude-phase curve can be fitted with a standard H, G function with H=24.8 +/- 0.5 mag and an assumed G=0.24. The asteroid colour indices are g-i=0.79 +/- 0.01 mag, and i-z=0.01 +/- 0.02 mag which indicates an S taxonomic class, with an average geometric albedo p_V=0.23 +/- 0.02. The asteroid effective diameter, derived from H and p_V, is D=30 +/- 10 m. It was found that the inclination of the spin axis of 2021 DW1 is not perpendicular to the orbital plane (obliquity epsilon=54 +/- 10 or epsilon=123 +/- 10). More spin axes of VSAs should be determined to check, if 2021 DW1 is an exception or a typical case.


title: Discovering and Validating AI Errors With Crowdsourced Failure Reports

abstract: AI systems can fail to learn important behaviors, leading to real-world issues like safety concerns and biases. Discovering these systematic failures often requires significant developer attention, from hypothesizing potential edge cases to collecting evidence and validating patterns. To scale and streamline this process, we introduce crowdsourced failure reports, end-user descriptions of how or why a model failed, and show how developers can use them to detect AI errors. We also design and implement Deblinder, a visual analytics system for synthesizing failure reports that developers can use to discover and validate systematic failures. In semi-structured interviews and think-aloud studies with 10 AI practitioners, we explore the affordances of the Deblinder system and the applicability of failure reports in real-world settings. Lastly, we show how collecting additional data from the groups identified by developers can improve model performance.


title: Four consecutive primitive elements in a finite field

abstract: For $q$ an odd prime power, we prove that there are always four consecutive primitive elements in the finite field $\mathbb{F}_{q}$ when $q> 2401$.


title: Dimension-Free Rates for Natural Policy Gradient in Multi-Agent Reinforcement Learning

abstract: Cooperative multi-agent reinforcement learning is a decentralized paradigm in sequential decision making where agents distributed over a network iteratively collaborate with neighbors to maximize global (network-wide) notions of rewards. Exact computations typically involve a complexity that scales exponentially with the number of agents. To address this curse of dimensionality, we design a scalable algorithm based on the Natural Policy Gradient framework that uses local information and only requires agents to communicate with neighbors within a certain range. Under standard assumptions on the spatial decay of correlations for the transition dynamics of the underlying Markov process and the localized learning policy, we show that our algorithm converges to the globally optimal policy with a dimension-free statistical and computational complexity, incurring a localization error that does not depend on the number of agents and converges to zero exponentially fast as a function of the range of communication.

categories: cs.LG cs.MA cs.SY eess.SY math.ST stat.ML stat.TH


title: Updating the Theory of Buffer Sizing

abstract: Routers have packet buffers to reduce packet drops during times of congestion. It is important to correctly size the buffer: make it too small, and packets are dropped unnecessarily and the link may be underutilized; make it too big, and packets may wait for a long time, and the router itself may be more expensive to build. Despite its importance, there are few guidelines for picking the buffer size. The two most well-known rules only apply to long-lived TCP Reno flows; either for a network carrying a single TCP Reno flow (the buffer size should equal the bandwidth-delay product, or $BDP$) or for a network carrying $n$ TCP Reno flows (the buffer size should equal $BDP/\sqrt{n}$). Since these rules were introduced, TCP Reno has been replaced by newer algorithms as the default congestion control algorithm in all major operating systems, yet little has been written about how the rules need to change. This paper revisits both rules. For the single flow case, we generalize the $BDP$ rule to account for changes to TCP, such as Proportional Rate Reduction (PRR), and the introduction of new algorithms including Cubic and BBR. We find that buffers can be made 60-75% smaller for newer algorithms. For the multiple flow case, we show that the square root of $n$ rule holds under a broader set of assumptions than previously known, including for these new congestion control algorithms. We also demonstrate situations where the square root of $n$ rule does not hold, including for unfair flows and certain settings with ECN. We validate our results by precisely measuring the time series of buffer occupancy in a real network, and comparing it to the per-packet window size.


title: Component-by-component construction of randomized rank-1 lattice rules achieving almost the optimal randomized error rate

abstract: We study a randomized quadrature algorithm to approximate the integral of periodic functions defined over the high-dimensional unit cube. Recent work by Kritzer, Kuo, Nuyens and Ullrich (2019) shows that rank-1 lattice rules with a randomly chosen number of points and good generating vector achieve almost the optimal order of the randomized error in weighted Korobov spaces, and moreover, that the error is bounded independently of the dimension if the weight parameters, $\gamma_j$, satisfy the summability condition $\sum_{j=1}^{\infty}\gamma_j^{1/\alpha}<\infty$, where $\alpha$ is a smoothness parameter. The argument is based on the existence result that at least half of the possible generating vectors yield almost the optimal order of the worst-case error in the same function spaces. In this paper we provide a component-by-component construction algorithm of such randomized rank-1 lattice rules, without any need to check whether the constructed generating vectors satisfy a desired worst-case error bound. Similarly to the above-mentioned work, we prove that our algorithm achieves almost the optimal order of the randomized error and that the error bound is independent of the dimension if the same condition $\sum_{j=1}^{\infty}\gamma_j^{1/\alpha}<\infty$ holds. We also provide analogous results for tent-transformed lattice rules for weighted half-period cosine spaces and for polynomial lattice rules in weighted Walsh spaces, respectively.


title: Generalized transitionless quantum driving for open quantum systems

abstract: A general approach for transitionless quantum driving in open quantum systems is introduced. Under the assumption of adiabatic evolution for time-local master equations, we derive the generalized transitionless Lindbladian required to implement a shortcut to adiabaticity in an open system scenario. The general counter-diabatic Lindbladian obtained accounts for a phase freedom, which translates into a set of free parameters throughout the dynamics. We then discuss how our generalized approach allows us to recover the transitionless Lindbladian introduced by G. Vacanti et al. [New J. Phys. 16, 053017 (2014)]. We then show how to engineer time-independent master equations that provide the same dynamics as the time-dependent master equation provided by the standard transitionless quantum driving in open systems. We illustrate our results by applying them both to the adiabatic Deutsch algorithm under dephasing and to the Landau-Zener Hamiltonian under bit-phase-flip.


title: Observational Evidence of Magnetic Reconnection in the Terrestrial Foreshock Region

abstract: Electron heating/acceleration in the foreshock, by which electrons may be energized beyond thermal energies prior to encountering the bow shock, is very important for the bow shock dynamics. And then these electrons would be more easily injected into a process like diffusive shock acceleration. Many mechanisms have been proposed to explain electrons heating/acceleration in the foreshock. Magnetic reconnection is one possible candidate. Taking advantage of the Magnetospheric Multiscale mission, we present two magnetic reconnection events in the dawn-side and dusk-side ion foreshock region, respectively. Super-Alfv\'enic electron outflow, demagnetization of the electrons and the ions, and crescent electron distributions in the plane perpendicular to the magnetic field are observed in the sub-ion-scale current sheets. Moreover, strong energy conversion from the fields to the plasmas and significant electron temperature enhancement are observed. Our observations provide direct evidence that magnetic reconnection could occur in the foreshock region and heat/accelerate the electrons therein.


title: In Situ Detection of Kinetic-Size Magnetic Holes in the Martian Magnetosheath

abstract: Depression in magnetic field strength with a scale below one proton gyroradius is referred to as kinetic-size magnetic hole (KSMH). KSMHs are frequently observed near terrestrial space environments and are thought to play an important role in electron energization and energy dissipation in space plasmas. Recently, KSMHs have been evidenced in the Venusian magnetosheath. However, observations of KSMHs in other planetary environments are still lacking. In this study, we present the in situ detection of KSMHs in Martian magnetosheath using Mars Atmosphere and Volatile EvolutioN (MAVEN) for the first time. The distribution of KSMHs is asymmetry in the southern northern hemisphere and no obvious asymmetry in the dawn dusk hemisphere. The observed KSMHs are accompanied by increases in the electron fluxes in the perpendicular direction, indicating the cues of trapped electrons and the formation of electron vortices inside KSMHs. These features are similar to the observations in the terrestrail magtosheath and magnetotail plasma sheet and the Venusian magnetosheath. This implies that KSMHs are a universal magnetic structure in space.


title: Saturation mechanism of the fluctuation dynamo in supersonic turbulent plasmas

abstract: Magnetic fields in several astrophysical objects are amplified and maintained by a dynamo mechanism, which is the conversion of the turbulent kinetic energy to magnetic energy. A dynamo that amplifies magnetic fields at scales $<$ the driving scale of turbulence is known as the fluctuation dynamo. We study the properties of the fluctuation dynamo in supersonic turbulent plasmas, which is of relevance to the ISM, structure formation, and lab experiments of laser-plasma turbulence. Using simulations, we explore the properties of the exponentially growing and saturated state of the fluctuation dynamo for subsonic and supersonic turbulence. We confirm that the fluctuation dynamo efficiency decreases with compressibility. We show that the fluctuation dynamo generated magnetic fields are spatially intermittent and the level of intermittency decreases as the field saturates. We find a stronger back reaction of the magnetic field on the velocity for the subsonic case as compared to the supersonic case. Locally, we find that the level of alignment between vorticity and velocity, velocity and magnetic field, and current density and magnetic field in the saturated stage is enhanced in comparison to the exponentially growing phase for the subsonic case, but only the current density and magnetic field alignment is enhanced for the supersonic case. We show that both the magnetic field amplification (due to weaker stretching of field lines) and diffusion decreases when the field saturates, but the diffusion is enhanced relative to amplification. This occurs throughout the volume in the subsonic turbulence, but primarily in the strong-field regions for the supersonic case. This leads to the saturation of the fluctuation dynamo. Overall, both the amplification and diffusion of magnetic fields are affected and thus a drastic change in either of them is not required for the saturation. [Abstract abridged]


title: Broadband complementary vibrational spectroscopy with cascaded intra-pulse difference frequency generation

abstract: One of the essential goals of molecular spectroscopy is to measure all fundamental molecular vibrations simultaneously. To this end, one needs to measure broadband infrared (IR) absorption and Raman scattering spectra, which provide complementary vibrational information. A recently demonstrated technique called complementary vibrational spectroscopy (CVS) enables simultaneous measurements of IR and Raman spectra with a single device based on a single laser source. However, the spectral coverage was limited to ~1000 cm-1, which partially covers the spectral regions of the fundamental vibrations. In this work, we demonstrate a simple method to expand the spectral bandwidth of the CVS with a cascaded intra-pulse difference-frequency generation (IDFG). Using the system, we measure broadband CVS spectra of organic liquids spanning over 2000 cm-1, more than double the previous study.


title: Untrained Graph Neural Networks for Denoising

abstract: A fundamental problem in signal processing is to denoise a signal. While there are many well-performing methods for denoising signals defined on regular supports, such as images defined on two-dimensional grids of pixels, many important classes of signals are defined over irregular domains such as graphs. This paper introduces two untrained graph neural network architectures for graph signal denoising, provides theoretical guarantees for their denoising capabilities in a simple setup, and numerically validates the theoretical results in more general scenarios. The two architectures differ on how they incorporate the information encoded in the graph, with one relying on graph convolutions and the other employing graph upsampling operators based on hierarchical clustering. Each architecture implements a different prior over the targeted signals. To numerically illustrate the validity of the theoretical results and to compare the performance of the proposed architectures with other denoising alternatives, we present several experimental results with real and synthetic datasets.


title: Microphysical manifestations of viscosity and consequences for anisotropies in the very early universe

abstract: It has been known that a non-perfect fluid that accounts for dissipative viscous effects can evade a highly anisotropic chaotic mixmaster approach to a singularity. Viscosity is often simply parameterised in this context, so it remains unclear whether isotropisation can really occur in physically motivated contexts. We present a few examples of microphysical manifestations of viscosity in fluids that interact either gravitationally or, for a scalar field for instance, through a self-coupling term in the potential. In each case, we derive the viscosity coefficient and comment on the applicability of the approximations involved when dealing with dissipative non-perfect fluids. Upon embedding the fluids in a cosmological context, we then show the extent to which these models allow for isotropisation of the universe in the approach to a singularity. We first do this in the context of expansion anisotropy only, i.e., in the case of a Bianchi type-I universe. We then include anisotropic 3-curvature modelled by the Bianchi type-IX metric. It is found that a self-interacting scalar field at finite temperature allows for efficient isotropisation, whether in a Bianchi type-I or type-IX spacetime, although the model is not tractable all the way to a singularity. Mixmaster chaotic behaviour, which is well known to arise in anisotropic models including anisotropic 3-curvature, is found to be suppressed in the latter case as well. We find that the only model permitting an isotropic singularity is that of a dense gas of black holes.


title: Stable representation theory: beyond the classical groups

abstract: The orthogonal groups are a series of simple Lie groups associated to symmetric bilinear forms. There is no analogous series associated to symmetric trilinear forms. We introduce an infinite dimensional group-like object that can be viewed as the limit of this non-existent series, were it to exist. We show that the representation theory of this object is well-behaved, and similar to the stable representation theory of orthogonal groups. Our theory is not specific to symmetric trilinear forms, and applies to any kind of tensorial forms. Our results can be also be viewed from the perspective of semi-linear representations of the infinite general linear group, and are closely related to twisted commutative algebras.


title: An Improved Frequent Directions Algorithm for Low-Rank Approximation via Block Krylov Iteration

abstract: Frequent Directions, as a deterministic matrix sketching technique, has been proposed for tackling low-rank approximation problems. This method has a high degree of accuracy and practicality, but experiences a lot of computational cost for large-scale data. Several recent works on the randomized version of Frequent Directions greatly improve the computational efficiency, but unfortunately sacrifice some precision. To remedy such issue, this paper aims to find a more accurate projection subspace to further improve the efficiency and effectiveness of the existing Frequent Directions techniques. Specifically, by utilizing the power of Block Krylov Iteration and random projection technique, this paper presents a fast and accurate Frequent Directions algorithm named as r-BKIFD. The rigorous theoretical analysis shows that the proposed r-BKIFD has a comparable error bound with original Frequent Directions, and the approximation error can be arbitrarily small when the number of iterations is chosen appropriately. Extensive experimental results on both synthetic and real data further demonstrate the superiority of r-BKIFD over several popular Frequent Directions algorithms, both in terms of computational efficiency and accuracy.


title: A Parallel Tempering Approach for Efficient Exploration of the Verification Tradespace in Engineered Systems

abstract: Verification is a critical process in the development of engineered systems. Through verification, engineers gain confidence in the correct functionality of the system before it is deployed into operation. Traditionally, verification strategies are fixed at the beginning of the system's development and verification activities are executed as the development progresses. Such an approach appears to give inferior results as the selection of the verification activities does not leverage information gained through the system's development process. In contrast, a set-based design approach to verification, where verification activities are dynamically selected as the system's development progresses, has been shown to provide superior results. However, its application under realistic engineering scenarios remains unproven due to the large size of the verification tradespace. In this work, we propose a parallel tempering approach (PTA) to efficiently explore the verification tradespace. First, we formulate exploration of the verification tradespace as a tree search problem. Second, we design a parallel tempering (PT) algorithm by simulating several replicas of the verification process at different temperatures to obtain a near-optimal result. Third, We apply the PT algorithm to all possible verification states to dynamically identify near-optimal results. The effectiveness of the proposed PTA is evaluated on a partial model of a notional satellite optical instrument.


title: Dimension-Grouped Mixed Membership Models for Multivariate Categorical Data

abstract: Mixed Membership Models (MMMs) are a popular family of latent structure models for complex multivariate data. Instead of forcing each subject to belong to a single cluster, MMMs incorporate a vector of subject-specific weights characterizing partial membership across clusters. With this flexibility come challenges in uniquely identifying, estimating, and interpreting the parameters. In this article, we propose a new class of Dimension-Grouped MMMs (Gro-M$^3$s) for multivariate categorical data, which improve parsimony and interpretability. In Gro-M$^3$s, observed variables are partitioned into groups such that the latent membership is constant for variables within a group but can differ across groups. Traditional latent class models are obtained when all variables are in one group, while traditional MMMs are obtained when each variable is in its own group. The new model corresponds to a novel decomposition of probability tensors. Theoretically, we derive transparent identifiability conditions for both the unknown grouping structure and model parameters in general settings. Methodologically, we propose a Bayesian approach for Dirichlet Gro-M$^3$s to inferring the variable grouping structure and estimating model parameters. Simulation results demonstrate good computational performance and empirically confirm the identifiability results. We illustrate the new methodology through applications to a functional disability survey dataset and a personality test dataset.


title: Indoor Navigation Algorithm Based on a Smartphone Inertial Measurement Unit and Map Matching

abstract: We propose an indoor navigation algorithm based on pedestrian dead reckoning (PDR) using an inertial measurement unit in a smartphone and map matching. The proposed indoor navigation system is user-friendly and convenient because it requires no additional device except a smartphone and works with a pedestrian in a casual posture who is walking with a smartphone in their hand. Because the performance of the PDR decreases over time, we greatly reduced the position error of the trajectory estimated by PDR using a map matching method with a known indoor map. To verify the proposed indoor navigation algorithm, we conducted an experiment in a real indoor environment using a commercial Android smartphone. The performance of our algorithm was demonstrated through the results of the experiment.


title: A Decomposition Augmented Lagrangian Method for Low-rank Semidefinite Programming

abstract: We develop a decomposition method based on the augmented Lagrangian framework to solve a broad family of semidefinite programming problems, possibly with nonlinear objective functions, nonsmooth regularization, and general linear equality/inequality constraints. In particular, the positive semidefinite variable along with a group of linear constraints can be transformed into a variable on a smooth manifold via matrix factorization. The nonsmooth regularization and other general linear constraints are handled by the augmented Lagrangian method. Therefore, each subproblem can be solved by a semismooth Newton method on a manifold. Theoretically, we show that the first and second-order necessary optimality conditions for the factorized subproblem are also sufficient for the original subproblem under certain conditions. Convergence analysis is established for the Riemannian subproblem and the augmented Lagrangian method. Extensive numerical experiments on large-scale semidefinite programming problems such as max-cut, nearest correlation estimation, clustering, and sparse principal component analysis demonstrate the strength of our proposed method compared to other state-of-the-art methods.


title: Detect and Perturb: Neutral Rewriting of Biased and Sensitive Text via Gradient-based Decoding

abstract: Written language carries explicit and implicit biases that can distract from meaningful signals. For example, letters of reference may describe male and female candidates differently, or their writing style may indirectly reveal demographic characteristics. At best, such biases distract from the meaningful content of the text; at worst they can lead to unfair outcomes. We investigate the challenge of re-generating input sentences to 'neutralize' sensitive attributes while maintaining the semantic meaning of the original text (e.g. is the candidate qualified?). We propose a gradient-based rewriting framework, Detect and Perturb to Neutralize (DEPEN), that first detects sensitive components and masks them for regeneration, then perturbs the generation model at decoding time under a neutralizing constraint that pushes the (predicted) distribution of sensitive attributes towards a uniform distribution. Our experiments in two different scenarios show that DEPEN can regenerate fluent alternatives that are neutral in the sensitive attribute while maintaining the semantics of other attributes.


title: User-Defined Functions for HDF5

abstract: Scientific datasets are known for their challenging storage demands and the associated processing pipelines that transform their information. Some of those processing tasks include filtering, cleansing, aggregation, normalization, and data format translation -- all of which generate even more data. In this paper, we present an infrastructure for the HDF5 file format that enables dataset values to be populated on the fly: task-related scripts can be attached into HDF5 files and only execute when the dataset is read by an application. We provide details on the software architecture that supports user-defined functions (UDFs) and how it integrates with hardware accelerators and computational storage. Moreover, we describe the built-in security model that limits the system resources a UDF can access. Last, we present several use cases that show how UDFs can be used to extend scientific datasets in ways that go beyond the original scope of this work.


title: Gauge invariant input to neural network for path optimization method

abstract: We investigate the efficiency of a gauge invariant input to a neural network for the path optimization method. While the path optimization with a completely gauge-fixed link-variable input has successfully tamed the sign problem in a simple gauge theory, the optimization does not work well when the gauge degrees of freedom remain. We propose to employ a gauge invariant input, such as plaquette, to overcome this problem. The efficiency of the gauge invariant input to the neural network is evaluated for the 2-dimensional $U(1)$ gauge theory with a complex coupling. The average phase factor is significantly enhanced by the path optimization with the plaquette input, indicating good control of the sign problem. It opens a possibility that the path optimization is available to complicated gauge theories, including Quantum Chromodynamics, in a realistic setup.


title: Stable Volumes for Persistent Homology

abstract: This paper proposes a stable volume and a stable volume variant, referred to as a stable sub-volume, for more reliable data analysis using persistent homology. In prior research, an optimal cycle and similar ideas have been proposed to identify the homological structure corresponding to each birth-death pair in a persistence diagram. While this is helpful for data analysis using persistent homology, the results are sensitive to noise. In this paper, stable volumes and stable sub-volumes are proposed to solve this problem. For a special case, we prove that a stable volume is the robust part of an optimal volume against noise. We implemented stable volumes and sub-volumes on HomCloud, a data analysis software package based on persistent homology, and show examples of stable volumes and sub-volumes.


title: Feasibility study of urban flood mapping using traffic signs for route optimization

abstract: Water events are the most frequent and costliest climate disasters around the world. In the U.S., an estimated 127 million people who live in coastal areas are at risk of substantial home damage from hurricanes or flooding. In flood emergency management, timely and effective spatial decision-making and intelligent routing depend on flood depth information at a fine spatiotemporal scale. In this paper, crowdsourcing is utilized to collect photos of submerged stop signs, and pair each photo with a pre-flood photo taken at the same location. Each photo pair is then analyzed using deep neural network and image processing to estimate the depth of floodwater in the location of the photo. Generated point-by-point depth data is converted to a flood inundation map and used by an A* search algorithm to determine an optimal flood-free path connecting points of interest. Results provide crucial information to rescue teams and evacuees by enabling effective wayfinding during flooding events.


title: Design of arbitrarily shaped acoustic cloaks through PDE-constrained optimization satisfying sonic-metamaterial design requirements

abstract: We develop an optimization framework for the design of acoustic cloaks, with the aim of overcoming the limitations of usual transformation-based cloaks in terms of microstructure complexity and shape arbitrarity of the obstacle. This is achieved by recasting the acoustic cloaking design as a nonlinear optimal control problem constrained by a linear elliptic partial differential equation. In this setting, isotropic material properties' distributions realizing the cloak take the form of control functions and a system of first-order optimality conditions is derived accordingly. Such isotropic media can then be obtained in practice with simple hexagonal lattices of inclusions in water. For this reason, the optimization problem is directly formulated to take into account suitable partitions of the control domain Two types of inclusions are considered, and long-wavelength homogenization is used to define the feasible set of material properties that is employed as a constraint in the optimization problem. In this manner, we link the stage of material properties optimization with that of microstructure design, aiming at finding the optimal implementable solution. As a test benchmark, cloaking of the silhouette of a ship is considered, for various frequencies and directions of incidence. The resulting cloak is numerically tested via coupled structural/acoustic simulations.


title: Numerical investigation of non-Gaussianities in the phase and modulus of density Fourier modes

abstract: We numerically investigate non-Gaussianities in the late-time cosmological density field in Fourier space. We explore various statistics, including the two-point and three-point probability distribution function (PDF) of phase and modulus, and two \& three-point correlation function of of phase and modulus. We detect significant non-Gaussianity for certain configurations. We compare the simulation results with the theoretical expansion series of \citet{2007ApJS..170....1M}. We find that the $\mathcal{O}(V^{-1/2})$ order term alone is sufficiently accurate to describe all the measured non-Gaussianities in not only the PDFs, but also the correlations. We also numerically find that the phase-modulus cross-correlation contributes $\sim 50\%$ to the bispectrum, further verifying the accuracy of the $\mathcal{O}(V^{-1/2})$ order prediction. This work demonstrates that non-Gaussianity of the cosmic density field is simpler in Fourier space, and may facilitate the data analysis in the era of precision cosmology.