Search

We consider the problem of private computation (PC) in a distributed storage system. In such a setting a user wishes to compute a function of $f$ messages replicated across $n$ noncolluding databases, while revealing no information about the desired function to the databases. We provide an information-theoretically accurate achievable PC rate, which is the ratio of the smallest desired amount of information and the total amount of downloaded information, for the scenario of nonlinear computation. For a large message size the rate equals the PC capacity, i.e., the maximum achievable PC rate, when the candidate functions are the $f$ independent messages and one arbitrary nonlinear function of these. When the number of messages grows, the PC rate approaches an outer bound on the PC capacity. As a special case, we consider private monomial computation (PMC) and numerically compare the achievable PMC rate to the outer bound for a finite number of messages., Comment: 5 pages, 1 figure, 1 table. Presented at the 2019 IEEE Information Theory Workshop (ITW). Figure 1 is updated as it contained incorrect data-points for $f=2$ and $g=3$. arXiv admin note: text overlap with arXiv:2003.10007

We present a general lower bound for the fundamental tone for the p-Laplacian on Riemannian manifolds carrying a special kind of function. We then apply our result to the cases of negatively curved simply connected manifolds, a class of warped product manifolds and for a class of Riemannian submersions.

The object of study is the problem of testing for uniformity of the multinomial distribution. We consider tests based on symmetric statistics, defined as the sum of some function of cell-frequencies. Mainly, attention is focused on the class of power divergence statistics, in particular, on the chi-square and log-likelihood ratio statistics. The main issue of the article is to study the asymptotic properties of tests at the concept of an intermediate setting in terms of so called -intermediate asymptotic efficiency due to Ivchenko and Mirakhmedov (1995), when the asymptotic power of tests are bounded away from zero and one, while sequences of alternatives converge to the hypothesis, but not too fast.

We put forward a model that describes a decaying and evanescent point source of non-interacting quantum waves in 1D. This point-source assumption allows for a simple description that captures the essential aspects of the dynamics of a wave traveling through a classically forbidden region without the need to specify the details of the inner region. The dynamics of the resulting wave is examined and several characteristic times are identified. One of them generalizes the tunneling time-scale introduced by Büttiker and Landauer and it characterizes the arrival of the maximum of the wave function. Diffraction in time and deviations from exponential decay are also studied. Here we show that there exists an optimal injection frequency and detection point for the observation of these two quantum phenomena., This work was supported by the Basque Country Government (Grant no. IT472-10), Ministerio de Economía y Competitividad (Grant no. FIS2012-36673-C03-01), and the program UFI 11/55 of UPV/EHU

The deep convolutional neural network has achieved significant progress for single image rain streak removal. However, most of the data-driven learning methods are full-supervised or semi-supervised, unexpectedly suffering from significant performance drops when dealing with real rain. These data-driven learning methods are representative yet generalize poor for real rain. The opposite holds true for the model-driven unsupervised optimization methods. To overcome these problems, we propose a unified unsupervised learning framework which inherits the generalization and representation merits for real rain removal. Specifically, we first discover a simple yet important domain knowledge that directional rain streak is anisotropic while the natural clean image is isotropic, and formulate the structural discrepancy into the energy function of the optimization model. Consequently, we design an optimization model-driven deep CNN in which the unsupervised loss function of the optimization model is enforced on the proposed network for better generalization. In addition, the architecture of the network mimics the main role of the optimization models with better feature representation. On one hand, we take advantage of the deep network to improve the representation. On the other hand, we utilize the unsupervised loss of the optimization model for better generalization. Overall, the unsupervised learning framework achieves good generalization and representation: unsupervised training (loss) with only a few real rainy images (input) and physical meaning network (architecture). Extensive experiments on synthetic and real-world rain datasets show the superiority of the proposed method., Comment: Accept to 2021ACMMM

We consider an automatic construction of locally optimal preconditioners for positive definite linear systems. To achieve this goal, we introduce a differentiable loss function that does not explicitly include the estimation of minimal eigenvalue. Nevertheless, the resulting optimization problem is equivalent to a direct minimization of the condition number. To demonstrate our approach, we construct a parametric family of modified BPX preconditioners. Namely, we define a set of empirical basis functions for coarse finite element spaces and tune them to achieve better condition number. For considered model equations (that includes Poisson, Helmholtz, Convection–diffusion, Biharmonic, and others), we achieve from two to twenty times smaller condition numbers for symmetric positive definite linear systems.

We consider the platoon matching problem for a set of trucks with the same origin, but different destinations. It is assumed that the vehicles benefit from traveling in a platoon for instance through reduced fuel consumption. The vehicles belong to different fleet owners and their strategic interaction is modeled as a non-cooperative game where the vehicle actions are their departure times. Each truck has a preferred departure time and its utility function is defined as the difference between its benefit from platooning and the cost of deviating from its preferred departure time. We show that the platoon matching game is an exact potential game. An algorithm based on best response dynamics is proposed for finding a Nash equilibrium of the game. At a Nash equilibrium, vehicles with the same departure time are matched to form a platoon. Finally, the total fuel reduction at the Nash equilibrium is studied and compared with that of a cooperative matching solution where a common utility function for all vehicles is optimized., Comment: ITSC 2018

