Search

The inference of causal relationships among observed variables is a pivotal, longstanding problem in the scientific community. An intuitive method for quantifying these causal links involves examining the response of one variable to perturbations in another. The fluctuation-dissipation theorem elegantly connects this response to the correlation functions of the unperturbed system, thereby bridging the concepts of causality and correlation. However, this relationship becomes intricate in nonlinear systems, where knowledge of the invariant measure is required but elusive, especially in high-dimensional spaces. In this study, we establish a novel link between the Koopman operator of nonlinear stochastic systems and the response function. This connection provides an alternative method for computing the response function using generalized correlation functions, even when the invariant measure is unknown. We validate our theoretical framework by applying it to a nonlinear high-dimensional system amenable to exact solutions, demonstrating convergence and consistency with established results. Finally, we discuss a significant interplay between the resulting causal network and the relevant time scales of the system.

In their comment (1), Hawreliak et al. claims that our observation of stacking fault formation and transonic dislocation propagation in diamond (2) is not valid as they interpret the observed features as cracks. In this response letter, we describe our rationale for interpreting the observed features as stacking faults. We also address other points raised in their comments, including the clarifications of how the results of Makarov et al. (3) are not in conflict with our study., Comment: Comment on arXiv:2401.04213. 10 pages, 2 figures

We present measurements on Fe$_2$O$_3$ amorphization and melt under laser-driven shock compression up to 209(10) GPa via time-resolved in situ x-ray diffraction. At 122(3) GPa, a diffuse signal is observed indicating the presence of a non-crystalline phase. Structure factors have been extracted up to 182(6) GPa showing the presence of two well-defined peaks. A rapid change in the intensity ratio of the two peaks is identified between 145(10) and 151(10) GPa, indicative of a phase change. Present DFT+$U$ calculations of temperatures along Fe$_2$O$_3$ Hugoniot are in agreement with SESAME 7440 and indicate relatively low temperatures, below 2000 K, up to 150 GPa. The non-crystalline diffuse scattering is thus consistent with the - as yet unreported - shock amorphization of Fe$_2$O$_3$ between 122(3) and 145(10) GPa, followed by an amorphous-to-liquid transition above 151(10) GPa. Upon release, a non-crystalline phase is observed alongside crystalline $\alpha$-Fe$_2$O$_3$. The extracted structure factor and pair distribution function of this release phase resemble those reported for Fe$_2$O$_3$ melt at ambient pressure., Comment: 11 pages, 4 figures, under review

Congestion games are attractive because they can model many concrete situations where some competing entities interact through the use of some shared resources, and also because they always admit pure Nash equilibria which correspond to the local minima of a potential function. We explore the problem of computing a state of minimum potential in this setting. Using the maximum number of resources that a player can use at a time, and the possible symmetry in the players' strategy spaces, we settle the complexity of the problem for instances having monotone (i.e., either non-decreasing or non-increasing) latency functions on their resources. The picture, delineating polynomial and NP-hard cases, is complemented with tight approximation algorithms.

We report here the first equation of state measurements of Fe$_2$O$_3$ obtained with laser-driven shock compression. The data are in excellent agreement with previous dynamic and static compression measurements at low pressure, and extend the known Hugoniot up to 700 GPa. We observe a large volume drop of $\sim$10% at 86 GPa, which could be associated, according to static compression observations, with the iron spin transition. Our measurements also suggest a change of the Hugoniot curve between 150 and 250 GPa. Above 250 GPa and within our error bars, we do not observe significant modifications up to the maximum pressure of 700 GPa reached in our experiment., Comment: 8 pages, 3 figures, 2 tables

In this work, we study the Stochastic Budgeted Multi-round Submodular Maximization (SBMSm) problem, where we aim to adaptively maximize the sum, over multiple rounds, of a monotone and submodular objective function defined on subsets of items. The objective function also depends on the realization of stochastic events, and the total number of items we can select over all rounds is bounded by a limited budget. This problem extends, and generalizes to multiple round settings, well-studied problems such as (adaptive) influence maximization and stochastic probing. We show that, if the number of items and stochastic events is somehow bounded, there is a polynomial time dynamic programming algorithm for SBMSm. Then, we provide a simple greedy $1/2(1-1/e-\epsilon)\approx 0.316$-approximation algorithm for SBMSm, that first non-adaptively allocates the budget to be spent at each round, and then greedily and adaptively maximizes the objective function by using the budget assigned at each round. Finally, we introduce the {\em budget-adaptivity gap}, by which we measure how much an adaptive policy for SBMSm is better than an optimal partially adaptive one that, as in our greedy algorithm, determines the budget allocation in advance. We show that the budget-adaptivity gap lies between $e/(e-1)\approx 1.582$ and $2$.

Local networks of neurons are nonlinear systems driven by synaptic currents elicited by its own spiking activity and the input received from other brain areas. Synaptic currents are well approximated by correlated Gaussian noise. Besides, the population dynamics of neuronal networks is often found to be multistable, allowing the noise source to induce state transitions. State changes in neuronal systems underlies the way information is encoded and transformed. The characterization of the escape time from metastable states is then a cornerstone to understand how information is processed in the brain. The effects of correlated input forcing bistable systems have been studied for over half a century, nonetheless most results are perturbative or valid only when a separation of time scales is present. Here, we present a novel and exact result holding when the correlation time of the noise source is identical to that of the neural population, hence solving in a very general setting the mean escape time problem.

Covariance matrix estimation is an important task in the analysis of multivariate data in disparate scientific fields, including neuroscience, genomics, and astronomy. However, modern scientific data are often incomplete due to factors beyond the control of researchers, and data missingness may prohibit the use of traditional covariance estimation methods. Some existing methods address this problem by completing the data matrix, or by filling the missing entries of an incomplete sample covariance matrix by assuming a low-rank structure. We propose a novel approach that exploits auxiliary variables to complete covariance matrix estimates. An example of auxiliary variable is the distance between neurons, which is usually inversely related to the strength of neuronal correlation. Our method extracts auxiliary information via regression, and involves a single tuning parameter that can be selected empirically. We compare our method with other matrix completion approaches via simulations, and apply it to the analysis of large-scale neuroscience data., Comment: 16 pages, 6 figures

