Search

Analysis-suitable T-splines (AST-splines) are a promising candidate to achieve a seamless integration between the design and the analysis of thin-walled structures in industrial settings. In this work, we allow multiple extraordinary points per face, i.e., we remove the restriction of preceding works that required extraordinary points to be at least four rings apart from each other. We do so by mathematically showing that AST-splines with multiple extraordinary points per face are linearly independent and their polynomial basis functions form a non-negative partition of unity. This extension of the subset of AST-splines drastically increases the flexibility to build geometries using AST-splines; e.g., much coarser meshes can be constructed around small holes. The AST-spline spaces detailed in this work have C-1 inter-element continuity near extraordinary points and C-2 inter-element continuity elsewhere. For the convergence studies performed in this paper involving second-and fourth-order linear elliptic problems with manufactured solutions, we have not found any drawback caused by allowing multiple EPs per face in either the first refinement levels or the asymptotic behavior. To illustrate a possible isogeometric framework that is already available, we design the B-pillar and the side outer panel of a car using T-splines with the commercial software Autodesk Fusion360, import the control nets into our in-house code to build AST-splines, and import the Bezier extraction information into the commercial software LS-DYNA to solve eigenvalue problems. The results are compared with conventional finite elements and good agreement is found between AST-splines and conventional finite elements. (C)& nbsp;2021 Elsevier B.V. All rights reserved.

A functional k-batch code of dimension $s$ consists of $n$ servers storing linear combinations of $s$ linearly independent information bits. Any multiset request of size $k$ of linear combinations (or requests) of the information bits can be recovered by $k$ disjoint subsets of the servers. The goal under this paradigm is to find the minimum number of servers for given values of $s$ and $k$ . A recent conjecture states that for any $k=2^{s-1}$ requests the optimal solution requires $2^{s}-1$ servers. This conjecture is verified for $s\leqslant 5$ but previous work could only show that codes with $n=2^{s}-1$ servers can support a solution for $k=2^{s-2}+2^{s-4}+ \left\lfloor\frac{2^{s/2}}{\sqrt{24}} \right\rfloor$ requests. This paper reduces this gap and shows the existence of codes for $k= \lfloor\frac{2}{3}2^{s-1}\rfloor$ requests with the same number of servers. Another construction in the paper provides a code with $n=2^{s+1}-2$ servers and $k=2^{s}$ requests, which is an optimal result. These constructions are mainly based on Hadamard codes and equivalently provide constructions for parallel Random I/O (RIO) codes.

Consider a finite collection $\{T_1, \ldots, T_J\}$ of differential operators with constant coefficients on $\mathbb{T}^n$ ($n\geq 2$) and the space of smooth functions generated by this collection, namely, the space of functions $f$ such that $T_j f \in C(\mathbb{T}^n)$, $1\leq j\leq J$. We prove that if there are at least two linearly independent operators among their senior parts (relative to some mixed pattern of homogeneity), then this space does not have local unconditional structure. This fact generalizes the previously known result that such spaces are not isomorphic to a complemented subspace of $C(S)$, Comment: 22 pages

We study the storage capacity of quantum neural networks (QNNs), described by completely positive trace preserving (CPTP) maps acting on an N-dimensional Hilbert space. We demonstrate that attractor QNNs can store in a non-trivial manner up to N linearly independent pure states. For n qubits, QNNs can reach an exponential storage capacity, O ( 2 n ) , clearly outperforming standard classical neural networks whose storage capacity scales linearly with the number of neurons n. We estimate, employing the Gardner program, the relative volume of CPTP maps with M ⩽ N stationary states and show that this volume decreases exponentially with M and shrinks to zero for M ⩾ N + 1. We generalize our results to QNNs storing mixed states as well as input–output relations for feed-forward QNNs. Our approach opens the path to relate storage properties of QNNs to the quantum features of the input–output states. This paper is dedicated to the memory of Peter Wittek.

We prove that these Cantor sets are made up of transcendental numbers, apart from their endpoints $0$ and $1$, under some arithmetical assumptions on the data. To that purpose, we establish a criterion of linear independence over the field of algebraic numbers for the three numbers $1$, a characteristic Sturmian number, and an arbitrary Sturmian number with the same slope.

We investigate the ghostfree scalar-tensor theory with a timelike scalar field, with derivatives of the scalar field up to the third order and with the Riemann tensor up to the quadratic order. We build two types of linear spaces. One is the set of linearly independent generally covariant scalar-tensor monomials, the other is the set of linearly independent spatially covariant gravity monomials. We argue that these two types of linear space are isomorphic to each other in the sense of gauge fixing/recovering procedures. We then identify the subspaces in the spatially covariant gravity, which are spanned by linearly independent monomials built of the extrinsic and intrinsic curvature, the lapse function as well as their spatial derivatives, up to the fourth order in the total number of derivatives. The vectors in these subspaces, i.e., spatially covariant polynomials, automatically propagate at most three degrees of freedom. As a result, their images under the gauge recovering mappings are automatically the subspaces of scalar-tensor theory that propagate up to three degrees of freedom as long as the scalar field is timelike. The mappings from the spaces of spatially covariant gravity to the spaces of scalar-tensor theory are encoded in the projection matrices, of which we also derived the expressions explicitly. Our formalism and results can be useful in deriving the generally covariant higher derivative scalar-tensor theory without ghost(s)., Comment: 20 pages, 1 figure