This paper is concerned with the variational problem for the elastic energy defined on symmetric graphs under the unilateral constraint. Assuming that the obstacle function satisfies the symmetric cone condition, we prove (i) uniqueness of minimizers, (ii) loss of regularity of minimizers, and give (iii) complete classification of existence and non-existence of minimizers in terms of the size of obstacle. As an application, we characterize the solution of the obstacle problem as equilibrium of the corresponding dynamical problem., Comment: 30 pages, 5 figures

Unsupervised domain adaptation (UDA) requires source domain samples with clean ground truth labels during training. Accurately labeling a large number of source domain samples is time-consuming and laborious. An alternative is to utilize samples with noisy labels for training. However, training with noisy labels can greatly reduce the performance of UDA. In this paper, we address the problem that learning UDA models only with access to noisy labels and propose a novel method called robust local preserving and global aligning network (RLPGA). RLPGA improves the robustness of the label noise from two aspects. One is learning a classifier by a robust informative-theoretic-based loss function. The other is constructing two adjacency weight matrices and two negative weight matrices by the proposed local preserving module to preserve the local topology structures of input data. We conduct theoretical analysis on the robustness of the proposed RLPGA and prove that the robust informative-theoretic-based loss and the local preserving module are beneficial to reduce the empirical risk of the target domain. A series of empirical studies show the effectiveness of our proposed RLPGA., Comment: Accepted by IEEE Transactions on Knowledge and Data Engineering (TKDE) 2022; Refer to https://ieeexplore.ieee.org/document/9540279

Recently, the hyperspectral sensors have improved our ability to monitor the earth surface with high spectral resolution. However, the high dimensionality of spectral data brings challenges for the image processing. Consequently, the dimensionality reduction is a necessary step in order to reduce the computational complexity and increase the classification accuracy. In this paper, we propose a new filter approach based on information gain for dimensionality reduction and classification of hyperspectral images. A special strategy based on hyperspectral bands selection is adopted to pick the most informative bands and discard the irrelevant and noisy ones. The algorithm evaluates the relevancy of the bands based on the information gain function with the support vector machine classifier. The proposed method is compared using two benchmark hyperspectral datasets (Indiana, Pavia) with three competing methods. The comparison results showed that the information gain filter approach outperforms the other methods on the tested datasets and could significantly reduce the computation cost while improving the classification accuracy. Keywords: Hyperspectral images; dimensionality reduction; information gain; classification accuracy. Keywords: Hyperspectral images; dimensionality reduction; information gain; classification accuracy.

Usually, opinion formation models assume that individuals have an opinion about a given topic which can change due to interactions with others. However, individuals can have different opinions in different topics and therefore n-dimensional models are best suited to deal with these cases. While there have been many efforts to develop analytical models for one dimensional opinion models, less attention has been paid to multidimensional ones. In this work, we develop an analytical approach for multidimensional models of continuous opinions where dimensions can be correlated or uncorrelated. We show that for any generic reciprocal interactions between agents, the mean value of initial opinion distribution is conserved. Moreover, for positive social influence interaction mechanisms, the variance of opinion distributions decreases with time and the system converges to a delta distributed function. In particular, we calculate the convergence time when agents get closer in a discrete quantity after interacting, showing a clear difference between correlated and uncorrelated cases.

This paper investigates repeated win-lose coordination games (WLC-games). We analyse which protocols are optimal for these games, covering both the worst case and average case scenarios, i,e., optimizing the guaranteed and expected coordination times. We begin by analysing Choice Matching Games (CM-games) which are a simple yet fundamental type of WLC-games, where the goal of the players is to pick the same choice from a finite set of initially indistinguishable choices. We give a fully complete classification of optimal expected and guaranteed coordination times in two-player CM-games and show that the corresponding optimal protocols are unique in every case - except in the CM-game with four choices, which we analyse separately. Our results on CM-games are essential for proving a more general result on the difficulty of all WLC-games: we provide a complete analysis of least upper bounds for optimal expected coordination times in all two-player WLC-games as a function of game size. We also show that CM-games can be seen as the most difficult games among all two-player WLC-games, as they turn out to have the greatest optimal expected coordination times., Comment: In Proceedings GandALF 2021, arXiv:2109.07798. A preprint with full proofs and further examples can be found at arXiv:2004.07381v1

Transition-metal dichalcogenides (TMDs) are promising for two-dimensional (2D) semiconducting devices and novel phenomena. For 2D applications, their work function, ionization energy, and electron affinity are required as a function of thickness, but research on this is yet to cover the full family of compounds. Here, we present the work function, ionization energy, and electron affinity of few-layer and bulk MX2 (M = Mo, W and X = S, Se, Te) in 2H phase obtained accurately by the density functional theory and GW calculations. For each compound, we consider one-, two-, three-, four-layer, and bulk geometry. In GW calculations, accurate results are obtained by nonuniform q sampling for two-dimensional geometry. From band energies including the GW self-energy correction, we estimate the work function, band gap, ionization energy, and electron affinity as functions of the number of layers. We compare our results with available theoretical and experimental reports, and we discuss types of band alignments in in-plane and out-of-plane junctions of these few-layer and bulk TMDs., Comment: 10 pages, 5 figures, 7 tables