Researchers continue exploring neurons' intricate patterns of activity in the cerebral visual cortex in response to visual stimuli. The way neurons communicate and optimize their interactions with each other under different experimental conditions remains a topic of active investigation. Probabilistic Graphical Models are invaluable tools in neuroscience research, as they let us identify the functional connections, or conditional statistical dependencies, between neurons. Graphical models represent these connections as a graph, where nodes represent neurons and edges indicate the presence of functional connections between them. We developed the R package neurofuncon for the computation and visualization of functional connectivity graphs from large-scale data based on the Graphical lasso. We illustrate the use of this package with publicly available two-photon calcium microscopy imaging data from approximately 10000 neurons in a 1mm cubic section of a mouse visual cortex., Comment: 7 pages, 5 figures

Modern Cloud/Edge architectures need to orchestrate multiple layers of heterogeneous computing nodes, including pervasive sensors/actuators, distributed Edge/Fog nodes, centralized data centers and quantum devices. The optimal assignment and scheduling of computation on the different nodes is a very difficult problem, with NP-hard complexity. In this paper, we explore the possibility of solving this problem with Variational Quantum Algorithms, which can become a viable alternative to classical algorithms in the near future. In particular, we compare the performances, in terms of success probability, of two algorithms, i.e., Quantum Approximate Optimization Algorithm (QAOA) and Variational Quantum Eigensolver (VQE). The simulation experiments, performed for a set of simple problems, %CM230124 that involve a Cloud and two Edge nodes, show that the VQE algorithm ensures better performances when it is equipped with appropriate circuit \textit{ansatzes} that are able to restrict the search space. Moreover, experiments executed on real quantum hardware show that the execution time, when increasing the size of the problem, grows much more slowly than the trend obtained with classical computation, which is known to be exponential., Comment: 14 pages, 13 figures

Factor models are widely used in the analysis of high-dimensional data in several fields of research. Estimating a factor model, in particular its covariance matrix, from partially observed data vectors is very challenging. In this work, we show that when the data are structurally incomplete, the factor model likelihood function can be decomposed into the product of the likelihood functions of multiple partial factor models relative to different subsets of data. If these multiple partial factor models are linked together by common parameters, then we can obtain complete maximum likelihood estimates of the factor model parameters and thereby the full covariance matrix. We call this framework Linked Factor Analysis (LINFA). LINFA can be used for covariance matrix completion, dimension reduction, data completion, and graphical dependence structure recovery. We propose an efficient Expectation-Maximization algorithm for maximum likelihood estimation, accelerated by a novel group vertex tessellation (GVT) algorithm which identifies a minimal partition of the vertex set to implement an efficient optimization in the maximization steps. We illustrate our approach in an extensive simulation study and in the analysis of calcium imaging data collected from mouse visual cortex., Comment: 21 page, 9 figures

The persistent issue of wrongful convictions in the United States emphasizes the need for scrutiny and improvement of the criminal justice system. While statistical methods for the evaluation of forensic evidence, including glass, fingerprints, and DNA, have significantly contributed to solving intricate crimes, there is a notable lack of national-level standards to ensure the appropriate application of statistics in forensic investigations. We discuss the obstacles in the application of statistics in court, and emphasize the importance of making statistical interpretation accessible to non-statisticians, especially those who make decisions about potentially innocent individuals. We investigate the use and misuse of statistical methods in crime investigations, in particular the likelihood ratio approach. We further describe the use of graphical models, where hypotheses and evidence can be represented as nodes connected by arrows signifying association or causality. We emphasize the advantages of special graph structures, such as object-oriented Bayesian networks and chain event graphs, which allow for the concurrent examination of evidence of various nature., Comment: 20 pages, 8 figures

The efficient management of energy communities relies on the solution of the "prosumer problem", i.e., the problem of scheduling the household loads on the basis of the user needs, the electricity prices, and the availability of local renewable energy, with the aim of reducing costs and energy waste. Quantum computers can offer a significant breakthrough in treating this problem thanks to the intrinsic parallel nature of quantum operations. The most promising approach is to devise variational hybrid algorithms, in which quantum computation is driven by parameters that are optimized classically, in a cycle that aims at finding the best solution with a significant speed-up with respect to classical approaches. This paper provides a reformulation of the prosumer problem, allowing to address it with a hybrid quantum algorithm, namely, Quantum Approximate Optimization Algorithm (QAOA), and with a recent variant, the Recursive QAOA. We report on an extensive set of experiments, on simulators and real quantum hardware, for different problem sizes. Results are encouraging in that Recursive QAOA is able, for problems involving up to 10 qubits, to provide optimal and admissible solutions with good probabilities, while the computation time is nearly independent of the system size, Comment: 14 pages, 13 figures. IEEE Transactions on Smart Grid (2023)

Data standardization has become one of the leading methods neuroimaging researchers rely on for data sharing and reproducibility. Data standardization promotes a common framework through which researchers can utilize others' data. Yet, as of today, formatting datasets that adhere to community best practices requires technical expertise involving coding and considerable knowledge of file formats and standards. We describe ezBIDS, a tool for converting neuroimaging data and associated metadata to the Brain Imaging Data Structure (BIDS) standard. ezBIDS provides four unique features: (1) No installation or programming requirements. (2) Handling of both imaging and task events data and metadata. (3) Automated inference and guidance for adherence to BIDS. (4) Multiple data management options: download BIDS data to local system, or transfer to OpenNeuro.org or brainlife.io. In sum, ezBIDS requires neither coding proficiency nor knowledge of BIDS and is the first BIDS tool to offer guided standardization, support for task events conversion, and interoperability with OpenNeuro and brainlife.io.