Hermite interpolation property is desired in applied and computational mathematics. Hermite and vector subdivision schemes are of interest in CAGD for generating subdivision curves and in computational mathematics for building Hermite wavelets to numerically solve partial differential equations. In contrast to well-studied scalar subdivision schemes, Hermite and vector subdivision schemes employ matrix-valued masks and vector input data, which make their analysis much more complicated and difficult than their scalar counterparts. Under the spectral condition or the spectral chain, analysis of Hermite subdivision schemes through factorization of matrix-valued masks has been extensively studied in the literature and sufficient conditions have been given for the convergence of Hermite subdivision schemes through the contractivity of their derived subdivision schemes. We contribute to the study of Hermite subdivision schemes from a different perspective by investigating vector subdivision operators acting on vector polynomials and by establishing connections among Hermite subdivision schemes, vector cascade algorithms, and refinable vector functions. This approach allows us to characterize and construct all masks for Hermite subdivision schemes, to explain the spectral condition and spectral chain in the literature, to characterize convergence and smoothness of Hermite subdivision schemes using vector cascade algorithms, and to provide simple factorizations of Hermite masks through the normal form of matrix-valued masks such that the Hermite subdivision scheme is convergent if and only if its derived subdivision scheme is contractive. We also constructively prove that there always exist arbitrarily smooth convergent Hermite subdivision schemes, whose basis vector functions are splines and have linearly independent shifts. Several examples of Hermite subdivision schemes with short support and high smoothness are presented to illustrate the results in this paper.

Let $L$ be a lattice of full rank in $n$-dimensional real space. A vector in $L$ is called $i$-sparse if it has no more than $i$ nonzero coordinates. We define the $i$-th successive sparsity level of $L$, $s_i(L)$, to be the minimal $s$ so that $L$ has $s$ linearly independent $i$-sparse vectors, then $s_i(L) \leq n$ for each $1 \leq i \leq n$. We investigate sufficient conditions for $s_i(L)$ to be smaller than $n$ and obtain explicit bounds on the sup-norms of the corresponding linearly independent sparse vectors in~$L$. This result can be viewed as a partial sparse analogue of Minkowski's successive minima theorem. We then use this result to study virtually rectangular lattices, establishing conditions for the lattice to be virtually rectangular and determining the index of a rectangular sublattice. We further investigate the $2$-dimensional situation, showing that virtually rectangular lattices in the plane correspond to elliptic curves isogenous to those with real $j$-invariant. We also identify planar virtually rectangular lattices in terms of a natural rationality condition of the geodesics on the modular curve carrying the corresponding points., Comment: 18 pages; to appear in Research in the Mathematical Sciences

We present the new version 2.0 of the Feynman integral reduction program Kira and describe the new features. The primary new feature is the reconstruction of the final coefficients in integration-by-parts reductions by means of finite field methods with the help of FireFly. This procedure can be parallelized on computer clusters with MPI. Furthermore, the support for user-provided systems of equations has been significantly improved. This mode provides the flexibility to integrate Kira into projects that employ specialized reduction formulas, direct reduction of amplitudes, or to problems involving linear system of equations not limited to relations among standard Feynman integrals. We show examples from state-of-the-art Feynman integral reduction problems and provide benchmarks of the new features, demonstrating significantly reduced main memory usage and improved performance w.r.t. previous versions of Kira., Comment: 32 pages, v2: minor changes and clarifications, matches published version except for another fixed typo in Eq. (7)

Prethermalization, by introducing emergent quasiconserved observables, plays a crucial role in protecting Floquet many-body phases over exponentially long time, while the ultimate fate of such quasiconserved operators can signal thermalization to infinite temperature. To elucidate the properties of prethermal quasiconservation in many-body Floquet systems, here we systematically analyze infinite temperature correlations between observables. We numerically show that the late-time behavior of the autocorrelations unambiguously distinguishes quasiconserved observables from non-conserved ones, allowing to single out a set of linearly-independent quasiconserved observables. By investigating two Floquet spin models, we identify two different mechanism underlying the quasi-conservation law. First, we numerically verify energy quasiconservation when the driving frequency is large, so that the system dynamics is approximately described by a static prethermal Hamiltonian. More interestingly, under moderate driving frequency, another quasiconserved observable can still persist if the Floquet driving contains a large global rotation. We show theoretically how to calculate this conserved observable and provide numerical verification. Having systematically identified all quasiconserved observables, we can finally investigate their behavior in the infinite-time limit and thermodynamic limit, using autocorrelations obtained from both numerical simulation and experiments in solid state nuclear magnetic resonance systems., Comment: 12 pages, 9 figures. arXiv admin note: substantial text overlap with arXiv:1912.05799

Quasicrystals can be described as projections of sections of higher dimensional periodic lattices into real space. The image of the lattice points in the projected out dimensions, called the perpendicular space, carries valuable information about the local structure of the real space lattice. In this paper, we show that perpendicular space projection can be used to analyze the elementary excitations of a quasicrystal. In particular, we consider the vertex tight binding model on the two dimensional Penrose lattice and investigate the properties of strictly localized states using perpendicular space images. Our method reproduces the previously reported frequencies for the six types of localized states in this model. We also calculate the overlaps between different localized states and show that the number of type-5 and type-6 localized states which are independent from the four other types is a factor of golden ratio $\tau=(1+\sqrt{5})/2$ higher than previously reported values. Two orientations of the same type-5 or type-6 which are supported around the same site are shown to be linearly dependent with the addition of other types. We also show through exhaustion of all lattice sites in perpendicular space that any point in the Penrose lattice is either in the support of at least one localized state or is forbidden by local geometry to host a strictly localized state., Comment: 28 Pages, 39 Figures

Performing active quantum error correction to protect fragile quantum states highly depends on the correctness of error information--error syndromes. To obtain reliable error syndromes using imperfect physical circuits, we propose the idea of quantum data-syndrome (DS) codes that are capable of correcting both data qubits and syndrome bits errors. We study fundamental properties of quantum DS codes, including split weight enumerators, generalized MacWilliams identities, and linear programming bounds. In particular, we derive Singleton and Hamming-type upper bounds on degenerate quantum DS codes. Then we study random DS codes and show that random DS codes with a relatively small additional syndrome measurements achieve the Gilbert-Varshamov bound of stabilizer codes. Constructions of quantum DS codes are also discussed. A family of quantum DS codes is based on classical linear block codes, called syndrome measurement codes, so that syndrome bits are encoded in additional redundant stabilizer measurements. Another family of quantum DS codes is CSS-type quantum DS codes based on classical cyclic codes, and this includes the Steane code and the quantum Golay code., Comment: 18 pages, 4 figures. This version comprehensively extends our previous paper (arXiv:1602.01545)