The wavelet analysis technique is a powerful tool and is widely used in broad disciplines of engineering, technology, and sciences. In this work, we present a novel scheme of constructing continuous wavelet functions, in which the wavelet functions are obtained by taking the first derivative of smoothing functions with respect to the scale parameter. Due to this wavelet constructing scheme, the inverse transforms are only one-dimensional integrations with respect to the scale parameter, and hence the continuous wavelet transforms constructed in this way are more ready to use than the usual scheme. We then apply the Gaussian-derived wavelet constructed by our scheme to computations of the density power spectrum for dark matter, the velocity power spectrum and the kinetic energy spectrum for baryonic fluid. These computations exhibit the convenience and strength of the continuous wavelet transforms. The transforms are very easy to perform, and we believe that the simplicity of our wavelet scheme will make continuous wavelet transforms very useful in practice., Comment: 9 Pages, 4 Figures, comments welcome

Lebesgue integration of derivatives of strongly-oscillatory functions is a recurring challenge in computational science and engineering. Integration by parts is an effective remedy for huge computational costs associated with Monte Carlo integration schemes. In case of Lebesgue integrals over a smooth manifold, however, integration by parts gives rise to a derivative of the density implied by charts describing the domain manifold. This paper focuses on the computation of that derivative, which we call the density gradient function, on general smooth manifolds. We analytically derive formulas for the density gradient and present examples of manifolds determined by popular differential equation-driven systems. We highlight the significance of the density gradient by demonstrating a numerical example of Monte Carlo integration involving oscillatory integrands., Comment: 24 pages, 11 figures, submitted to journal

We study from the proof complexity perspective the (informal) proof search problem: Is there an optimal way to search for propositional proofs? We note that for any fixed proof system there exists a time-optimal proof search algorithm. Using classical proof complexity results about reflection principles we prove that a time-optimal proof search algorithm exists without restricting proof systems iff a p-optimal proof system exists. To characterize precisely the time proof search algorithms need for individual formulas we introduce a new proof complexity measure based on algorithmic information concepts. In particular, to a proof system $P$ we attach {\bf information-efficiency function} $i_P(\tau)$ assigning to a tautology a natural number, and we show that: - $i_P(\tau)$ characterizes time any $P$-proof search algorithm has to use on $\tau$ and that for a fixed $P$ there is such an information-optimal algorithm, - a proof system is information-efficiency optimal iff it is p-optimal, - for non-automatizable systems $P$ there are formulas $\tau$ with short proofs but having large information measure $i_P(\tau)$. We isolate and motivate the problem to establish unconditional super-logarithmic lower bounds for $i_P(\tau)$ where no super-polynomial size lower bounds are known. We also point out connections of the new measure with some topics in proof complexity other than proof search., Comment: Preliminary version February 2021

Overlap functions are a class of aggregation functions that measure the overlapping degree between two values. Interval-valued overlap functions were defined as an extension to express the overlapping of interval-valued data, and they have been usually applied when there is uncertainty regarding the assignment of membership degrees. The choice of a total order for intervals can be significant, which motivated the recent developments on interval-valued aggregation functions and interval-valued overlap functions that are increasing to a given admissible order, that is, a total order that refines the usual partial order for intervals. Also, width preservation has been considered on these recent works, in an intent to avoid the uncertainty increase and guarantee the information quality, but no deeper study was made regarding the relation between the widths of the input intervals and the output interval, when applying interval-valued functions, or how one can control such uncertainty propagation based on this relation. Thus, in this paper we: (i) introduce and develop the concepts of width-limited interval-valued functions and width limiting functions, presenting a theoretical approach to analyze the relation between the widths of the input and output intervals of bivariate interval-valued functions, with special attention to interval-valued aggregation functions; (ii) introduce the concept of $(a,b)$-ultramodular aggregation functions, a less restrictive extension of one-dimension convexity for bivariate aggregation functions, which have an important predictable behaviour with respect to the width when extended to the interval-valued context; (iii) define width-limited interval-valued overlap functions, taking into account a function that controls the width of the output interval; (iv) present and compare three construction methods for these width-limited interval-valued overlap functions., Comment: submitted

We introduce a machine learning model, the q-CNN model, sharing key features with convolutional neural networks and admitting a tensor network description. As examples, we apply q-CNN to the MNIST and Fashion MNIST classification tasks. We explain how the network associates a quantum state to each classification label, and study the entanglement structure of these network states. In both our experiments on the MNIST and Fashion-MNIST datasets, we observe a distinct increase in both the left/right as well as the up/down bipartition entanglement entropy (EE) during training as the network learns the fine features of the data. More generally, we observe a universal negative correlation between the value of the EE and the value of the cost function, suggesting that the network needs to learn the entanglement structure in order the perform the task accurately. This supports the possibility of exploiting the entanglement structure as a guide to design the machine learning algorithm suitable for given tasks.

Resolvent analysis has demonstrated encouraging results for modeling coherent structures in jets when compared against their data-educed counterparts from high-fidelity large-eddy simulations (LES). We formulate resolvent analysis as an acoustic analogy that relates the near-field resolvent forcing to the near- and far-field pressure. We use an LES database of round, isothermal, Mach 0.9 and 1.5 jets to produce an ensemble of realizations for the acoustic field that we project onto a limited set of resolvent modes. In the near-field, we perform projections on a restricted acoustic output domain, $r/D = [5,6]$, while the far-field projections are performed on a Kirchhoff surface comprising a 100-diameter arc centered at the nozzle. This allows the LES realizations to be expressed in the resolvent basis via a data-deduced, low-rank, cross-spectral density matrix. We find that a single resolvent mode reconstructs the most energetic regions of the acoustic field across Strouhal numbers, $St = [0-1]$, and azimuthal wavenumbers, $m=[0,2]$. Finally, we present a simple function that results in a rank-1 resolvent model agreeing within 2dB of the peak noise for both jets., Comment: 14 pages, 11 figures, Accepted to the Journal of the Acoustical Society of America