The classification of galaxy morphologies is an important step in the investigation of theories of hierarchical structure formation. While human expert visual classification remains quite effective and accurate, it cannot keep up with the massive influx of data from emerging sky surveys. A variety of approaches have been proposed to classify large numbers of galaxies; these approaches include crowdsourced visual classification, and automated and computational methods, such as machine learning methods based on designed morphology statistics and deep learning. In this work, we develop two novel galaxy morphology statistics, descent average and descent variance, which can be efficiently extracted from telescope galaxy images. We further propose simplified versions of the existing image statistics concentration, asymmetry, and clumpiness, which have been widely used in the literature of galaxy morphologies. We utilize the galaxy image data from the Sloan Digital Sky Survey to demonstrate the effective performance of our proposed image statistics at accurately detecting spiral and elliptical galaxies when used as features of a random forest classifier., Comment: 11 pages, 6 figures

Neuroscience research has expanded dramatically over the past 30 years by advancing standardization and tool development to support rigor and transparency. Consequently, the complexity of the data pipeline has also increased, hindering access to FAIR (Findable, Accessible, Interoperabile, and Reusable) data analysis to portions of the worldwide research community. brainlife.io was developed to reduce these burdens and democratize modern neuroscience research across institutions and career levels. Using community software and hardware infrastructure, the platform provides open-source data standardization, management, visualization, and processing and simplifies the data pipeline. brainlife.io automatically tracks the provenance history of thousands of data objects, supporting simplicity, efficiency, and transparency in neuroscience research. Here brainlife.io's technology and data services are described and evaluated for validity, reliability, reproducibility, replicability, and scientific utility. Using data from 4 modalities and 3,200 participants, we demonstrate that brainlife.io's services produce outputs that adhere to best practices in modern neuroscience research.

Taking informed decisions, namely acting rationally, is an individual attitude of paramount relevance in nature and human societies. In this work, we study how rationality spreads in a community. To this end, through an agent-based model, we analyse the dynamics of a population whose individuals, endowed with a rational attitude controlled by a numerical parameter, play a simple game. The latter consists of multiple strategies, each associated with a given reward. The proposed model is then used as a benchmark for studying the behaviour of Bitcoin users, inferred by analysing transactions recorded in the Blockchain. Remarkably, a population undergoing a sharp transition from irrational to rational attitudes shows a behavioural pattern similar to that of Bitcoin users, whose rationality showed up as soon as their cryptocurrency became worth just a few cents (USD). To conclude, a behavioural analysis that relies on an entropy measure combined with a simple agent-based model allows us to detect the rise of rationality across a community. Although further investigations are essential to corroborate our results, we deem the proposed approach could also get used for studying other social phenomena and behaviours., Comment: 10 pages, 5 figures

Deep generative learning cannot only be used for generating new data with statistical characteristics derived from input data but also for anomaly detection, by separating nominal and anomalous instances based on their reconstruction quality. In this paper, we explore the performance of three unsupervised deep generative models -- variational autoencoders (VAEs) with Gaussian, Bernoulli, and Boltzmann priors -- in detecting anomalies in flight-operations data of commercial flights consisting of multivariate time series. We devised two VAE models with discrete latent variables (DVAEs), one with a factorized Bernoulli prior and one with a restricted Boltzmann machine (RBM) as prior, because of the demand for discrete-variable models in machine-learning applications and because the integration of quantum devices based on two-level quantum systems requires such models. The DVAE with RBM prior, using a relatively simple -- and classically or quantum-mechanically enhanceable -- sampling technique for the evolution of the RBM's negative phase, performed better than the Bernoulli DVAE and on par with the Gaussian model, which has a continuous latent space. Our studies demonstrate the competitiveness of a discrete deep generative model with its Gaussian counterpart on anomaly-detection tasks. Moreover, the DVAE model with RBM prior can be easily integrated with quantum sampling by outsourcing its generative process to measurements of quantum states obtained from a quantum annealer or gate-model device., Comment: 25 pages, 7 figures, 3 tables, appendix, supplementary material

The motion of line defects (dislocations) has been studied for over 60 years but the maximum speed at which they can move is unresolved. Recent models and atomistic simulations predict the existence of a limiting velocity of dislocation motions between the transonic and subsonic ranges at which the self-energy of dislocation diverges, though they do not deny the possibility of the transonic dislocations. We use femtosecond x-ray radiography to track ultrafast dislocation motion in shock-compressed single-crystal diamond. By visualizing stacking faults extending faster than the slowest sound wave speed of diamond, we show the evidence of partial dislocations at their leading edge moving transonically. Understanding the upper limit of dislocation mobility in crystals is essential to accurately model, predict, and control the mechanical properties of materials under extreme conditions.

This paper proposes an optimal management strategy for a Renewable Energy Community defined according to the Italian legislation. The specific case study is composed by a set of houses equipped with smart appliances, that share a PV plant. The objective is to minimize the cost of electrical energy use for each member of the community, taking into account the discount achievable from government incentives with proper shaping of the community daily consumption. Such incentives are indeed proportional to the shared energy, i.e. the portion of the renewable energy consumed at each hour by community members. The management algorithm allows an optimal coordination of houses power demands, according to the degree of flexibility granted by users. Moreover, a policy to fairly distribute the obtained discount is introduced. Simulation results show the potentialities of the approach.