We investigate the quantum state discrimination task for sets of linear independent pure states with an intrinsic ordering. This structured discrimination problems allow for a novel scheme that provides a certified level of error, that is, answers that never deviate from the true value more than a specified distance and hence a control of the desired quality of the results. We obtain an efficient semidefinite program and also find a general lower bound valid for any error distance that only requires the knowledge of optimal minimum error scheme. We apply our results to the quantum change point and quantum state anomaly detection cases., Comment: 9 pages, 10 figures

For every associative algebra $A$ and every class $\mathcal{C}$ of representations of $A$ the following question (related to nullstellensatz) makes sense: Characterize all tuples of elements $a_1,\ldots,a_n \in A$ such that vectors $\pi(a_1)v,\ldots,\pi(a_n)v$ are linearly dependent for every $\pi \in \mathcal{C}$ and every $v$ from the representation space of $\pi$. We answer this question in the following cases: (1) $A=U(L)$ is the enveloping algebra of a finite-dimensional complex Lie algebra $L$ and $\mathbb{C}$ is the class of all finite-dimensional representations of $A$. (2) $A=U(\mathfrak{sl}_2(\mathbb{C}))$ and $\mathbb{C}$ is the class of all finite-dimensional irreducible representations of $A$. (3) $A=U(\mathfrak{sl}_3(\mathbb{C}))$ and $\mathbb{C}$ is the class of all finite-dimensional irreducible representations of $A$ with sufficiently high weights. In case (1) the answer is: tuples that are linearly dependent over $\mathbb{C}$ while in cases (2) and (3) the answer is: tuples that are linearly dependent over the center of $A$. Similar results have been proved before for free algebras and Weyl algebras., Comment: 21 pages

Trirefringence is related to the existence of three distinct phase velocity solutions (and polarizations) for light propagation in a same wave-vector direction. This implies that when a trirefringent medium refracts a light ray, it is split into three rays with different velocities and linearly independent polarizations. Here a previous investigation [1] is revisited and its results are generalized to include a broader class of magnetoelectric materials. Moreover, it is argued that trirefringent media could already be achieved using present day technologies. Examples are given to support that focused on built layered media under the influence of applied electromagnetic fields., Comment: 9 pages, 2 figures

We investigate quantum information masking for arbitrary dimensional quantum states. We show that mutually orthogonal quantum states can always be served for deterministic masking of quantum information. We further construct a probabilistic masking machine for linearly independent states. It is shown that a set of d dimensional states, $\{ |a_1 \rangle_A, |t a_2 \rangle_A, \dots, |a_n \rangle_A \}$, $n \leq d$, can be probabilistically masked by a general unitary-reduction operation if they are linearly independent. The maximal successful probability of probabilistic masking is analyzed and derived for the case of two initial states., Comment: 5 pages, 1 figures

Interferometric visibilities, reflecting the complex correlations between signals recorded at antennas in an interferometric array, carry information about the angular structure of a distant source. While unknown antenna gains in both amplitude and phase can prevent direct interpretation of these measurements, certain combinations of visibilities called closure phases and closure amplitudes are independent of antenna gains and provide a convenient set of robust observables. However, these closure quantities have subtle noise properties and are generally both linearly and statistically dependent. These complications have obstructed the proper use of closure quantities in interferometric analysis, and they have obscured the relationship between analysis with closure quantities and other analysis techniques such as self calibration. We review the statistics of closure quantities, noting common pitfalls that arise when approaching low signal-to-noise due to the nonlinear propagation of statistical errors. We then develop a strategy for isolating and fitting to the independent degrees of freedom captured by the closure quantities through explicit construction of linearly independent sets of quantities along with their noise covariance in the Gaussian limit, valid for moderate signal-to-noise, and we demonstrate that model fits have biased posteriors when this covariance is ignored. Finally, we introduce a unified procedure for fitting to both closure information and partially calibrated visibilities, and we demonstrate both analytically and numerically the direct equivalence of inference based on closure quantities to that based on self calibration of complex visibilities with unconstrained antenna gains., Comment: 31 pages, 17 figures

Many rotational invariants for crystal structure representations have been used to describe the structure-property relationship by machine learning. The machine learning interatomic potential (MLIP) is one of the applications of rotational invariants, which provides the relationship between the energy and the crystal structure. Since the MLIP requires the highest accuracy among machine learning estimations of the structure-property relationship, the enumeration of rotational invariants is useful for constructing MLIPs with the desired accuracy. In this study, we introduce high-order linearly independent rotational invariants up to the sixth order based on spherical harmonics and apply them to linearized MLIPs for elemental aluminum. A set of rotational invariants is derived by the general process of reducing the Kronecker products of irreducible representations (Irreps) for the SO(3) group using a group-theoretical projector method. A high predictive power for a wide range of structures is accomplished by using high-order invariants with low-order invariants equivalent to pair and angular structural features., Comment: 15 pages, 4 figures

This paper proposes a new evaluation metric and boosting method for weight separability in neural network design. In contrast to general visual recognition methods designed to encourage both intra-class compactness and inter-class separability of latent features, we focus on estimating linear independence of column vectors in weight matrix and improving the separability of weight vectors. To this end, we propose an evaluation metric for weight separability based on semi-orthogonality of a matrix and Frobenius distance, and the feed-backward reconstruction loss which explicitly encourages weight separability between the column vectors in the weight matrix. The experimental results on image classification and face recognition demonstrate that the weight separability boosting via minimization of feed-backward reconstruction loss can improve the visual recognition performance, hence universally boosting the performance on various visual recognition tasks., Comment: 8 pages, 6 figures