We study the motion of an impurity in a two-leg ladder interacting with two fermionic baths along each leg, a simple model bridging cold atom quantum simulators with an idealised description of the basic transport processes in a layered heterostructure. Using the linked-cluster expansion we obtain exact analytical results for the single-particle Green's function and find that the long-time behaviour is dominated by an intrinsic orthogonality catastrophe associated to the motion of the impurity in each one-dimensional chain. We explore both the case of two identical legs as well as the case where the legs are characterised by different interaction strengths: in the latter case we observe a subleading correction which can be relevant for intermediate-time transport at an interface between different materials. In all the cases we do not find significant differences between the intra- and inter-leg Green's functions in the long-time limit., Comment: 16 pages, 7 figures

The viscous dissipation between rigid, randomly rough indenters and linearly elastic counter bodies sliding past them is investigated using Green's function molecular dynamics. The study encompasses a variety of models differing in the height spectra properties of the rigid indenter, in the viscoelasticity of the elastomer, and in their interaction. All systems reveal the expected damping linear in sliding velocity $v$ at small $v$ and a pronounced maximum at intermediate $v$. Persson's theory of rubber friction, which is adopted to the studied model systems, reflects all observed trends. However, close quantitative agreement is only found up to intermediate sliding velocities. Relative errors in the friction force become significant once the contact area is substantially reduced by sliding., Comment: 13 pages, 9 figures

The paper analyzes the rotation averaging problem as a minimization problem for a potential function of the corresponding gradient system. This dynamical system is one generalization of the famous Kuramoto model on special orthogonal group SO(3), which is known as the non-Abelian Kuramoto model. We have proposed a novel method for finding weighted and unweighted rotation average. In order to verify the correctness of our algorithms, we have compared the simulation results with geometric and projected average using real and random data sets. In particular, we have discovered that our method gives approximately the same results as geometric average.

This paper introduces a collision avoidance system for navigating a multicopter in cluttered outdoor environments based on the recent memory-less motion planner, rectangular pyramid partitioning using integrated depth sensors (RAP-PIDS). The RAPPIDS motion planner generates collision-free flight trajectories at high speed with low computational cost using only the latest depth image. In this work we extend it to improve the performance of the planner by taking the following issues into account. (a) Changes in the dynamic characteristics of the multicopter that occur during flight, such as changes in motor input/output characteristics due to battery voltage drop. (b) The noise of the flight sensor, which can cause unwanted control input components. (c) Planner utility function which may not be suitable for the cluttered environment. Therefore, in this paper we introduce solutions to each of the above problems and propose a system for the successful operation of the RAPPIDS planner in an outdoor cluttered flight environment. At the end of the paper, we validate the proposed method's effectiveness by presenting the flight experiment results in a forest environment. A video can be found at www.youtube.com/watch?v=3av5xEuKg2w&ab_channel=HiPeRLab

Considering an energy density of the form [Formula: see text] (where [Formula: see text] is an arbitrary positive constant with dimension of energy density and [Formula: see text]), a shape function is obtained by using field equations of braneworld gravity theory in this paper. Under isotropic scenario wormhole solutions are obtained considering six different redshift functions along with the obtained new shape function. For anisotropic case, wormhole solutions are obtained under the consideration of five different shape functions along with the redshift function [Formula: see text], where [Formula: see text] is an arbitrary constant. In each case all energy conditions are examined and it is found that for some cases all energy conditions are satisfied in the vicinity of the wormhole throat and for the rest of the cases all energy conditions are satisfied except strong energy condition.

We present an extensive comparison between the statistical properties of non-hierarchical three-body systems and the corresponding three-body theoretical predictions. We perform and analyze 1 million realizations for each different initial condition considering equal and unequal mass three-body systems to provide high accuracy statistics. We measure 4 quantities characterizing the statistical distribution of ergodic disintegrations: escape probability of each body, the characteristic exponent for escapes by a narrow margin, predicted absorptivity as a function of binary energy and binary angular momentum, and, finally, the lifetime distribution. The escape probabilities are shown to be in agreement down to the 1% level with the emissivity-blind, flux-based theoretical prediction. This represents a leap in accuracy compared to previous three-body statistical theories. The characteristic exponent at the threshold for marginally unbound escapes is an emissivity-independent flux-based prediction, and the measured values are found to agree well with the prediction. We interpret both tests as strong evidence for the flux-based three-body statistical formalism. The predicted absorptivity and lifetime distributions are measured to enable future tests of statistical theories., Comment: 17 pages, 11 figures

We combine Newton's variational method with ideas from eigenvector continuation to construct a fast & accurate emulator for two-body scattering observables. The emulator will facilitate the application of rigorous statistical methods for interactions that depend smoothly on a set of free parameters. Our approach begins with a trial $K$ or $T$ matrix constructed from a small number of exact solutions to the Lippmann--Schwinger equation. Subsequent emulation only requires operations on small matrices. We provide several applications to short-range potentials with and without the Coulomb interaction and partial-wave coupling. It is shown that the emulator can accurately extrapolate far from the support of the training data. When used to emulate the neutron-proton cross section with a modern chiral interaction as a function of 26 free parameters, it reproduces the exact calculation with negligible error and provides an over 300x improvement in CPU time., 9 pages, 6 figures; minor updates, close to published version