The energy sustainability of blockchains, whose consensus protocol rests on the Proof-of-Work, nourishes a heated debate. The underlying issue lies in a highly energy-consuming process, defined as mining, required to validate crypto-asset transactions. Mining is the process of solving a cryptographic puzzle, incentivised by the possibility of gaining a reward. The higher the number of users performing mining, i.e. miners, the higher the overall electricity consumption of a blockchain. For that reason, mining constitutes a negative environmental externality. Here, we study whether miners' interests can meet the collective need to curb energy consumption. To this end, we introduce the Crypto-Asset Game, namely a model based on the framework of Evolutionary Game Theory devised for studying the dynamics of a population whose agents can play as crypto-asset users or as miners. The energy consumption of mining impacts the payoff of both strategies, representing a direct cost for miners and an environmental factor for crypto-asset users. The proposed model, studied via numerical simulations, shows that, in some conditions, the agent population can reach a strategy profile that optimises global energy consumption, i.e. composed of a low density of miners. To conclude, can a Proof-of-Work-based blockchain become energetically sustainable? Our results suggest that blockchain protocol parameters could have a relevant role in the global energy consumption of this technology., Comment: 7 pages, 5 figures

The behaviour of Bitcoin owners is reflected in the structure and the number of bitcoin transactions encoded in the Blockchain. Likewise, the behaviour of Bitcoin traders is reflected in the formation of bullish and bearish trends in the crypto market. In light of these observations, we wonder if human behaviour underlies some relationship between the Blockchain and the crypto market. To address this question, we map the Blockchain to a spin-lattice problem, whose configurations form ordered and disordered patterns, representing the behaviour of Bitcoin owners. This novel approach allows us to obtain time series suitable to detect a causal relationship between the dynamics of the Blockchain and market trends of the Bitcoin and to find that disordered patterns in the Blockchain precede Bitcoin panic selling. Our results suggest that human behaviour underlying Blockchain evolution and the crypto market brings out a fascinating connection between disorder and panic in Bitcoin dynamics., Comment: 6 pages, 3 figures

Composite materials achieved by including transition-metal oxides with different structures and morphologies in sulfur are suggested as scalable cathodes for high-energy lithium-sulfur (Li-S) batteries. The composites contain 80 wt.% sulfur and 20 wt.% of either MnO2 or TiO2, leading to a sulfur content in the electrode of 64 wt.% and revealing a reversible, fast, and lowly polarized conversion process in the cell with limited interphase resistance. The S-TiO2 composite exhibits an excellent rate capability between C/10 and 2C, and a cycle life extended over 400 cycles at 2C, owing to the effects of the nanometric TiO2 additive in boosting the reaction kinetics. Instead, the micrometric sized particles of MnO2 partially limit the electrochemical activity of S-MnO2 to the current rate of 1C. Nevertheless, both S-MnO2 and S-TiO2 withstand a sulfur loading up to values approaching 6 mgcm-2, and deliver an areal capacit ranging from about 4.5 to 5.5 mAhcm-2 at C/5. The excellent performances of the metal oxide-sulfur electrodes, even at high active material loading, and the possible scalability of the synthetic pathway adopted in the work suggest that the composites are viable cathodes for next-generation Li-S batteries with high energy density and efficient electrochemical process.

An implementation of the electromagnetic Particle in Cell loop in the code Smilei using task programming is presented. Through OpenMP, the macro-particles operations are formulated in terms of tasks. This formulation allows asynchronous execution respecting the data dependencies of the macro-particle operations, the most time-consuming part of the code in simulations of interest for plasma physics. Through some benchmarks it is shown that this formulation can help mitigating the load imbalance of these operations at the OpenMP thread level. The improvements in strong scaling for load-imbalanced physical cases are discussed., Comment: 14 pages, 15 figures, conference: Platform for Advanced Scientific Computing 2022 (PASC22)

A backside-illuminated (BSI) near-infrared enhanced silicon single-photon avalanche diode (SPAD) for time-of-flight (ToF) light detection and ranging applications is presented. The detector contains a 2 $\mu$m wide multiplication region with a spherically-uniform electric field peak enforced by field-line crowding. A charge-focusing electric field extends into a 10 $\mu$m deep absorption volume, whereby electrons generated in all corners of the device can move efficiently towards the multiplication region. The SPAD is integrated with a customized 130 nm CMOS technology and a dedicated BSI process. The device has a pitch of 15 $\mu$m, which has the potential to be scaled down without significant performance loss. Furthermore, the detector achieves a photon detection efficiency of 27% at 905 nm, with an excess bias of 3.5 V that is controlled by integrated CMOS electronics, and a timing resolution of 240 ps. By virtue of these features, the device architecture is well-suited for large format ToF imaging arrays with integrated electronics.

To what extent does the structure of the players' strategy space influence the efficiency of decentralized solutions in congestion games? In this work, we investigate whether better performance are possible when restricting to load balancing games in which players can only choose among single resources. We consider three different solutions concepts, namely, approximate pure Nash equilibria, approximate one-round walks generated by selfish players aiming at minimizing their personal cost and approximate one-round walks generated by cooperative players aiming at minimizing the marginal increase in the sum of the players' personal costs. The last two concepts can be interpreted as solutions of greedy online algorithms for the related resource selection problem. We show that, under fairly general latency functions on the resources, better bounds cannot be achieved if players are either weighted or asymmetric. On the positive side, we prove that, under mild assumptions on the latency functions, improvements on the performance of approximate pure Nash equilibria are possible for load balancing games with weighted and symmetric players in the case of identical resources. We also design lower bounds on the performance of one-round walks in load balancing games with unweighted players and identical resources., Comment: To cite this paper, consider the preliminary version appeared in the proceedings of ESA 2017: "Vittorio Bil\`o, Cosimo Vinci: On the Impact of Singleton Strategies in Congestion Games. The 25th Annual European Symposium on Algorithms (ESA), 17:1-17:14 (2017)"