Let $R_k(x)$ denote the error incurred by approximating the number of $k$-free integers less than $x$ by $x/\zeta(k)$. It is well known that $R_k(x)=\Omega(x^{\frac{1}{2k}})$, and widely conjectured that $R_k(x)=O(x^{\frac{1}{2k}+\epsilon})$. By establishing weak linear independence of some subsets of zeros of the Riemann zeta function, we establish an effective proof of the lower bound, with significantly larger bounds on the constant compared to those obtained in prior work. For example, we show that $R_k(x)/x^{1/2k} > 3$ infinitely often and that $R_k(x)/x^{1/2k} < -3$ infinitely often, for $k=2$, $3$, $4$, and $5$. We also investigate $R_2(x)$ and $R_3(x)$ in detail and establish that our bounds far exceed the oscillations exhibited by these functions over a long range: for $0

Mathematical programs with disjunctive constraints (MPDCs for short) cover several different problem classes from nonlinear optimization including complementarity-, vanishing-, cardinality-, and switching-constrained optimization problems. In this paper, we introduce an abstract but reasonable version of the prominent linear independence constraint qualification which applies to MPDCs. Afterwards, we derive first- and second-order optimality conditions for MPDCs under validity of this constraint qualification based on so-called strongly stationary points. Finally, we apply our findings to some popular classes of disjunctive programs and compare the obtained results to those ones available in the literature. Particularly, new second-order optimality conditions for mathematical programs with switching constraints are by-products of our approach., Comment: 35 pages, 3 figures

For the same potential as originally studied by Ma [Phys. Rev. {\bf 71}, 195 (1947)] we obtain analytic expressions for the Jost functions and the residui of the S-matrix of both (i) redundant poles and (ii) the poles corresponding to true bound states. This enables us to demonstrate that the Heisenberg condition is valid in spite of the presence of redundant poles and singular behaviour of the S-matrix for $k\to \infty$. In addition, we analytically determine the overall contribution of redundant poles to the asymptotic completeness relation, provided that the residuum theorem can be applied. The origin of redundant poles and zeros is shown to be related to peculiarities of analytic continuation of a parameter of two linearly independent analytic functions., Comment: The Heisenberg condition has been reduced down to analytic relation (35) for the cylindrical Bessel functions (7 pages, 1 figure)

A functional k-Private Information Retrieval (k -PIR) code of dimension s consists of $n$ servers storing linear combinations of s linearly independent information symbols. Any linear combination of the s information symbols can be recovered by k disjoint subsets of servers. The goal is to find the minimum number of servers for given k and s . We provide lower bounds on the minimum number of servers and constructions which yield upper bounds on this number. For k ≤ 4, exact bounds on this number are proved. Furthermore, we provide some asymptotic bounds. The problem coincides with the well known PIR problem based on a coded database to reduce the storage overhead, when each linear combination contains exactly one information symbol. If any multiset of size k of linear combinations from the linearly independent information symbols can be recovered by k disjoint subset of servers, then the servers form a functional k -batch code . A functional k -batch code is a functional k -PIR code, where all the k linear combinations in the multiset are equal. We provide some bounds on the minimum number of servers for functional k -batch codes. In particular we present a random construction and a construction based on simplex codes, Write-Once Memory (WOM) codes, and Random I/O (RIO) codes.

Let $n$ be a positive integer and $I$ a $k$-subset of integers in $[0,n-1]$. Given a $k$-tuple $A=(\alpha_0, \cdots, \alpha_{k-1})\in \mathbb{F}^k_{q^n}$, let $M_{A,I}$ denote the matrix $(\alpha_i^{q^j})$ with $0\leq i\leq k-1$ and $j\in I$. When $I=\{0,1,\cdots, k-1\}$, $M_{A,I}$ is called a Moore matrix which was introduced by E. H. Moore in 1896. It is well known that the determinant of a Moore matrix equals $0$ if and only if $\alpha_0,\cdots, \alpha_{k-1}$ are $\mathbb{F}_q$-linearly dependent. We call $I$ that satisfies this property a Moore exponent set. In fact, Moore exponent sets are equivalent to maximum rank-distance (MRD) code with maximum left and right idealisers over finite fields. It is already known that $I=\{0,\cdots, k-1\}$ is not the unique Moore exponent set, for instance, (generalized) Delsarte-Gabidulin codes and the MRD codes recently discovered by Csajb\'ok, Marino, Polverino and the second author both give rise to new Moore exponent sets. By using algebraic geometry approach, we obtain an asymptotic classification result: for $q>5$, if $I$ is not an arithmetic progression, then there exist an integer $N$ depending on $I$ such that $I$ is not a Moore exponent set provided that $n>N$., Comment: 15 pages,typos are corrected, to appear in JCTA

A simple binary matroid is called claw-free if none of its rank-3 flats are independent sets. These objects can be equivalently defined as the sets E of points in PG ( n − 1 , 2 ) for which E ∩ P is not a basis of P for any plane P, or as the subsets X of F 2 n containing no linearly independent triple x , y , z for which x + y , y + z , x + z , x + y + z ∉ X . We prove a decomposition theorem that exactly determines the structure of all claw-free matroids. The theorem states that claw-free matroids either belong to one of three particular basic classes of claw-free matroids, or can be constructed from these basic classes using a certain ‘join’ operation.