We investigate a possibility to describe the non-Debye relaxation processes using the Volterra-type equations with kernels given by the Prabhakar functions with the upper parameter ν being negative. Proposed integro-differential equations mimic the fading memory effects and are explicitly solved using the umbral calculus and the Laplace transform methods. Both approaches lead to the same results valid for admissible domain of the parameters α, μ and ν characterizing the Prabhakar function. For the special case α ∈ (0, 1], μ = 0 and ν = − 1 we recover the Cole-Cole model, in general having a residual polarization. We also show that our scheme gives results equivalent to those obtained using the stochastic approach to relaxation phenomena merged with integral equations involving kernels given by the Prabhakar functions with the positive upper parameter.

We propose a novel approach to obtain the growth rate of cosmic structures, $f(z)$, from the evolution of the cosmic homogeneity scale, $R_{\text{H}}(z)$. Our methodology needs two ingredients in a specific functional form: $R_{\text{H}}(z)$ data and the matter two-point correlation function today, i.e., $\xi(r, z=0)$. We use a Gaussian Process approach to reconstruct the function $R_{\text{H}}$. In the absence of suitable observational information of the matter correlation function in the local Universe, $z \simeq 0$, we assume a fiducial cosmology to obtain $\xi(r, z=0)$. For this reason, our final result turns out to be a consistency test of the cosmological model assumed. Our results show a good agreement between: (i) the growth rate $f^{R_{\text{H}}}(z)$ obtained through our approach, (ii) the $f^{\Lambda\text{CDM}}(z)$ expected in the fiducial model, and (iii) the best-fit $f(z)$ from data compiled in the literature. Moreover, using this data compilation, we perform a Gaussian Process to reconstruct the growth rate function $f^{\text{data}}(z)$ and compare it with the function $f^{R_{\text{H}}}(z)$ finding a concordance of $< \!2 \,\sigma$, a good result considering the few data available for both reconstruction processes. With more accurate $R_{\text{H}}(z)$ data, from forthcoming surveys, the homogeneity scale function might be better determined and would have the potential to discriminate between $\Lambda$CDM and alternative scenarios as a new cosmological observable., Comment: Accepted for publication in MNRAS

[Formula: see text] gravity is a novel extension of Einsteinian cubic gravity in which the Ricci scalar in the action is replaced by a function of the curvature invariant [Formula: see text] which represents the contractions of the Riemann tensor at the cubic order [C. Erices, E. Papantonopoulos and E. N. Saridakis, Phys. Rev. D 99 (2019) 123527]. This work is concentrated on bounding some [Formula: see text] gravity models using the concept of energy conditions where the functional forms of [Formula: see text] are represented as (a) [Formula: see text] and (b) [Formula: see text], where [Formula: see text] is the sole model parameter. Energy conditions are interesting linear relationships between pressure and density and have been extensively employed to derive interesting results in Einstein’s gravity, and are also an excellent tool to impose constraints on any cosmological model. To place the bounds, we ensured that the energy density must remain positive, the pressure must remain negative, and the equation of state parameter must attain a value close to [Formula: see text] to make sure that the bounds respect the accelerated expansion of the universe and are also in harmony with the latest observational data. We report that for both the models, suitable parameter spaces exist which satisfy the aforementioned conditions, and therefore, posit the [Formula: see text] theory of gravity to be a promising modified theory of gravitation.

Generative Adversarial Networks are a popular method for learning distributions from data by modeling the target distribution as a function of a known distribution. The function, often referred to as the generator, is optimized to minimize a chosen distance measure between the generated and target distributions. One commonly used measure for this purpose is the Wasserstein distance. However, Wasserstein distance is hard to compute and optimize, and in practice entropic regularization techniques are used to improve numerical convergence. The influence of regularization on the learned solution, however, remains not well-understood. In this paper, we study how several popular entropic regularizations of Wasserstein distance impact the solution in a simple benchmark setting where the generator is linear and the target distribution is high-dimensional Gaussian. We show that entropy regularization promotes the solution sparsification, while replacing the Wasserstein distance with the Sinkhorn divergence recovers the unregularized solution. Both regularization techniques remove the curse of dimensionality suffered by Wasserstein distance. We show that the optimal generator can be learned to accuracy $\epsilon$ with $O(1/\epsilon^2)$ samples from the target distribution. We thus conclude that these regularization techniques can improve the quality of the generator learned from empirical data for a large class of distributions., Comment: 29 pages, 7 figures

This paper presents a generative model method for multispectral image fusion in remote sensing which involves training without supervision. This method eases the supervision of learning and also uses a multi-objective loss function to achieve image fusion. The loss function used incorporates both spectral and spatial distortions. Two discriminators are designed to minimize the spectral and spatial distortions of the generative output. Extensive experimentations are conducted using three public domain datasets. The comparison results across four reduced-resolution and three full-resolution objective metrics show the superiority of the developed method over several recently developed methods.

In this work, an adaptive edge element method is developed for an H(curl)-elliptic constrained optimal control problem. We use the lowest-order Nedelec's edge elements of first family and the piecewise (element-wise) constant functions to approximate the state and the control, respectively, and propose a new adaptive algorithm with error estimators involving both residual-type error estimators and lower-order data oscillations. By using a local regular decomposition for H(curl)-functions and the standard bubble function techniques, we derive the a posteriori error estimates for the proposed error estimators. Then we exploit the convergence properties of the orthogonal $L^2$-projections and the mesh-size functions to demonstrate that the sequences of the discrete states and controls generated by the adaptive algorithm converge strongly to the exact solutions of the state and control in the energy-norm and $L^2$-norm, respectively, by first achieving the strong convergence towards the solution to a limiting control problem. Three-dimensional numerical experiments are also presented to confirm our theoretical results and the quasi-optimality of the adaptive edge element method., Comment: submitted