Despite the huge number of neurons composing a brain network, ongoing activity of local cell assemblies composing cortical columns is intrinsically stochastic. Fluctuations in their instantaneous rate of spike firing $\nu(t)$ scale with the size of the assembly and persist in isolated network, i.e., in absence of external source of noise. Although deterministic chaos due to the quenched disorder of the synaptic couplings likely underlies this seemingly stochastic dynamics, an effective theory for the network dynamics of a finite ensemble of spiking neurons is lacking. Here, we fill this gap by extending the so-called population density approach including an activity- and size-dependent stochastic source in the Fokker-Planck equation for the membrane potential density. The finite-size noise embedded in this stochastic partial derivative equation is analytically characterized leading to a self-consistent and non-perturbative description of $\nu(t)$ valid for a wide class of spiking neuron networks. Its power spectra of $\nu(t)$ are found in excellent agreement with those from detailed simulations both in the linear regime and across a synchronization phase transition, when a size-dependent smearing of the critical dynamics emerges., Comment: 6 pages, 4 figures

Populations of spiking neuron models have densities of their microscopic variables (e.g., single-cell membrane potentials) whose evolution fully capture the collective dynamics of biological networks, even outside equilibrium. Despite its general applicability, the Fokker-Planck equation governing such evolution is mainly studied within the borders of the linear response theory, although alternative spectral expansion approaches offer some advantages in the study of the out-of-equilibrium dynamics. This is mainly due to the difficulty in computing the state-dependent coefficients of the expanded system of differential equations. Here, we address this issue by deriving analytic expressions for such coefficients by pairing perturbative solutions of the Fokker-Planck approach with their counterparts from the spectral expansion. A tight relationship emerges between several of these coefficients and the Laplace transform of the inter-spike interval density (i.e., the distribution of first-passage times). `Coefficients' like the current-to-rate gain function, the eigenvalues of the Fokker-Planck operator and its eigenfunctions at the boundaries are derived without resorting to integral expressions. For the leaky integrate-and-fire neurons, the coupling terms between stationary and nonstationary modes are also worked out paving the way to accurately characterize the critical points and the relaxation timescales in networks of interacting populations., Comment: 11 pages, no figures

Magnetized laser-produced plasmas are central to many new studies in laboratory astrophysics, inertial confinement fusion, and industrial applications. Here we present the results of large-scale, three-dimensional magneto-hydrodynamic simulations of the dynamics of a laser-produced plasma expanding into a transverse magnetic field with a strength of tens of Tesla. The simulations show the plasma is confined by the strong magnetic field into a slender slab structured by the magnetized Rayleigh-Taylor instability that develops at the plasma-vacuum interface. We find that by perturbing the initial velocity of the plume the slab can develop kink-like motion which disrupts its propagation.

Type-III-burst radio signals can be mimicked in the laboratory via laser-plasma interaction. Instead of an electron beam generating Langmuir waves (LW) in the interplanetary medium, the LWs are created by a laser interacting with a millimeter-sized plasma through the stimulated Raman instability. In both cases, the LWs feed the Langmuir decay instability which scatters them in several directions. The resulting LWs may couple to form electromagnetic emission at twice the plasma frequency, which has been detected in the interplanetary medium, and recently in a laboratory laser experiment [Marqu\`es et al. Phys. Rev. Lett. 124, 135001 (2020)]. This article presents the first numerical analysis of this laser configuration using particle-in-cell simulations, providing details on the wave spectra that are too difficult to measure in experiments. The role of some parameters is addressed, with a focus on laser intensity, in order to illustrate the behavior of the electromagnetic emission's angular distribution and polarization.

Ammonia is predicted to be one of the major components in the depths of the ice giant planets Uranus and Neptune. Their dynamics, evolution, and interior structure are insufficiently understood and models rely imperatively on data for equation of state and transport properties. Despite its great significance, the experimentally accessed region of the ammonia phase diagram today is still very limited in pressure and temperature. Here we push the probed regime to unprecedented conditions, up to $\sim$350 GPa and $\sim$40000 K. Along the Hugoniot, the temperature measured as a function of pressure shows a subtle change in slope at $\sim$7000 K and $\sim$90 GPa, in agreement with ab initio simulations we have performed. This feature coincides with the gradual transition from a molecular liquid to a plasma state. Additionally, we performed reflectivity measurements, providing the first experimental evidence of electronic conduction in high-pressure ammonia. Shock reflectance continuously rises with pressure above 50 GPa and reaches saturation values above 120 GPa. Corresponding electrical conductivity values are up to 1 order of magnitude higher than in water in the 100 GPa regime, with possible significant contributions of the predicted ammonia-rich layers to the generation of magnetic dynamos in ice giant interiors., Comment: Please visit publisher's website for supplementary information

Efficient sampling of complex data distributions can be achieved using trained invertible flows (IF), where the model distribution is generated by pushing a simple base distribution through multiple non-linear bijective transformations. However, the iterative nature of the transformations in IFs can limit the approximation to the target distribution. In this paper we seek to mitigate this by implementing RBM-Flow, an IF model whose base distribution is a Restricted Boltzmann Machine (RBM) with a continuous smoothing applied. We show that by using RBM-Flow we are able to improve the quality of samples generated, quantified by the Inception Scores (IS) and Frechet Inception Distance (FID), over baseline models with the same IF transformations, but with less expressive base distributions. Furthermore, we also obtain D-Flow, an IF model with uncorrelated discrete latent variables. We show that D-Flow achieves similar likelihoods and FID/IS scores to those of a typical IF with Gaussian base variables, but with the additional benefit that global features are meaningfully encoded as discrete labels in the latent space.

We investigate traffic routing both from the perspective of theory as well as real world data. First, we introduce a new type of games: $\theta$-free flow games. Here, commuters only consider, in their strategy sets, paths whose free-flow costs (informally their lengths) are within a small multiplicative $(1+\theta)$ constant of the optimal free-flow cost path connecting their source and destination, where $\theta\geq0$. We provide an exhaustive analysis of tight bounds on PoA($\theta$) for arbitrary classes of cost functions, both in the case of general congestion/routing games as well as in the special case of path-disjoint networks. Second, by using a large mobility dataset in Singapore, we inspect minute-by-minute decision-making of thousands of commuters, and find that $\theta=1$ is a good estimate of agents' route (pre)selection mechanism. In contrast, in Pigou networks, the ratio of the free-flow costs of the routes, and thus $\theta$, is \textit{infinite}; so, although such worst case networks are mathematically simple, they correspond to artificial routing scenarios with little resemblance to real world conditions, opening the possibility of proving much stronger Price of Anarchy guarantees by explicitly studying their dependency on $\theta$. For example, in the case of the standard Bureau of Public Roads (BPR) cost model, where$c_e(x)= a_e x^4+b_e$, and for quartic cost functions in general, the standard PoA bound for $\theta=\infty$ is $2.1505$, and this is tight both for general networks as well as path-disjoint and even parallel-edge networks. In comparison, for $\theta=1$, the PoA in the case of general networks is only $1.6994$, whereas for path-disjoint/parallel-edge networks is even smaller ($1.3652$), showing that both the route geometries as captured by the parameter $\theta$ as well as the network topology have significant effects on PoA., Comment: Extended version of the conference paper presented in WINE 2020: The 16th Conference on Web and Internet Economics

In load balancing problems there is a set of clients, each wishing to select a resource from a set of permissible ones, in order to execute a certain task. Each resource has a latency function, which depends on its workload, and a client's cost is the completion time of her chosen resource. Two fundamental variants of load balancing problems are {\em selfish load balancing} (aka. {\em load balancing games}), where clients are non-cooperative selfish players aimed at minimizing their own cost solely, and {\em online load balancing}, where clients appear online and have to be irrevocably assigned to a resource without any knowledge about future requests. We revisit both selfish and online load balancing under the objective of minimizing the {\em Nash Social Welfare}, i.e., the geometric mean of the clients' costs. To the best of our knowledge, despite being a celebrated welfare estimator in many social contexts, the Nash Social Welfare has not been considered so far as a benchmarking quality measure in load balancing problems. We provide tight bounds on the price of anarchy of pure Nash equilibria and on the competitive ratio of the greedy algorithm under very general latency functions, including polynomial ones. For this particular class, we also prove that the greedy strategy is optimal as it matches the performance of any possible online algorithm.

In the adaptive influence maximization problem, we are given a social network and a budget $k$, and we iteratively select $k$ nodes, called seeds, in order to maximize the expected number of nodes that are reached by an influence cascade that they generate according to a stochastic model for influence diffusion. Differently from the non-adaptive influence maximization problem, where all the seeds must be selected beforehand, here nodes are selected sequentially one by one, and the decision on the $i$th seed is based on the observed cascade produced by the first $i-1$ seeds. We focus on the myopic feedback model, in which we can only observe which neighbors of previously selected seeds have been influenced and on the independent cascade model, where each edge is associated with an independent probability of diffusing influence. Previous works showed that the adaptivity gap is at most $4$, which implies that the non-adaptive greedy algorithm guarantees an approximation factor of $\frac{1}{4}\left(1-\frac{1}{e}\right)$ for the adaptive problem. In this paper, we improve the bounds on both the adaptivity gap and on the approximation factor. We directly analyze the approximation factor of the non-adaptive greedy algorithm, without passing through the adaptivity gap, and show that it is at least $\frac{1}{2}\left(1-\frac{1}{e}\right)$. Therefore, the adaptivity gap is at most $\frac{2e}{e-1}\approx 3.164$. To prove these bounds, we introduce a new approach to relate the greedy non-adaptive algorithm to the adaptive optimum. The new approach does not rely on multi-linear extensions or random walks on optimal decision trees, which are commonly used techniques in the field. We believe that it is of independent interest and may be used to analyze other adaptive optimization problems., Comment: arXiv admin note: text overlap with arXiv:2006.15374

In the influence maximization (IM) problem, we are given a social network and a budget $k$, and we look for a set of $k$ nodes in the network, called seeds, that maximize the expected number of nodes that are reached by an influence cascade generated by the seeds, according to some stochastic model for influence diffusion. In this paper, we study the adaptive IM, where the nodes are selected sequentially one by one, and the decision on the $i$th seed can be based on the observed cascade produced by the first $i-1$ seeds. We focus on the full-adoption feedback in which we can observe the entire cascade of each previously selected seed and on the independent cascade model where each edge is associated with an independent probability of diffusing influence. Our main result is the first sub-linear upper bound that holds for any graph. Specifically, we show that the adaptivity gap is upper-bounded by $\lceil n^{1/3}\rceil $, where $n$ is the number of nodes in the graph. Moreover, we improve over the known upper bound for in-arborescences from $\frac{2e}{e-1}\approx 3.16$ to $\frac{2e^2}{e^2-1}\approx 2.31$. Finally, we study $\alpha$-bounded graphs, a class of undirected graphs in which the sum of node degrees higher than two is at most $\alpha$, and show that the adaptivity gap is upper-bounded by $\sqrt{\alpha}+O(1)$. Moreover, we show that in 0-bounded graphs, i.e. undirected graphs in which each connected component is a path or a cycle, the adaptivity gap is at most $\frac{3e^3}{e^3-1}\approx 3.16$. To prove our bounds, we introduce new techniques to relate adaptive policies with non-adaptive ones that might be of their own interest., Comment: 18 pages

A finite length ferromagnetic chain with opposite spin polarisation imposed at its two ends is one of the simplest frustrated spin models. In the clean classical limit the domain wall inserted on account of the boundary conditions resides with equal probability on any one of the bonds, and the degeneracy is precisely equal to the number of bonds. If quantum mechanics is introduced via a transverse field, the domain wall will behave as a particle in a box, and prefer to be nearer the middle of the chain rather than the ends. A simple characteristic of a real quantum annealer is therefore which of these limits obtains in practice. Here we have used the ferromagnetic chain with antiparallel boundary spins to test a real flux qubit quantum annealer and discover that contrary to both expectations, the domain walls found are non-uniformly distributed on account of effective random longitudinal fields present notwithstanding tuning carried out to zero out such fields when the couplings between qubits are nominally zero. We present a simple derivation of the form of the distribution function for the domain walls, and show also how the effect we have discovered can be used to determine the strength of the effective random fields (noise) characterising the annealer. The noise measured in this fashion is smaller than what is seen during the single qubit tuning process, but nonetheless qualitatively affects the outcome of the simulation performed by the annealer., Comment: 14 pages, 9 figures plus supplemental material, title change in v2 and v3, accepted in npj Quantum Information, arXiv matches accepted version

Recent progress in quantum algorithms and hardware indicates the potential importance of quantum computing in the near future. However, finding suitable application areas remains an active area of research. Quantum machine learning is touted as a potential approach to demonstrate quantum advantage within both the gate-model and the adiabatic schemes. For instance, the Quantum-assisted Variational Autoencoder has been proposed as a quantum enhancement to the discrete VAE. We extend on previous work and study the real-world applicability of a QVAE by presenting a proof-of-concept for similarity search in large-scale high-dimensional datasets. While exact and fast similarity search algorithms are available for low dimensional datasets, scaling to high-dimensional data is non-trivial. We show how to construct a space-efficient search index based on the latent space representation of a QVAE. Our experiments show a correlation between the Hamming distance in the embedded space and the Euclidean distance in the original space on the Moderate Resolution Imaging Spectroradiometer (MODIS) dataset. Further, we find real-world speedups compared to linear search and demonstrate memory-efficient scaling to half a billion data points.

The shaping of astrophysical outflows into bright, dense and collimated jets due to magnetic pressure is here investigated using laboratory experiments. We notably look at the impact on jet collimation of a misalignment between the outflow, as it stems from the source, and the magnetic field. For small misalignments, a magnetic nozzle forms and redirects the outflow in a collimated jet. For growing misalignments, this nozzle becomes increasingly asymmetric, disrupting jet formation. Our results thus suggest outflow/magnetic field misalignment to be a plausible key process regulating jet collimation in a variety of objects from our Sun's outflows to extragalatic jets. Furthermore, they provide a possible interpretation for the observed structuring of astrophysical jets. Jet modulation could be interpreted as the signature of changes over time in the outflow/ambient field angle, and the change in the direction of the jet could be the signature of changes in the direction of the ambient field.

A new experimental platform based on laser-plasma interaction is proposed to explore the fundamental processes of wave coupling at the origin of interplanetary radio emissions. It is applied to the study of electromagnetic (EM) emission at twice the plasma frequency ($2\omega_p$) observed during solar bursts and thought to result from the coalescence of two Langmuir waves (LWs). In the interplanetary medium, the first LW is excited by electron beams, while the second is generated by electrostatic decay of Langmuir waves. In the present experiment, instead of an electron beam, an energetic laser propagating through a plasma excites the primary LW, with characteristics close to those at near-Earth orbit. The EM radiation at $2\omega_p$ is observed at different angles. Its intensity, spectral evolution and polarization confirm the LW-coalescence scenario., Comment: Accepted for publication in Phys. Rev. Lett

Graphical model selection is a seemingly impossible task when many pairs of variables are never jointly observed; this requires inference of conditional dependencies with no observations of corresponding marginal dependencies. This under-explored statistical problem arises in neuroimaging, for example, when different partially overlapping subsets of neurons are recorded in non-simultaneous sessions. We call this statistical challenge the "Graph Quilting" problem. We study this problem in the context of sparse inverse covariance learning, and focus on Gaussian graphical models where we show that missing parts of the covariance matrix yields an unidentifiable precision matrix specifying the graph. Nonetheless, we show that, under mild conditions, it is possible to correctly identify edges connecting the observed pairs of nodes. Additionally, we show that we can recover a minimal superset of edges connecting variables that are never jointly observed. Thus, one can infer conditional relationships even when marginal relationships are unobserved, a surprising result! To accomplish this, we propose an $\ell_1$-regularized partially observed likelihood-based graph estimator and provide performance guarantees in population and in high-dimensional finite-sample settings. We illustrate our approach using synthetic data, as well as for learning functional neural connectivity from calcium imaging data., Comment: 6 figures

The development of quantum-classical hybrid (QCH) algorithms is critical to achieve state-of-the-art computational models. A QCH variational autoencoder (QVAE) was introduced in Ref. [1] by some of the authors of this paper. QVAE consists of a classical auto-encoding structure realized by traditional deep neural networks to perform inference to, and generation from, a discrete latent space. The latent generative process is formalized as thermal sampling from either a quantum or classical Boltzmann machine (QBM or BM). This setup allows quantum-assisted training of deep generative models by physically simulating the generative process with quantum annealers. In this paper, we have successfully employed D-Wave quantum annealers as Boltzmann samplers to perform quantum-assisted, end-to-end training of QVAE. The hybrid structure of QVAE allows us to deploy current-generation quantum annealers in QCH generative models to achieve competitive performance on datasets such as MNIST. The results presented in this paper suggest that commercially available quantum annealers can be deployed, in conjunction with well-crafted classical deep neutral networks, to achieve competitive results in unsupervised and semisupervised tasks on large-scale datasets. We also provide evidence that our setup is able to exploit large latent-space (Q)BMs, which develop slowly mixing modes. This expressive latent space results in slow and inefficient classical sampling, and paves the way to achieve quantum advantage with quantum annealing in realistic sampling applications., Comment: 20 pages, 14 figures

Due to superstition, license plates with desirable combinations of characters are highly sought after in China, fetching prices that can reach into the millions in government-held auctions. Despite the high stakes involved, there has been essentially no attempt to provide price estimates for license plates. We present an end-to-end neural network model that simultaneously predict the auction price, gives the distribution of prices and produces latent feature vectors. While both types of neural network architectures we consider outperform simpler machine learning methods, convolutional networks outperform recurrent networks for comparable training time or model complexity. The resulting model powers our online price estimator and search engine.

Constructing powerful generative models for natural images is a challenging task. PixelCNN models capture details and local information in images very well but have limited receptive field. Variational autoencoders with a factorial decoder can capture global information easily, but they often fail to reconstruct details faithfully. PixelVAE combines the best features of the two models and constructs a generative model that is able to learn local and global structures. Here we introduce PixelVAE++, a VAE with three types of latent variables and a PixelCNN++ for the decoder. We introduce a novel architecture that reuses a part of the decoder as an encoder. We achieve the state of the art performance on binary data sets such as MNIST and Omniglot and achieve the state of the art performance on CIFAR-10 among latent variable models while keeping the latent variables informative.

For the last few years, the NASA Quantum Artificial Intelligence Laboratory (QuAIL) has been performing research to assess the potential impact of quantum computers on challenging computational problems relevant to future NASA missions. A key aspect of this research is devising methods to most effectively utilize emerging quantum computing hardware. Research questions include what experiments on early quantum hardware would give the most insight into the potential impact of quantum computing, the design of algorithms to explore on such hardware, and the development of tools to minimize the quantum resource requirements. We survey work relevant to these questions, with a particular emphasis on our recent work in quantum algorithms and applications, in elucidating mechanisms of quantum mechanics and their uses for quantum computational purposes, and in simulation, compilation, and physics-inspired classical algorithms. To our early application thrusts in planning and scheduling, fault diagnosis, and machine learning, we add thrusts related to robustness of communication networks and the simulation of many-body systems for material science and chemistry. We provide a brief update on quantum annealing work, but concentrate on gate-model quantum computing research advances within the last couple of years., Comment: 20 pages plus extensive references, 3 figures

Particle-In-Cell (PIC) codes are broadly applied to the kinetic simulation of plasmas, from laser-matter interaction to astrophysics. Their heavy simulation cost can be mitigated by using the Single Instruction Multiple Data (SIMD) capibility, or vectorization, now available on most architectures. This article details and discusses the vectorization strategy developed in the code Smilei which takes advantage from an efficient, systematic, cell-based sorting of the particles. The PIC operators on particles (projection, push, deposition) have been optimized to benefit from large SIMD vectors on both recent and older architectures. The efficiency of these vectorized operations increases with the number of particles per cell (PPC), typically speeding up three-dimensional simulations by a factor 2 with 256 PPC. Although this implementation shows acceleration from as few as 8 PPC, it can be slower than the scalar version in domains containing fewer PPC as usually observed in vectorization attempts. This issue is overcome with an adaptive algorithm which switches locally between scalar (for few PPC) and vectorized operators (otherwise). The newly implemented methods are benchmarked on three different, large-scale simulations considering configurations frequently studied with PIC codes.

Characterising materials at pressures of several megabar and temperatures of a few thousand Kelvin is critical for the understanding of the Warm Dense Matter regime and to improve planetary models as these conditions are typical of planets' interiors. Today, laser-driven shock compression is the only technique to achieve multimegabar pressures, but the associated temperatures are too high to be representative of planetary states. Double-shock compression represents an alternative to explore lower temperatures. Here we present a method to create well-controlled double-shocked states and measure their thermodynamic state and optical reflectivity using standard optical diagnostics (Doppler velocimetry and optical pyrometry) in a laser-driven shock experiment. This method, which does not require the support of hydrodynamical simulations, is based on the application of generalised Rankine-Hugoniot relations together with a self impedance mismatch technique. A validation experiment has been performed at the LULI2000 facility (\'Ecole Polytechnique, France) on a sample of $\alpha$-quartz. A temperature $60 \%$ lower than along the principal Hugoniot has been obtained at 7.5 Mbar., Comment: 10 pages, 7 figures

We study the existence of allocations of indivisible goods that are envy-free up to one good (EF1), under the additional constraint that each bundle needs to be connected in an underlying item graph. If the graph is a path and the utility functions are monotonic over bundles, we show the existence of EF1 allocations for at most four agents, and the existence of EF2 allocations for any number of agents; our proofs involve discrete analogues of the Stromquist's moving-knife protocol and the Su--Simmons argument based on Sperner's lemma. For identical utilities, we provide a polynomial-time algorithm that computes an EF1 allocation for any number of agents. For the case of two agents, we characterize the class of graphs that guarantee the existence of EF1 allocations as those whose biconnected components are arranged in a path; this property can be checked in linear time., Comment: Accepted journal version

Water, ethanol, and ammonia are the key components of the mantles of Uranus and Neptune. To improve structure and evolution models and give an explanation of the magnetic fields and luminosities of the icy giants, those components need to be characterised at planetary conditions (some Mbar and a few $10^3$ K). Those conditions are typical of the Warm Dense Matter regime, which exhibits a rich phase diagram, with the coexistence of many states of matter and a large variety of chemical processes. H$_2$O, C:H:O, and C:H:N:O mixtures have been compressed up to 2.8 Mbar along the principal Hugoniot using laser-driven decaying shocks. The experiments were performed at the GEKKO XII and LULI 2000 laser facilities using standard optical diagnostics (Doppler velocimetry and pyrometry) to characterise equation of state and optical reflectivity of the shocked states. The results show that H$_2$O and the C:H:N:O mixture share the same equation of state with a density scaling, while the reflectivity behaves differently by what concerns both the onset pressures and the saturation values. The reflectivity measurement at two frequencies allows to estimate the conductivity and the complex refractive index using a Drude model.

Quantum annealing in a real device is necessarily susceptible to errors due to diabatic transitions and thermal noise. Nested quantum annealing correction is a method to suppress errors by using an all-to-all penalty coupling among a set of physical qubits representing a logical qubit. We show analytically that nested quantum annealing correction can suppress errors effectively in ferromagnetic and antiferromagnetic Ising models with infinite-range interactions. Our analysis reveals that the nesting structure can significantly weaken or even remove first-order phase transitions, in which the energy gap closes exponentially. The nesting structure also suppresses thermal fluctuations by reducing the effective temperature., Comment: 28 pages