We consider sets of trace-normalized non-negative operators in Hilbert-Schmidt balls that maximize their mutual Hilbert-Schmidt distance; these are optimal arrangements in the sets of purity-limited classical or quantum states on a finite-dimensional Hilbert space. Classical states are understood to be represented by diagonal matrices, with the diagonal entries forming a probability vector. We also introduce the concept of spectrahedron arrangements which provides a unified framework for classical and quantum arrangements and the flexibility to define new types of optimal packings. Continuing a prior work, we combine combinatorial structures and line packings associated with frames to arrive at optimal arrangements of higher-rank quantum states. One new construction that is presented involves generating an optimal arrangement we call a Gabor-Steiner equiangular tight frame as the orbit of a projective representation of the Weyl-Heisenberg group over any finite abelian group. The minimal sets of linearly dependent vectors, the so-called binder, of the Gabor-Steiner equiangular tight frames are then characterized; under certain conditions these form combinatorial block designs and in one case generate a new class of block designs. The projections onto the span of minimal linearly dependent sets in the Gabor-Steiner equiangular tight frame are then used to generate further optimal spectrahedron arrangements., Comment: 39 pages

Let $E$ be an elliptic curve over $\mathbb{Q}$ described by $y^2= x^3+ Kx+ L$ where $K, L \in \mathbb{Q}$. A set of rational points $(x_i,y_i) \in E(\mathbb{Q})$ for $i=1, 2, \cdots, k$, is said to be a sequence of consecutive cubes on $E$ if the $x-$coordinates of the points $x_i$'s for $i=1, 2, \cdots$ form consecutive cubes. In this note, we show the existence of an infinite family of elliptic curves containing a length-$5$-term sequence of consecutive cubes. Morever, these five rational points in $E (\mathbb{Q})$ are linearly independent and the rank $r$ of $E(\mathbb{Q})$ is at least $5$., Comment: 10 pages, to appear, Rocky Mountain Journal of Mathematics (2018)

Unsupervised machine learning methods are used to identify structural changes using the melting point transition in classical molecular dynamics simulations as an example application of the approach. Dimensionality reduction and clustering methods are applied to instantaneous radial distributions of atomic configurations from classical molecular dynamics simulations of metallic systems over a large temperature range. Principal component analysis is used to dramatically reduce the dimensionality of the feature space across the samples using an orthogonal linear transformation that preserves the statistical variance of the data under the condition that the new feature space is linearly independent. From there, $k$-means clustering is used to partition the samples into solid and liquid phases through a criterion motivated by the geometry of the reduced feature space of the samples, allowing for an estimation of the melting point transition. This pattern criterion is conceptually similar to how humans interpret the data but with far greater throughput, as the shapes of the radial distributions are different for each phase and easily distinguishable by humans. The transition temperature estimates derived from this machine learning approach produce comparable results to other methods on similarly small system sizes. These results show that machine learning approaches can be applied to structural changes in physical systems.

Rubinstein and Sarnak have shown, conditional on the Riemann hypothesis (RH) and the linear independence hypothesis (LI) on the non-real zeros of $\zeta(s)$, that the set of real numbers $x\ge2$ for which $\pi(x)>$ li$(x)$ has a logarithmic density, which they computed to be about $2.6\times10^{-7}$. A natural problem is to examine the actual primes in this race. We prove, assuming RH and LI, that the logarithmic density of the set of primes $p$ for which $\pi(p)>$ li$(p)$ relative to the prime numbers exists and is the same as the Rubinstein-Sarnak density. We also extend such results to a broad class of prime number races, including the "Mertens race" between $\prod_{p< x}(1-1/p)^{-1}$ and $e^{\gamma}\log x$ and the "Zhang race" between $\sum_{p\ge x}1/(p\log p)$ and $1/\log x$. These latter results resolve a question of the first and third author from a previous paper, leading to further progress on a 1988 conjecture of Erd\H{o}s on primitive sets., Comment: 14 pages

New techniques based on weak measurements have recently been introduced to the field of quantum state reconstruction. Some of them allow the direct measurement of each matrix element of an unknown density operator and need only $O(d)$ different operations, compared to $d^2$ linearly independent projectors in the case of standard quantum state tomography, for the reconstruction of an arbitrary mixed state. However, due to the weakness of these couplings, these protocols are approximated and prone to large statistical errors. We propose a method which is similar to the weak measurement protocols but works regardless of the coupling strength: our protocol is not approximated and thus improves the accuracy and precision of the results with respect to weak measurement schemes. We experimentally apply it to the polarization state of single photons and compare the results to those of preexisting methods for different values of the coupling strength. Our results show that our method outperforms previous proposals in terms of accuracy and statistical errors., Comment: RevTex, 6 pages

Forking is a central notion of model theory, generalizing linear independence in vector spaces and algebraic independence in fields. We develop the theory of forking in abstract, category-theoretic terms, for reasons both practical (we require a characterization suitable for work in $\mu$-abstract elementary classes, i.e. accessible categories with all morphisms monomorphisms) and expository (we hope, with this account, to make forking accessible - and useful - to a broader mathematical audience). In particular, we present an axiomatic definition of what we call a stable independence notion on a category and show that this is in fact a purely category-theoretic axiomatization of the properties of model-theoretic forking in a stable first-order theory., Comment: 50 pages

We propose a scheme to enhance the fidelity of a symmetric quantum cloning machine via weak measurements. By adjusting the intensity of the weak measurement parameter $p$, we obtain copies of initial states with different values of fidelity. Choosing a proper value of $p$, we can obtain perfect copies. In this paper, we focus on $1\ensuremath{-}2$ quantum cloning for two nonorthogonal states. Sets containing more than two linear independent states are also briefly discussed. Due to the probabilistic nature of weak measurements, we obtain high-fidelity copies probabilistically. If the weak measurement is a success, we do subsequent operations to obtain high-fidelity copies; otherwise, the cloning process fails and we quit. From this perspective, the scheme we propose is economical for saving quantum resources and time. Our scheme may be very useful in quantum information processing.

A lattice is a set of all the integer linear combinations of certain linearly independent vectors. One of the most important concepts on lattice is the successive minima which is of vital importance from both theoretical and practical applications points of view. In this paper, we first study some properties of successive minima and then employ some of them to improve the suboptimal algorithm for solving an optimization problem about maximizing the achievable rate of the integer-forcing strategy for cloud radio access networks in [1]., Comment: to appear in ISIT 2018