A famous conjecture of Parkin-Shanks predicts that $p(n)$ is odd with density $1/2$. Despite the remarkable amount of work of the last several decades, however, even showing this density is positive seems out of reach. In a 2018 paper with Judge, we introduced a different approach and conjectured the "striking" fact that, if for any $A \equiv \pm 1\ (\bmod 6)$ the multipartition function $p_A(n)$ has positive odd density, then so does $p(n)$. Similarly, the positive odd density of any $p_{A}(n)$ with $A\equiv 3\ (\bmod 6)$ would imply that of $p_3(n)$. Our conjecture was shown to be a corollary of an earlier conjecture of the same paper. In this brief note, we provide an unconditional proof of it. An important tool will be Chen's recent breakthrough on a special case of our earlier conjecture., Comment: Minor updates. To appear in the J. of Number Theory

We derive generalized Fluctuation-Dissipation Relations (FDR) holding for a general stochastic dynamics that includes as subcases both equilibrium models for passive colloids and non-equilibrium models used to describe active particles. The relations reported here differ from previous formulations of the FDR because of their simplicity: they require only the microscopic knowledge of the dynamics instead of the whole expression of the steady-state probability distribution function that, except for linear interactions, is unknown for systems displaying non-vanishing currents. From the response function, we can extrapolate generalized versions of the Mesoscopic Virial equation and the equipartition theorem, which still holds far from equilibrium. Our results are tested in the case of equilibrium colloids described by underdamped or overdamped Langevin equations and for models describing the non-equilibrium behavior of active particles. Both the Active Brownian Particle and the Active Ornstein-Uhlenbeck particle models are compared in the case of a single particle confined in an external potential.

Since most industrial control applications use PID controllers, PID tuning and anti-windup measures are significant problems. This paper investigates tuning the feedback gains of a PID controller via back-calculation and automatic differentiation tools. In particular, we episodically use a cost function to generate gradients and perform gradient descent to improve controller performance. We provide a theoretical framework for analyzing this non-convex optimization and establish a relationship between back-calculation and disturbance feedback policies. We include numerical experiments on linear systems with actuator saturation to show the efficacy of this approach., Comment: Extension of paper in 2021 55th Annual Conference on Information Sciences and Systems (CISS). IEEE, 2021

We study the surface elastic response of pure Ni, the random alloy FeNiCr and an average FeNiCr alloy in terms of the surface lattice Green's function. We propose a scheme for computing per-site Green's function and study their per-site variations. The average FeNiCr alloy accurately reproduces the mean Green's function of the full random alloy. Variation around this mean is largest near the edge of the surface Brillouin-zone and decays as $q^{-2}$ with wavevector $q$ towards the $\Gamma$-point. We also present expressions for the continuum surface Green's function of anisotropic solids of finite and infinite thickness and show that the atomistic Green's function approaches continuum near the $\Gamma$-point. Our results are a first step towards efficient contact calculations and Peierls-Nabarro type models for dislocations in high-entropy alloys., Comment: 15 pages, 7 figures

We compare three first-principles methods of calculating the Curie temperature in two-dimensional (2D) ferromagnetic materials (FM), modeled using the Heisenberg model, and propose a simple formula for estimating the Curie temperature with high accuracy that works for all common 2D lattice types. First, we study the effect of exchange anisotropy on the Curie temperature calculated using the Monte Carlo (MC), the Green's function, and the renormalized spin-wave (RNSW) methods. We find that the Green's function method overestimates the Curie temperature in high-anisotropy regimes compared to the MC method, whereas the RNSW method underestimates the Curie temperature compared to the MC and the Green's function methods. Next, we propose a closed-form formula for calculating the Curie temperature of 2D FMs, which provides an estimate of the Curie temperature that is greatly improved over the mean-field expression for magnetic material screening. We apply the closed-form formula to predict the Curie temperature 2D magnets screened from the C2DB database and discover several high Curie temperature FMs, with ${\mathrm{Fe}}_{2}{\mathrm{F}}_{2}$ and ${\mathrm{MoI}}_{2}$ emerging as the most promising 2D ferromagnets. Finally, by comparing to experimental results for ${\mathrm{CrI}}_{3}$, ${\mathrm{CrCl}}_{3}$, and ${\mathrm{CrBr}}_{3}$, we conclude that for small effective anisotropies, the Green's-function-based equations are preferable, while for larger anisotropies, MC-based results are more predictive.

The cyber security of smart grids has become one of key problems in developing reliable modern power and energy systems. This paper introduces a non-stationary adversarial cost with a variation constraint for smart grids and enables us to investigate the problem of optimal smart grid protection against cyber attacks in a relatively practical scenario. In particular, a Bayesian multi-node bandit (MNB) model with adversarial costs is constructed and a new regret function is defined for this model. An algorithm called Thompson–Hedge algorithm is presented to solve the problem and the superior performance of the proposed algorithm is proven in terms of the convergence rate of the regret function. The applicability of the algorithm to real smart grid scenarios is verified and the performance of the algorithm is also demonstrated by numerical examples.

In the bulk reconstruction program, one constructs boundary representations of bulk fields. We investigate the relation between the global/Poincare and AdS-Rindler representations for AdS$_2$. We obtain the AdS-Rindler smearing function for massive and massless fields and show that the global and AdS-Rindler boundary representations are related by conformal transformations. We also use the boundary representations of creation and annihilation operators to compute the Bogoliubov transformation relating global modes to AdS-Rindler modes for both massive and massless particles., Comment: 20 pages

Shape reconstruction of deformable organs from two-dimensional X-ray images is a key technology for image-guided intervention. In this paper, we propose an image-to-graph convolutional network (IGCN) for deformable shape reconstruction from a single-viewpoint projection image. The IGCN learns relationship between shape/deformation variability and the deep image features based on a deformation mapping scheme. In experiments targeted to the respiratory motion of abdominal organs, we confirmed the proposed framework with a regularized loss function can reconstruct liver shapes from a single digitally reconstructed radiograph with a mean distance error of 3.6 mm., 24th International Conference, Strasbourg, France, September 27–October 1, 2021, Part of the Lecture Notes in Computer Science book series (LNIP, volume 12904)

Normal-ordering provides an approach to approximate three-body forces as effective two-body operators and it is therefore an important tool in many-body calculations with realistic nuclear interactions. The corresponding neglect of certain three-body terms in the normal-ordered Hamiltonian is known to influence translational invariance, although the magnitude of this effect has not yet been systematically quantified. In this work we study in particular the normal-ordering two-body approximation applied to a single harmonic-oscillator reference state. We explicate the breaking of translational invariance and demonstrate the magnitude of the approximation error as a function of model space parameters for $^4\rm{He}$ and $^{16}\rm{O}$ by performing full no-core shell-model calculations with and without three-nucleon forces. We combine two different diagnostics to better monitor the breaking of translational invariance. While the center-of-mass effect is shown to become potentially very large for $^4\rm{He}$, it is also shown to be much smaller for $^{16}\rm{O}$ although full convergence is not reached. These tools can be easily implemented in studies using other many-body frameworks and bases., Comment: 13 pages, 10 figures

The energy spectrum of electrons produced in molecular gas by interstellar cosmic rays (CRs) is rigorously calculated as a function of gas column density $N$ traversed by the CRs. This allows us to accurately compute the local value of the secondary ionization rate of molecular hydrogen, $\zeta_{\rm sec}(N)$, as a function of the local primary ionization rate, $\zeta_p(N)$. The ratio $\zeta_{\rm sec}/\zeta_p$ increases monotonically with $N$, and can considerably exceed the value of $\approx0.67$ commonly adopted in the literature. For sufficiently soft interstellar spectra, the dependence $\zeta_{\rm sec}/\zeta_p$ versus $N$ is practically insensitive to their particular shape and thus is a general characteristic of the secondary CR ionization in dense gas., Comment: Accepted to ApJ

We present long-duration numerical simulations of the tidal disruption of stars modelled with accurate stellar structures and spanning a range of pericentre distances, corresponding to cases where the stars are partially and completely disrupted. We substantiate the prediction that the late-time power-law index of the fallback rate $n_{\infty} \simeq -5/3$ for full disruptions, while for partial disruptions---in which the central part of the star survives the tidal encounter intact---we show that $n_{\infty} \simeq -9/4$. For the subset of simulations where the pericenter distance is close to that which delineates full from partial disruption, we find that a stellar core can reform after the star has been completely destroyed; for these events the energy of the zombie core is slightly positive, which results in late-time evolution from $n \simeq -9/4$ to $n \simeq -5/3$. We find that self-gravity can generate an $n(t)$ that deviates from $n_{\infty}$ by a small but significant amount for several years post-disruption. In one specific case with the stellar pericenter near the critical value, we find self-gravity also drives the re-collapse of the central regions of the debris stream into a collection of several cores while the rest of the stream remains relatively smooth. We also show that it is possible for the surviving stellar core in a partial disruption to acquire a circumstellar disc that is shed from the rapidly rotating core. Finally, we provide a novel analytical fitting function for the fallback rates that may also be useful in a range of contexts beyond TDEs., Comment: 21 pages, 12 figures, accepted for publication in ApJ

We introduce the Voronoi fundamental zone octonion interpolation framework for grain boundary (GB) structure-property models and surrogates. The VFZO framework offers an advantage over other five degree-of-freedom based property interpolation methods because it is constructed as a point set in a manifold. This means that directly computed Euclidean distances approximate the original octonion distance with significantly reduced computation runtime (~7 CPU minutes vs. 153 CPU days for a 50000x50000 pairwise-distance matrix). This increased efficiency facilitates lower interpolation error through the use of significantly more input data. We demonstrate grain boundary energy interpolation results for a non-smooth validation function and simulated bi-crystal datasets for Fe and Ni using four interpolation methods: barycentric interpolation, Gaussian process regression (GPR), inverse-distance weighting, and nearest-neighbor interpolation. These are evaluated for 50000 random input GBs and 10 000 random prediction GBs. The best performance was achieved with GPR, which resulted in a reduction of the root mean square error (RMSE) by 83.0% relative to RMSE of a constant, average model. Likewise, interpolation on a large, noisy, molecular statics Fe simulation dataset improves performance by 34.4% compared to 21.2% in prior work. Interpolation on a small, low-noise MS Ni simulation dataset is similar to interpolation results for the original octonion metric (57.6% vs. 56.4%). A vectorized, parallelized, MATLAB interpolation function (interp5DOF.m) and related routines are available in our VFZO repository (github.com/sgbaird-5dof/interp) which can be applied to other crystallographic point groups. The VFZO framework offers advantages for computing distances between GBs, estimating property values for arbitrary GBs, and modeling surrogates of computationally expensive 5DOF functions and simulations., Comment: main: 22 pages, 10 figures; appendices: 5 pages, 3 figures; supp: 13 pages, 12 figures