This paper considers a low-complexity Gaussian Message Passing (GMP) scheme for a coded massive Multiple-Input Multiple-Output (MIMO) systems with Non-Orthogonal Multiple Access (massive MIMO-NOMA), in which a base station with $N_s$ antennas serves $N_u$ sources simultaneously in the same frequency. Both $N_u$ and $N_s$ are large numbers, and we consider the overloaded cases with $N_u>N_s$. The GMP for MIMO-NOMA is a message passing algorithm operating on a fully-connected loopy factor graph, which is well understood to fail to converge due to the correlation problem. In this paper, we utilize the large-scale property of the system to simplify the convergence analysis of the GMP under the overloaded condition. First, we prove that the \emph{variances} of the GMP definitely converge to the mean square error (MSE) of Linear Minimum Mean Square Error (LMMSE) multi-user detection. Secondly, the \emph{means} of the traditional GMP will fail to converge when $ N_u/N_s< (\sqrt{2}-1)^{-2}\approx5.83$. Therefore, we propose and derive a new convergent GMP called scale-and-add GMP (SA-GMP), which always converges to the LMMSE multi-user detection performance for any $N_u/N_s>1$, and show that it has a faster convergence speed than the traditional GMP with the same complexity. Finally, numerical results are provided to verify the validity and accuracy of the theoretical results presented., Accepted by IEEE TWC, 16 pages, 11 figures