In this work we present the production of charged particles associated with high-$p_{\rm T}$ trigger particles ($8, Comment: 6 pages, 3 figures; to be published in the Proceedings of Science (PoS) of the 9th Edition of the Large Hadron Collider Physics Conference, 7-12 Jun 2021 (Online)

We present an energy scaling function to predict, in a specific range, the energy of bosonic trimers with large scattering lengths and finite range interactions, which is validated by quantum Monte Carlo calculations using microscopic Hamiltonians with two- and three-body potentials. The proposed scaling function depends on the scattering length, effective range, and a reference energy, which we chose as the trimer energy at unitarity. We obtained the scaling function as a limit cycle from the solution of the renormalized zero-range model with effective range corrections. We proposed a simple parametrization of the energy scaling function. Besides the intrinsic interest in theoretical and experimental investigations, this scaling function allows one to probe Efimov physics with only the trimer ground states, which may open opportunities to identify Efimov trimers whenever access to excited states is limited., Comment: 13 pages, 5 figures

A local UV cutoff $\Lambda(x)$ transforming under Weyl rescalings allows to construct Weyl invariant kinetic terms for scalar fields including Wilsonian cutoff functions. First we consider scalar fields in curved space-time with local bare couplings of any canonical dimension, and anomalous dimensions which describe their dependence on the UV cutoff. The local component of the UV cutoff plays the role of an additional coupling, albeit with a trivial constant $\beta$ function. This approach allows to derive Weyl consistency conditions for the corresponding anomalous dimensions which assume the form of an exact gradient flow. For renormalizable theories the Weyl consistency conditions are initially of the form of an approximate gradient flow for the $\beta$ functions, and we derive conditions under which it becomes the form of an exact gradient flow., Comment: 22 pages, Sections 3 and 4 re-written

International audience; The shape assumed by a slender elastic structure is a function both of the geometry of the space in which it exists and the forces it experiences. We explore, by experiments and theoretical analysis, the morphological phase space of a filament confined to the surface of a spherical bubble. The morphology is controlled by varying bending stiffness and weight of the filament, and its length relative to the bubble radius. When the dominant considerations are the geometry of confinement and elastic energy, the filament lies along a geodesic and when gravitational energy becomes significant, a bifurcation occurs, with a part of the filament occupying a longitude and the rest along a curve approximated by a latitude. Far beyond the transition, when the filament is much longer than the diameter, it coils around the selected latitudinal region. A simple model with filament shape as a composite of two arcs captures the transition well. For better quantitative agreement with the subcritical nature of bifurcation, we study the morphology by numerical energy minimization. Our analysis of the filament’s morphological space spanned by one geometric parameter, and one parameter that compares elastic energy with body forces, may provide guidance for packing slender structures on complex surfaces.

In this work, a novel multifidelity machine learning (ML) model, the gradient-enhanced multifidelity neural networks (GEMFNNs), is proposed. This model is a multifidelity version of gradient-enhanced neural networks (GENNs) as it uses both function and gradient information available at multiple levels of fidelity to make function approximations. Its construction is similar to multifidelity neural networks (MFNNs). This model is tested on three analytical function, a one, two, and a 20 variable function. It is also compared to neural networks (NNs), GENNs, and MFNNs, and the number of samples required to reach a global accuracy of 0.99 coefficient of determination (R^2) is measured. GEMFNNs required 18, 120, and 600 high-fidelity samples for the one, two, and 20 dimensional cases, respectively, to meet the target accuracy. NNs performed best on the one variable case, requiring only ten samples, while GENNs worked best on the two variable case, requiring 120 samples. GEMFNNs worked best for the 20 variable case, while requiring nearly eight times fewer samples than its nearest competitor, GENNs. For this case, NNs and MFNNs did not reach the target global accuracy even after using 10,000 high-fidelity samples. This work demonstrates the benefits of using gradient as well as multifidelity information in NNs for high-dimensional problems., Comment: Submitted to AMSE 2021 IDETC/CIE conference. 10 pages, 5 figures

We investigate how the properties of inhomogeneous patterns of activity, appearing in many natural and social phenomena, depend on the temporal resolution used to define individual bursts of activity. To this end, we consider time series of microscopic events produced by a self-exciting Hawkes process, and leverage a percolation framework to study the formation of macroscopic bursts of activity as a function of the resolution parameter. We find that the very same process may result in different distributions of avalanche size and duration, which are understood in terms of the competition between the 1D percolation and the branching process universality class. Pure regimes for the individual classes are observed at specific values of the resolution parameter corresponding to the critical points of the percolation diagram. A regime of crossover characterized by a mixture of the two universal behaviors is observed in a wide region of the diagram. The hybrid scaling appears to be a likely outcome for an analysis of the time series based on a reasonably chosen, but not precisely adjusted, value of the resolution parameter., Comment: 5 pages , 3 figures + Supplemental Material