A [well-known conjecture](https://en.wikipedia.org/wiki/Gy%C3%A1rf%C3%A1s%E2%80%93Sumner_conjecture) of [Andras Gyarfas](https://en.wikipedia.org/wiki/Andr%C3%A1s_Gy%C3%A1rf%C3%A1s) and [David Sumner](https://en.wikipedia.org/wiki/David_Sumner) states that for every positive integer $m$ and every finite tree $T$ there exists $k$ such that all graphs that do not contain the clique $K_m$ or an induced copy of $T$ have chromatic number at most $k$. The conjecture has been proved in many special cases, but the general case has been open for several decades. The main purpose of this paper is to consider a natural analogue of the conjecture for matroids, where it turns out, interestingly, to be false. Matroids are structures that result from abstracting the notion of independent sets in vector spaces: that is, a matroid is a set $M$ together with a nonempty hereditary collection $\mathcal I$ of subsets deemed to be independent where all maximal independent subsets of every set are equicardinal. They can also be regarded as generalizations of graphs, since if $G$ is any graph and $\mathcal I$ is the collection of all acyclic subsets of $E(G)$, then the pair $(E(G),\mathcal I)$ is a matroid. In fact, it is a _binary_ matroid, which means that it can be represented as a subset of a vector space over $\mathbb F_2$. To do this, we take the space of all formal sums of vertices and represent the edge $vw$ by the sum $v+w$. A set of edges is easily seen to be acyclic if and only if the corresponding set of sums is linearly independent. There is a natural analogue of an induced subgraph for matroids: an induced restriction of a matroid $M$ is a subset $M’$ of $M$ with the property that adding any element of $M-M’$ to $M’$ produces a matroid with a larger independent set than $M’$. The natural analogue of a tree with $m$ edges is the matroid $I_m$, where one takes a set of size $m$ and takes all its subsets to be independent. (Note, however, that unlike with graph-theoretic trees there is just one such matroid up to isomorphism for each $m$.) Every graph can be obtained by deleting edges from a complete graph. Analogously, every binary matroid can be obtained by deleting elements from a finite _binary projective geometry_, that is, the set of all one-dimensional subspaces in a finite-dimensional vector space over $\mathbb F_2$. Finally, the analogue of the chromatic number for binary matroids is a quantity known as the _critical number_ introduced by Crapo and Rota, which in the case of a graph $G$ turns out to be $\lceil\log_2(\chi(G))\rceil$ -- that is, roughly the logarithm of its chromatic number. One of the results of the paper is that a binary matroid can fail to contain $I_3$ or the Fano plane $F_7$ (which is the simplest projective geometry) as an induced restriction, but also have arbitrarily large critical number. By contrast, the critical number is at most two if one also excludes the matroid associated with $K_5$ as an induced restriction. The main result of the paper is a structural description of all simple binary matroids that have neither $I_3$ nor $F_7$ as an induced restriction.

We show that the bounded Borel class of any dense representation $\rho: G\to \PSL_n\bC$ is non-zero in degree three bounded cohomology and has maximal semi-norm, for any discrete group $G$. When $n=2$, the Borel class is equal to the $3$-dimensional hyperbolic volume class. Using tools from the theory of Kleinian groups, we show that the volume class of a dense representation $\rho: G\to \PSL_2\bC$ is uniformly separated in semi-norm from any other representation $\rho': G\to \G$ for which there is a subgroup $H\le G$ on which $\rho$ is still dense but $\rho'$ is discrete or indiscrete but stabilizes a point, line, or plane in $\bH^3\cup \partial \bH^3$. We exhibit a family of dense representations of a non-abelian free group on two letters and a family of discontinuous dense representations of $\PSL_2\bR$, whose volume classes are linearly independent and satisfy some additional properties; the cardinality of these families is that of the continuum. We explain how the strategy employed may be used to produce non-trivial volume classes in higher dimensions, contingent on the existence of a family of hyperbolic manifolds with certain topological and geometric properties., Comment: 36 pages, 1 figure; final version accepted for publication in IMRN

We present a mechanical proof of the Cauchy-Schwarz inequality in ACL2(r) and a formalisation of the necessary mathematics to undertake such a proof. This includes the formalisation of $\mathbb{R}^n$ as an inner product space. We also provide an application of Cauchy-Schwarz by formalising $\mathbb R^n$ as a metric space and exhibiting continuity for some simple functions $\mathbb R^n\to\mathbb R$. The Cauchy-Schwarz inequality relates the magnitude of a vector to its projection (or inner product) with another: \[|\langle u,v\rangle| \leq \|u\| \|v\|\] with equality iff the vectors are linearly dependent. It finds frequent use in many branches of mathematics including linear algebra, real analysis, functional analysis, probability, etc. Indeed, the inequality is considered to be among "The Hundred Greatest Theorems" and is listed in the "Formalizing 100 Theorems" project. To the best of our knowledge, our formalisation is the first published proof using ACL2(r) or any other first-order theorem prover., Comment: In Proceedings ACL2 2018, arXiv:1810.03762

Recently, path norm was proposed as a new capacity measure for neural networks with Rectified Linear Unit (ReLU) activation function, which takes the rescaling-invariant property of ReLU into account. It has been shown that the generalization error bound in terms of the path norm explains the empirical generalization behaviors of the ReLU neural networks better than that of other capacity measures. Moreover, optimization algorithms which take path norm as the regularization term to the loss function, like Path-SGD, have been shown to achieve better generalization performance. However, the path norm counts the values of all paths, and hence the capacity measure based on path norm could be improperly influenced by the dependency among different paths. It is also known that each path of a ReLU network can be represented by a small group of linearly independent basis paths with multiplication and division operation, which indicates that the generalization behavior of the network only depends on only a few basis paths. Motivated by this, we propose a new norm Basis-path Norm based on a group of linearly independent paths to measure the capacity of neural networks more accurately. We establish a generalization error bound based on this basis path norm, and show it explains the generalization behaviors of ReLU networks more accurately than previous capacity measures via extensive experiments. In addition, we develop optimization algorithms which minimize the empirical risk regularized by the basis-path norm. Our experiments on benchmark datasets demonstrate that the proposed regularization method achieves clearly better performance on the test set than the previous regularization approaches.

Given $k\in\mathbb N$, we study the vanishing of the Dirichlet series $$D_k(s,f):=\sum_{n\geq1} d_k(n)f(n)n^{-s}$$ at the point $s=1$, where $f$ is a periodic function modulo a prime $p$. We show that if $(k,p-1)=1$ or $(k,p-1)=2$ and $p\equiv 3\mod 4$, then there are no odd rational-valued functions $f\not\equiv 0$ such that $D_k(1,f)=0$, whereas in all other cases there are examples of odd functions $f$ such that $D_k(1,f)=0$. As a consequence, we obtain, for example, that the set of values $L(1,\chi)^2$, where $\chi$ ranges over odd characters mod $p$, are linearly independent over $\mathbb Q$., Comment: 16 pages

The aim of this paper is to study some properties of left translates of a square integrable function on the Heisenberg group. First, a necessary and sufficient condition for the existence of the canonical dual to a function $\varphi\in L^{2}(\mathbb{R}^{2n})$ is obtained in the case of twisted shift-invariant spaces. Further, characterizations of $\ell^{2}$-linear independence and the Hilbertian property of the twisted translates of a function $\varphi\in L^{2}(\mathbb{R}^{2n})$ are obtained. Later these results are shown in the case of the Heisenberg group., Comment: 13 pages

Let $S$ be a compact oriented surface. We construct homogeneous quasimorphisms on $Diff(S, area)$, on $Diff_0(S, area)$ and on $Ham(S)$ generalizing the constructions of Gambaudo-Ghys and Polterovich. We prove that there are infinitely many linearly independent homogeneous quasimorphisms on $Diff(S, area)$, on $Diff_0(S, area)$ and on $Ham(S)$ whose absolute values bound from below the topological entropy. In case when $S$ has a positive genus, the quasimorphisms we construct on $Ham(S)$ are $C^0$-continuous. We define a bi-invariant metric on these groups, called the entropy metric, and show that it is unbounded. In particular, we reprove the fact that the autonomous metric on $Ham(S)$ is unbounded., Comment: 23 pages, one figure. In this version the main results are proved for all compact oriented surfaces. To appear in Journal of Modern Dynamics

We show that for each Seifert form of an algebraically slice knot with nontrivial Alexander polynomial, there exists an infinite family of knots having the Seifert form such that the knots are linearly independent in the knot concordance group and not concordant to any knot with coprime Alexander polynomial. Key ingredients for the proof are Cheeger-Gromov-von Neumann rho-invariants for amenable groups developed by Cha and Orr and polynomial splittings of metabelian rho-invariants., Comment: 10 pages, 1 figure; typos corrected

Let $\Omega\subset \mathbb{R}^n$ be open and let $\mathcal{R}$ be a partial frame on $\Omega$, that is a set of $m$ linearly independent vector fields prescribed on $\Omega$ ($m\leq n$). We consider the issue of describing the set of all maps $F:\Omega\to\mathbb{R}^n$ with the property that each of the given vector fields is an eigenvector of the Jacobian matrix of $F$. By introducing a coordinate independent definition of the Jacobian, we obtain an intrinsic formulation of the problem, which leads to an overdetermined PDE system, whose compatibility conditions can be expressed in an intrinsic, coordinate independent manner. To analyze this system we formulate and prove a generalization of the classical Frobenius integrability theorems. The size and structure of the solution set of this system depends on the properties of the partial frame, in particular, whether or not it is in involution. A particularly nice subclass of involutive partial frames, called rich, can be completely analyzed. Involutive, but non-rich case is somewhat harder to handle. We provide a complete answer in the case of $m=3$ and arbitrary $n$, as well as some general results for arbitrary $m$. The non-involutive case is far more challenging, and we only obtain a comprehensive analysis in the case $n=3$, $m=2$. Finally, we provide explicit examples illustrating the various possibilities. Our initial motivation for considering this problem comes from the geometric study of hyperbolic conservative systems in one spatial dimension.

We present an efficient approach of Finite Element Method (FEM)-based nonrigid image registration, in which the spatial transformation is constructed using truncated hierarchical B-splines (THB-splines). The image registration framework minimizes an energy functional using an FEM-based method and thus involves solving a large system of linear equations. This framework is carried out on a set of successively refined grids. However, due to the increased number of control points during subdivision, large linear systems are generated which are generally demanding to solve. Instead of using uniform subdivision, an adaptive local refinement scheme is carried out, only refining the areas of large change in deformation of the image. By incorporating the key advantages of THB-spline basis functions such as linear independence, partition of unity and reduced overlap into the FEM-based framework, we improve the matrix sparsity and computational efficiency. The performance of the proposed method is demonstrated on 2D synthetic and medical images.

Bourgain and Chang recently showed that any subset of $\mathbb{F}_p$ of density $\gg p^{-1/15}$ contains a nontrivial progression $x,x+y,x+y^2$. We answer a question of theirs by proving that if $P_1,P_2\in\mathbb{Z}[y]$ are linearly independent and satisfy $P_1(0)=P_2(0)=0$, then any subset of $\mathbb{F}_p$ of density $\gg_{P_1,P_2}p^{-1/24}$ contains a nontrivial polynomial progression $x,x+P_1(y),x+P_2(y)$., Comment: 21 pages; v2: minor exposition changes based on referee comments

About two decades ago, Tsfasman and Boguslavsky conjectured a formula for the maximum number of common zeros that $r$ linearly independent homogeneous polynomials of degree $d$ in $m+1$ variables with coefficients in a finite field with $q$ elements can have in the corresponding $m$-dimensional projective space. Recently, it has been shown by Datta and Ghorpade that this conjecture is valid if $r$ is at most $m+1$ and can be invalid otherwise. Moreover a new conjecture was proposed for many values of $r$ beyond $m+1$. In this paper, we prove that this new conjecture holds true for several values of $r$. In particular, this settles the new conjecture completely when $d=3$. Our result also includes the positive result of Datta and Ghorpade as a special case. Further, we determine the maximum number of zeros in certain cases not covered by the earlier conjectures and results, namely, the case of $d=q-1$ and of $d=q$. All these results are directly applicable to the determination of the maximum number of points on sections of Veronese varieties by linear subvarieties of a fixed dimension, and also the determination of generalized Hamming weights of projective Reed-Muller codes., Comment: 15 pages

In this article, we present a new implementation of the Laporta algorithm to reduce scalar multi-loop integrals---appearing in quantum field theoretic calculations---to a set of master integrals. We extend existing approaches by using an additional algorithm based on modular arithmetic to remove linearly dependent equations from the system of equations arising from integration-by-parts and Lorentz identities. Furthermore, the algebraic manipulations required in the back substitution are optimized. We describe in detail the implementation as well as the usage of the program. In addition, we show benchmarks for concrete examples and compare the performance to Reduze 2 and FIRE 5. In our benchmarks we find that Kira is highly competitive with these existing tools., Comment: 37 pages, 3 figures

We prove that, under certain conditions on the function pair $\varphi_1$ and $\varphi_2$, bilinear average $p^{-1}\sum_{y\in \mathbb{F}_p}f_1(x+\varphi_1(y)) f_2(x+\varphi_2(y))$ along curve $(\varphi_1, \varphi_2)$ satisfies certain decay estimate. As a consequence, Roth type theorems hold in the setting of finite fields. In particular, if $\varphi_1,\varphi_2\in \mathbb{F}_p[X]$ with $\varphi_1(0)=\varphi_2(0)=0$ are linearly independent polynomials, then for any $A\subset \mathbb{F}_p, |A|=\delta p$ with $\delta>c p^{-\frac{1}{12}}$, there are $\gtrsim \delta^3p^2$ triplets $x,x+\varphi_1(y), x+\varphi_2(y)\in A$. This extends a recent result of Bourgain and Chang who initiated this type of problems, and strengthens the bound in a result of Peluse, who generalized Bourgain and Chang's work. The proof uses discrete Fourier analysis and algebraic geometry., Comment: The assumption "having distinct leading terms" is removed in the results and Peluse is acknowledged

The linear superposition principle in quantum mechanics is essential for several no-go theorems such as the no-cloning theorem, the no-deleting theorem and the no-superposing theorem. It remains an open problem of finding general forbidden principles to unify these results. In this paper, we investigate general quantum transformations forbidden or permitted by the superposition principle for various goals. First, we prove a no-encoding theorem that forbids linearly superposing of an unknown pure state and a fixed state in Hilbert space of finite dimension. Two general forms include the no-cloning theorem, the no-deleting theorem, and the no-superposing theorem as special cases. Second, we provide a unified scheme for presenting perfect and imperfect quantum tasks (cloning and deleting) in a one-shot manner. This scheme may yield to fruitful results that are completely characterized with the linear independence of the input pure states. The generalized upper bounds for the success probability will be proved. Third, we generalize a recent superposing of unknown states with fixed overlaps when multiple copies of the input states are available.

Condition for distinguishability of countably infinite number of pure states by a single measurement is given. Distinguishability is to be understood as possibility of an unambiguous measurement. For finite number of states, it is known that the necessary and sufficient condition of distinguishability is that the states are linearly independent. For infinite number of states, several natural classes of distinguishability can be defined. We give a necessary and sufficient condition for a system of pure states to be distinguishable. It turns out that each level of distinguishability naturally corresponds to one of the generalizations of linear independence to families of infinite vectors. As an important example, we apply the general theory to von Neumann's lattice, a subsystem of coherent states which corresponds to a lattice in the classical phase space. We prove that the condition for distinguishability is that the area of the fundamental region of the lattice is greater than the Planck constant, and also find subtle behavior on the threshold. These facts reveal the measurement theoretical meaning of the Planck constant and give a justification for the interpretation that it is the smallest unit of area in the phase space. The cases of uncountably many states and of mixed states are also discussed., Comment: 19 pages