Source author record

Paul B. Slater

Paul B. Slater appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
12topics
2close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2021arXiv

Hilbert-Schmidt Separability Probabilities from Bures Ensembles and vice versa: Applications to Quantum Steering Ellipsoids and Monotone Metrics

We reexamine a recent analysis in which, using the volume of the associated quantum steering ellipsoid (QES) as a measure, we sought to estimate the probability that a two-qubit state is separable. In the estimation process, we, in effect, sought to attach to states random with respect to Hilbert-Schmidt (HS) measure, the corresponding QES volumes. However, a study of the relations between HS and Bures ensembles and their well-supported separability probabilities of $\frac{8}{33}$ and $\frac{25}{341}$, respectively, now lead us to explore as a possible alternative measure, the QES volume divided by the $Π_{j<k}^{1...4} (λ_j-λ_k)^2 $ term of the HS volume element (the $λ$'s being the four eigenvalues of the associated $4 \times 4$ density matrix $ρ$). This measure is applied to the members of a HS ensemble of random two-qubit states, yielding a QES separability probability estimate of 0.105458. Alternatively, weighting members of a Bures ensemble by the QES volume divided by the eigenvalue part $\frac{1}{\sqrt{\mbox{det} ρ}} Π_{j<k}^{1...4} \frac{(λ_j-λ_k)^2}{λ_j+λ_k}$ of the Bures volume element, gives a close estimate of 0.100223. We also weight members of a HS ensemble by the QES volume divided not only by the indicated HS eigenvalue term, but also by the unitary component $|Π_{j<k}^{1...4} \mbox{Re} (U^{-1}) \mbox{Im} (U^{-1})|$ of the volume element. For one hundred thirty (rather variable) independent separability probability estimates, we, then, obtain median and mean estimates of 0.0447729 and 0.117485 with variance 0.0381468.

preprint2020arXiv

8 Boolean Atoms Spanning the 256-Dimensional Entanglement-Probability Three-Set Algebra of the Two-Qutrit Hiesmayr-Loffler Magic Simplex of Bell States

We obtain formulas (bot. p. 12)--including $\frac{2}{121}$ and $\frac{4 \left(242 \sqrt{3} π-1311\right)}{9801}$--for the eight atoms (Fig. 11), summing to 1, which span a 256-dimensional three-set (P, S, PPT) entanglement-probability boolean algebra for the two-qutrit Hiesmayr-Loffler states. PPT denotes positive partial transpose, while P and S provide the Li-Qiao necessary and}sufficient conditions for entanglement. The constraints ensuring entanglement are $s> \frac{16}{9} \approx 1.7777$ and $p> \frac{2^{27}}{3^{18} \cdot 7^{15} \cdot13} \approx 5.61324 \cdot 10^{-15}$. Here, $s$ is the square of the sum (Ky Fan norm) of the eight singular values of the $8 \times 8$ correlation matrix in the Bloch representation, and $p$, the square of the product of the singular values. In the two-ququart Hiesmayr-Loffler case, one constraint is $s>\frac{9}{4} \approx 2.25$, while $\frac{3^{24}}{2^{134}} \approx 1.2968528306 \cdot 10^{-29}$ is an upper bound on the appropriate $p$ value, with an entanglement probability $\approx 0.607698$. The $S$ constraints, in both cases, prove equivalent to the well-known CCNR/realignment criteria. Further, we detect and verify--using software of A. Mandilara--pseudo-one-copy undistillable (POCU) negative partial transposed two-qutrit states distributed over the surface of the separable states. Additionally, we study the best separable approximation problem within this two-qutrit setting, and obtain explicit decompositions of separable states into the sum of eleven product states. Numerous quantities of interest--including the eight atoms--were, first, estimated using a quasirandom procedure.

preprint2020arXiv

Archipelagos of Total Bound and Free Entanglement

First, we considerably simplify an initially quite complicated formula -- involving dilogarithms. It yields the total bound entanglement probability ($\approx 0.0865542$) for a qubit-ququart ($2 \times 4$) three-parameter model, recently analyzed for its separability properties by Li and Qiao. An "archipelago" of disjoint bound-entangled regions appears in the space of parameters, somewhat similarly to those recently found in our preprint, "Jagged Islands of Bound Entanglement and Witness-Parameterized Probabilities". There, two-qutrit and two-ququart Hiesmayr-L{ö}ffler "magic simplices" and generalized Horodecki states had been examined. However, contrastingly, in the present study, the entirety of bound entanglement--given by the formula obtained--is clearly captured in the archipelago found. Further, we "upgrade" the qubit-ququart model to a two-ququart one, for which we again find a bound-entangled archipelago, with its total probability simply being now $\frac{1}{729} \left(473-512 \log \left(\frac{27}{16}\right) \left(1+\log \left(\frac{27}{16}\right)\right)\right) \approx 0.0890496$. Then, "downgrading" the qubit-ququart model to a two-qubit one, we find an archipelago of total non-bound/free entanglement probability $\frac{1}{2}$.

preprint2020arXiv

Archipelagos of Total Bound and Free Entanglement. II

In the indicated preceding preprint (I), we reported the results of, in particular interest here, certain three-parameter qubit-ququart ($2 \times 4$) and two-ququart ($4 \times 4$) analyses. In them, we relied upon entanglement constraints given by Li and Qiao. However, further studies of ours conclusively show--using the well-known necessary and sufficient conditions for positive-semidefiniteness that all leading minors (of separable components, in this context) be nonnegative--that certain of the constraints given are flawed and need to be replaced (by weaker ones). Doing so, leads to a new set of results, somewhat qualitatively different and, in certain respects, simpler in nature. For example, bound-entanglement probabilities of $\frac{2}{3} \left(\sqrt{2}-1\right) \approx 0.276142$, $\frac{1}{4} \left(3-2 \log ^2(2)-\log (4)\right) \approx 0.1632$, $\frac{1}{2}-\frac{2}{3 π^2} \approx 0.432453$ and $\frac{1}{6}$, are reported for various implementations of constraints. We also adopt the Li-Qiao three-parameter framework to a two-parameter one, with interesting visual results.

preprint2020arXiv

Jagged Islands of Bound Entanglement and Witness-Parameterized Probabilities

We report several witness-parameterized families of bound-entangled probabilities. Two pertain to the $d=3$ (two-qutrit) and a third to the $d=4$ (two-ququart) subsets analyzed by Hiesmayr and L{ö}ffler of "magic" simplices of Bell states. The Hilbert-Schmidt probabilities of positive-partial-transpose (PPT) states--within which we search for bound-entangled states--are $\frac{8 π}{27 \sqrt{3}} \approx 0.537422$ ($d=3$) and $\frac{1}{2}+\frac{\log \left(2-\sqrt{3}\right)}{8 \sqrt{3}} \approx 0.404957$ ($d=4$). We obtain bound-entangled probabilities of $-\frac{4}{9}+\frac{4 π}{27 \sqrt{3}}+\frac{\log (3)}{6} \approx 0.00736862$ and $\frac{-204+7 \log (7)+168 \sqrt{3} \cos ^{-1}\left(\frac{11}{14}\right)}{1134} \approx 0.00325613$ ($d=3$) and $\frac{8 \log (2)}{27}-\frac{59}{288} \approx 0.00051583$ and $\frac{24 \text{csch}^{-1}\left(\frac{8}{\sqrt{17}}\right)}{17 \sqrt{17}}-\frac{91}{544} \approx 0.00218722$ ($d=4$). Thus, the total entanglement probability appears to equal $(1-\frac{8 π}{27 \sqrt{3}})+\frac{2}{81} \left(4 \sqrt{3} π-21\right) = \frac{13}{27} \approx 0.481481$.) The families, encompassing these results, are parameterized using generalized Choi and Jafarizadeh-Behzadi-Akbari witnesses. The same bound-entangled probability was achieved with both--the sets ("jagged islands") detected having void intersection. The entanglement (bound and "non-bound"/"free") probability for both was $\frac{1}{6} \approx 0.16667$, while their union and intersection gave $\frac{2}{9} \approx 0.22222$ and $\frac{1}{9} \approx 0.11111$. Further, we examine generalized Horodecki states, as well as estimating PPT-probabilities of approximately 0.39339 (very well-fitted by $\frac{7 π}{25 \sqrt{5}} \approx 0.39338962$) and 0.115732 ( for the original (8- [two-qutrit] and 15 [two-ququart]-dimensional) magic simplices themselves.

preprint2016arXiv

Bloch Radii Repulsion in Separable Two-Qubit Systems

Milz and Strunz recently reported substantial evidence to further support the previously conjectured separability probability of $\frac{8}{33}$ for two-qubit systems ($ρ$) endowed with Hilbert-Schmidt measure. Additionally, they found that along the radius ($r$) of the Bloch ball representing either of the two single-qubit subsystems, this value appeared constant (but jumping to unity at the locus of the pure states, $r=1$). Further, they also observed (personal communication) such separability probability $r$-invariance, when using, more broadly, random induced measure ($K=3,4,5,\ldots$), with $K=4$ corresponding to the (symmetric) Hilbert-Schmidt case. Among the findings here is that this invariance is maintained even after splitting the separability probabilities into those parts arising from the determinantal inequality $|ρ^{PT}| >|ρ|$ and those from $|ρ| > |ρ^{PT}| >0$, where the partial transpose is indicated. The nine-dimensional set of generic two-re[al]bit states endowed with Hilbert-Schmidt measure is also examined, with similar $r$-invariance conclusions. Contrastingly, two-qubit separability probabilities based on the Bures (minimal monotone) measure {\it diminish} with $r$. Moreover, we study the forms that the separability probabilities take as joint (bivariate) functions of the radii ($r_A, r_B$) of the Bloch balls of {\it both} single-qubit subsystems. Here, a form of Bloch radii {\it repulsion} for separable two-qubit systems emerges in {\it all} our several analyses. Separability probabilities tend to be smaller when the lengths of the two radii are closer. In Appendix A, we report certain companion analytic results for the much-investigated, more amenable (7-dimensional) $X$-states model.

preprint2016arXiv

Invariance of Bipartite Separability and PPT-Probabilities over Casimir Invariants of Reduced States

Milz and Strunz ({\it J. Phys. A}: {\bf{48}} [2015] 035306) recently studied the probabilities that two-qubit and qubit-qutrit states, randomly generated with respect to Hilbert-Schmidt (Euclidean/flat) measure, are separable. They concluded that in both cases, the separability probabilities (apparently exactly $\frac{8}{33}$ in the two-qubit scenario) hold {\it constant} over the Bloch radii ($r$) of the single-qubit subsystems, jumping to 1 at the pure state boundaries ($r=1$). Here, firstly, we present evidence that in the qubit-qutrit case, the separability probability is uniformly distributed, as well, over the {\it generalized} Bloch radius ($R$) of the qutrit subsystem. While the qubit (standard) Bloch vector is positioned in three-dimensional space, the qutrit generalized Bloch vector lives in eight-dimensional space. The radii variables $r$ and $R$ themselves are the lengths/norms (being square roots of {\it quadratic} Casimir invariants) of these ("coherence") vectors. Additionally, we find that not only are the qubit-qutrit separability probabilities invariant over the quadratic Casimir invariant of the qutrit subsystem, but apparently also over the {\it cubic} one--and similarly the case, more generally, with the use of random induced measure. We also investigate two-qutrit ($3 \times 3$) and qubit-{\it qudit} ($2 \times 4$) systems--with seemingly analogous {\it positive-partial-transpose}-probability invariances holding over what have been termed by Altafini, the {\it partial} Casimir invariants of these systems.

preprint2016arXiv

Two-Qubit Separability Probabilities as Joint Functions of the Bloch Radii of the Qubit Subsystems

We detect a certain pattern of behavior of separability probabilities $p(r_A,r_B)$ for two-qubit systems endowed with Hilbert-Schmidt, and more generally, random induced measures, where $r_A$ and $r_B$ are the Bloch radii ($0 \leq r_A,r_B \leq 1$) of the qubit reduced states ($A,B$). We observe a relative repulsion of radii effect, that is $p(r_A,r_A) < p(r_A,1-r_A)$, except for rather narrow "crossover" intervals $[\tilde{r}_A,\frac{1}{2}]$. Among the seven specific cases we study are, firstly, the "toy" seven-dimensional $X$-states model and, then, the fifteen-dimensional two-qubit states obtained by tracing over the pure states in $4 \times K$-dimensions, for $K=3, 4, 5$, with $K=4$ corresponding to Hilbert-Schmidt (flat/Euclidean) measure. We also examine the real (two-rebit) $K=4$, the $X$-states $K=5$, and Bures (minimal monotone)--for which no nontrivial crossover behavior is observed--instances. In the two $X$-states cases, we derive analytical results, for $K=3, 4$, we propose formulas that well-fit our numerical results, and for the other scenarios, rely presently upon large numerical analyses. The separability probability crossover regions found expand in length (lower $\tilde{r}_A$) as $K$ increases. This report continues our efforts (arXiv:1506.08739) to extend the recent work of Milz and Strunz (J. Phys. A: 48 [2015] 035306) from a univariate ($r_A$) framework---in which they found separability probabilities to hold constant with $r_A$---to a bivariate ($r_A,r_B$) one. We also analyze the two-qutrit and qubit-qutrit counterparts reported in arXiv:1512.07210 in this context, and study two-qubit separability probabilities of the form $p(r_A,\frac{1}{2})$.

preprint2015arXiv

Formulas for Rational-Valued Separability Probabilities of Random Induced Generalized Two-Qubit States

Previously, a formula, incorporating a $5F4$ hypergeometric function, for the Hilbert-Schmidt-averaged determinantal moments $\left\langle \left\vert ρ^{PT}\right\vert ^{n}\left\vert ρ\right\vert ^{k}\right\rangle /\left\langle \left\vert ρ\right\vert ^{k}\right\rangle$ of $4 \times 4$ density-matrices ($ρ$), and their partial transposes ($ρ^{PT}|$) was applied with $k=0$ to the generalized two-qubit separability-probability question. The formula can, further, be viewed we note here, as an averaging over "induced measures in the space of mixed quantum states". The associated induced-measure separability probabilities ($k =1, 2,\ldots$) are found--{\it via} a high-precision density approximation procedure--to assume interesting, relatively simple rational values in the two-re[al]bit ($α= \frac{1}{2}$), (standard) two-qubit ($α= 1$) and two-quater[nionic]bit ($α=2$) cases. We deduce rather simple companion (rebit, qubit, quaterbit,\ldots) formulas that successfully reproduce the rational values assumed for {\it general} $k$. These formulas are observed to share certain features, possibly allowing them to be incorporated into a single master formula.

preprint2015arXiv

Generalized Two-Qubit Whole and Half Hilbert-Schmidt Separability Probabilities

Compelling evidence-though yet no formal proof-has been adduced that the probability that a generic (standard) two-qubit state ($ρ$) is separable/disentangled is $\frac{8}{33}$ (arXiv:1301.6617, arXiv:1109.2560, arXiv:0704.3723). Proceeding in related analytical frameworks, using a further determinantal $4F3$-hypergeometric moment formula (Appendix A), we reach, {\it via} density-approximation procedures, the conclusion that one-half ($\frac{4}{33}$) of this probability arises when the determinantal inequality $|ρ^{PT}|>|ρ|$, where $PT$ denotes the partial transpose, is satisfied, and, the other half, when $|ρ|>|ρ^{PT}|$. These probabilities are taken with respect to the flat, Hilbert-Schmidt measure on the fifteen-dimensional convex set of $4 \times 4$ density matrices. We find fully parallel bisection/equipartition results for the previously adduced, as well, two-"re[al]bit" and two-"quater[nionic]bit"separability probabilities of $\frac{29}{64}$ and $\frac{26}{323}$, respectively. The new determinantal $4F3$-hypergeometric moment formula is, then, adjusted (Appendices B and C) to the boundary case of minimally degenerate states ($|ρ|=0$), and its consistency manifested-also using density-approximation-with a theorem of Szarek, Bengtsson and {Ż}yczkowski (arXiv:quant-ph/0509008). This theorem states that the Hilbert-Schmidt separability probabilities of generic minimally degenerate two-qubit states are (again) one-half those of the corresponding generic nondegenerate states.

preprint2015arXiv

Hypergeometric/Difference-Equation-Based Separability Probability Formulas and Their Asymptotics for Generalized Two-Qubit States Endowed with Random Induced Measure

We find equivalent hypergeometric- and difference-equation-based formulas, $Q(k,α)= G_1^k(α) G_2^k(α)$, for $k = -1, 0, 1,\ldots,9$, for that (rational-valued) portion of the total separability probability for generalized two-qubit states endowed with random induced measure, for which the determinantal inequality $|ρ^{PT}| >|ρ|$ holds. Here $ρ$ denotes a $4 \times 4$ density matrix and $ρ^{PT}$, its partial transpose, while $α$ is a Dyson-index-like parameter with $α= 1$ for the standard (15-dimensional) convex set of two-qubit states. The dimension of the space in which these density matrices is embedded is $4 \times (4 +k)$. For the symmetric case of $k=0$, we obtain the previously reported Hilbert-Schmidt formulas, with (the two-re[al]bit case) $Q(0,\frac{1}{2}) = \frac{29}{128}$, (the standard two-qubit case) $Q(0,1)=\frac{4}{33}$, and (the two-quater[nionic]bit case) $Q(0,2)= \frac{13}{323}$. The factors $G_2^k(α)$ can be written as the sum of weighted hypergeometric functions $_{p}F_{p-1}$, $p \geq 7$, all with argument $\frac{27}{64} =(\frac{3}{4})^3$. We find formulas for the upper and lower parameter sets of these functions and, then, equivalently express $G_2^k(α)$ in terms of first-order difference equations. The factors $G_1^k(α)$ are equal to $(\frac{27}{64})^{α-1}$ times ratios of products of six Pochhammer symbols involving the indicated parameters. Some remarkable $α-$ and $k$-specific invariant asymptotic properties (again, involving $\frac{27}{64}$ and related quantities) of separability probability formulas emerge.

preprint2015arXiv

Separability Probability Formulas and Their Proofs for Generalized Two-Qubit X-Matrices Endowed with Hilbert-Schmidt and Induced Measures

Two-qubit X-matrices have been the subject of considerable recent attention, as they lend themselves more readily to analytical investigations than two-qubit density matrices of arbitrary nature. Here, we maximally exploit this relative ease of analysis to formally derive an exhaustive collection of results pertaining to the separability probabilities of generalized two-qubit X-matrices endowed with Hilbert-Schmidt and, more broadly, induced measures. Further, the analytical results obtained exhibit interesting parallels to corresponding earlier (but, contrastingly, not yet fully rigorous) results for general 2-qubit states--deduced on the basis of determinantal moment formulas. Geometric interpretations can be given to arbitrary positive values of the random-matrix Dyson-index-like parameter $α$ employed.

preprint2014arXiv

A Balanced $_{10}F_{9}$ Hypergeometric Hybrid Hilbert-Schmidt/Bures Two-Qubit Function and Related Constructions

We seek to develop a Bures (minimal monotone/statistical distinguishability) metric-based series of formulas for the moments of probability distributions over the determinants $|ρ|$ and $|ρ^{PT}|$ of $4 \times 4$ density matrices, $ρ$, for generalized (rebit, quater[nionic]bit,\ldots) two-qubit systems, analogous to a series that has been obtained for the Hilbert-Schmidt (HS) metric. In particular, we desire--using moment-inversion procedures--to be able to closely test the previously-developed conjecture (J. Geom. Phys., 53, 74 [2005]) that the Bures separability probability over the (standard, fifteen-dimensional convex set of) two-qubit states is $\frac{1680 \left(\sqrt{2}-1\right)}{π^8} \approx 0.0733389$--while, in the HS context, strong evidence has been adduced, along the indicated analytical lines, that the counterpart of this value is $\frac{8}{33}$ (J. Phys. A}, 45, 095305 [2012]). Working within the "utility function" framework of Dunkl employed in that latter study, we obtain an interesting $_{10}F_{9}$ balanced hypergeometric function based on a "hybridization" of known Bures and HS terms. This exercise appears to provide an upper bound on the Bures two-qubit separability probability of 0.0798218. We also examine the yet unresolved HS qubit-qutrit scenario. Mathematica calculations indicate that if the same form of hypergeometric paradigm as has been established for the generalized two-qubit HS moments is followed in either the HS qubit-qutrit or Bures two-qubit cases, then a balanced hypergeometric function $_{p}F_{p-1}$ with $p>9$ would be required.

preprint2014arXiv

Symmetric Halves of the 8/33-Probability that the Joint State of Two Quantum Bits is Disentangled

Compelling evidence-though yet no formal proof--has been adduced that the probability that a generic two-qubit state ($ρ$) is separable is $\frac{8}{33}$ (arXiv:1301.6617, arXiv:1109.2560, arXiv:0704.3723). Proceeding in related analytical frameworks, using a further determinantal moment formula of C. Dunkl (Appendix), we reach the conclusion that one-half of this probability arises when the determinantal inequality $|ρ^{PT}|>|ρ|$, where $PT$ denotes the partial transpose, is satisfied, and, the other half, when $|ρ|>|ρ^{PT}|$. These probabilities are taken with respect to the flat, Hilbert-Schmidt measure on the fifteen-dimensional convex set of $4 \times 4$ density matrices. We find fully parallel bisection/equipartition results for the previously adduced, as well, two-"re[al]bit" and two-"quater[nionic]bit"separability probabilities of $\frac{29}{64}$ and $\frac{26}{323}$, respectively. The computational results reported lend strong support to those obtained earlier--including the "concise formula" that yields them--most conspicuously amongst those findings being the $\frac{29}{64}, \frac{8}{33}$ and $\frac{26}{323}$ probabilities noted.

preprint2013arXiv

A Concise Formula for Generalized Two-Qubit Hilbert-Schmidt Separability Probabilities

We report major advances in the research program initiated in "Moment-Based Evidence for Simple Rational-Valued Hilbert-Schmidt Generic 2 x 2 Separability Probabilities" (J. Phys. A, 45, 095305 [2012]). A highly succinct separability probability function P(alpha) is put forth, yielding for generic (9-dimensional) two-rebit systems, P(1/2) = 29/64, (15-dimensional) two-qubit systems, P(1) = 8/33 and (27-dimensional) two-quater(nionic)bit systems, P(2)=26/323. This particular form of P(alpha) was obtained by Qing-Hu Hou and colleagues by applying Zeilberger's algorithm ("creative telescoping") to a fully equivalent--but considerably more complicated--expression containing six 7F6 hypergeometric functions (all with argument 27/64 = (3/4)^3). That hypergeometric form itself had been obtained using systematic, high-accuracy probability-distribution-reconstruction computations. These employed 7,501 determinantal moments of partially transposed four-by-four density matrices, parameterized by alpha = 1/2, 1, 3/2,...,32. From these computations, exact rational-valued separability probabilities were discernible. The (integral/half-integral) sequences of 32 rational values, then, served as input to the Mathematica FindSequenceFunction command, from which the initially obtained hypergeometric form of P(alpha) emerged.

preprint2013arXiv

Two-Qubit Rational-Valued Entanglement-Boundary Probability Densities and a Fisher Information Equality Conjecture

We consider a pair of one-parameter (alpha) families of generalized two-qubit determinantal Hilbert-Schmidt probability distributions, p_{alpha}(|rho^{PT}|) and q_{alpha}(|rho|), where rho is a 4 x 4 density matrix, rho^{PT}, its partial transpose, with |rho^{PT}| \in [-1/16,1/256] and |rho| \in [0, 1/256]. The Dyson-index-like (random matrix) parameter alpha is 1/2 for the 9-dimensional generic two-rebit systems, 1 for the 15-dimensional generic two-qubit systems,... Numerical (moment-based probability-distribution-reconstruction) analyses suggest the conjecture that the Fisher information--a measure over alpha--is identical for the two distinct families. Further, we study the values of p_{alpha}(0), the probability densities at the separability-entanglement boundary, with evidence strongly indicating that p_2(0) =7425/34 and p_3(0)= 7696/69. Despite extensive results of such a nature, we have not yet succeeded--in contrast with the corresponding rational-valued separability probabilities (arXiv.org:1301.6617)--in generating an underlying, explanatory ("concise") formula for p_α(0), even though the corresponding denominators of both sets of rational-valued results are closely related, having almost identical (small) prime factors. The first derivatives p_{alpha}^{'}(0) are positive for alpha = 1/2 and 1, but negative for alpha > 1.

preprint2012arXiv

A Further (Itakura-Saito/beta=0) Bi-stochaticization and Associated Clustering/Regionalization of the 3,107-County 1995-2000 U. S. Migration Network

We extend to the beta-divergence (Itakura-Saito) case beta =0, the comparative bi-stochaticization analyses-previously conducted (arXiv:1208.3428) for the (Kullback-Leibler) beta=1 and (squared-Euclidean) beta = 2 cases -of the 3,107 - county 1995-2000 U. S. migration network. A heuristic, "greedy" algorithm is devised. While the largest 25,329 entries of the 735,531 non-zero entries of the bi-stochasticized table - in the beta=1 case - are required to complete the widely-applied two-stage (double-standardization and strong-component hierarchical clustering) procedure, 105,363 of the 735,531 are needed (reflective of greater uniformity of entries) in the beta=0 instance. The North Carolina counties of Mecklenburg (Charlotte) and Wake (Raleigh) are considerably relatively more cosmopolitan in the beta=0 study. The Colorado county of El Paso (Colorado Springs) replaces the Florida Atlantic county of Brevard (the "Space Coast") as the most cosmopolitan, with Brevard becoming the second-most. Honolulu County splinters away from the other four (still-grouped) Hawaiian counties, becoming the fifth most cosmopolitan county nation-wide. The five counties of Rhode Island remain intact as a regional entity, but the eight counties of Connecticut fragment, leaving only five counties clustered.

preprint2012arXiv

A Hypergeometric Formula Yielding Hilbert-Schmidt Generic 2 x 2 Generalized Separability Probabilities

We significantly advance the research program initiated in "Moment-Based Evidence for Simple Rational-Valued Hilbert-Schmidt Generic 2 x 2 Separability Probabilities" (J. Phys. A, 45, 095305 [2012]). A function P(alpha), incorporating a family of six hypergeometric functions, all with argument 27/64 =(3/4)^3, is obtained here. It reproduces a series, alpha = 1/2, 1, 3/2, 2,...,32$, of sixty-four Hilbert-Schmidt generic 2 x 2 generalized separability probabilities, advanced on the basis of systematic high-accuracy probability-distribution-reconstruction computations, employing 7,501 determinantal moments of partially transposed 4 x 4 density matrices. For generic (9-dimensional) two-rebit systems, P(1/2) = 29/64, (15-dimensional) two-qubit, P(1) = 8/33, and (27-dimensional) two-quat(ernionic)bit systems, P(2)= 26/323. The function P(alpha) is generated--applying the Mathematica command FindSequenceFunction--using solely either) the thirty-two integral or half-integral probabilities, yet successfully predicts the unused thirty-two.

preprint2012arXiv

Bures and Hilbert-Schmidt 2 x 2 Determinantal Moments

We seek to gain insight into the nature of the determinantal moments of generic (9-dimensional) two-rebit and (15-dimensional) two-qubit systems (rho). Such information-as it has proved to be in the Hilbert-Schmidt counterpart--should be useful, employing probability-distribution reconstruction (inverse) procedures, in obtaining improved, or possibly even exact Bures 2 x 2 separability probabilities for such systems. The (regularizing) strategy we first adopt is to plot the ratio of numerically-generated (Ginibre ensemble) estimates of the Bures moments to the corresponding (apparently) exactly-known Hilbert-Schmidt moments (J. Phys. A, 45, 095305 [2012]). Then, through a combination of symbolic and numerical computations, we obtain strong evidence as to the exact values (and underlying patterns) of certain Bures moments. In particular, the first moment (average) of |rho^{PT}| (where |rho^{PT}| in [-1/16, 1/256]), PT denoting partial transpose, for the two-qubit systems is, remarkably, -1 / 256. The analogous value for the two-rebit systems is -2663/(2^(13) x 3 x 5 x 7) approx -0.00309594. While, an important function in the Hilbert-Schmidt case is the ratio of 3n-degree polynomials, it appears to be the ratio of 5n-degree polynomials in the Bures case.

preprint2012arXiv

Comparative Bi-stochastizations and Associated Clusterings/Regionalizations of the 1995-2000 U. S. Intercounty Migration Network

Wang, Li and Konig have recently compared the cluster-theoretic properties of bi-stochasticized symmetric data similarity (e. g. kernel) matrices, produced by minimizing two different forms of Bregman divergences. We extend their investigation to non-symmetric matrices, specifically studying the 1995-2000 U. S. 3,107 x 3,107 intercounty migration matrix. A particular bi-stochastized form of it had been obtained (arXiv:1207.0437), using the well-established Sinkhorn-Knopp (SK) (biproportional) algorithm--which minimizes the Kullback-Leibler form of the divergence. This matrix has but a single entry equal to (the maximal possible value of) 1. Highly contrastingly, the bi-stochastic matrix obtained here, implementing the Wang-Li-Konig-algorithm for the minimum of the alternative, squared-norm form of the divergence, has 2,707 such unit entries. The corresponding 3,107-vertex, 2,707-link directed graph has 2,352 strong components. These consist of 1,659 single/isolated counties, 654 doublets (thirty-one interstate in nature), 22 triplets (one being interstate), 13 quartets (one being interstate), three quintets and one septet. Not manifest in these graph-theoretic results, however, are the five-county states of Hawaii and Rhode Island and the eight-county state of Connecticut. These--among other regional configurations--appealingly emerged as well-defined entities in the SK-based strong-component hierarchical clustering.

preprint2012arXiv

Moment-Based Evidence for Simple Rational-Valued Hilbert-Schmidt Generic 2 x 2 Separability Probabilities

Employing Hilbert-Schmidt measure, we explicitly compute and analyze a number of determinantal product (bivariate) moments |rho|^k |rho^{PT}|^n, k,n=0,1,2,3,..., PT denoting partial transpose, for both generic (9-dimensional) two-rebit (alpha = 1/2) and generic (15-dimensional) two-qubit (alpha=1) density matrices rho. The results are, then, incorporated by Dunkl into a general formula (Appendix D6), parameterized by k, n and alpha, with the case alpha=2, presumptively corresponding to generic (27-dimensional) quaternionic systems. Holding the Dyson-index-like parameter alpha fixed, the induced univariate moments (|rho| |rho^{PT}|)^n and |rho^{PT}|^n are inputted into a Legendre-polynomial-based (least-squares) probability-distribution reconstruction algorithm of Provost (Mathematica J., 9, 727 (2005)), yielding alpha-specific separability probability estimates. Since, as the number of inputted moments grows, estimates based on |rho| |rho^{PT}| strongly decrease, while ones employing |rho^{PT}| strongly increase (and converge faster), the gaps between upper and lower estimates diminish, yielding sharper and sharper bounds. Remarkably, for alpha = 2, with the use of 2,325 moments, a separability-probability lower-bound 0.999999987 as large as 26/323 = 0.0804954 is found. For alpha=1, based on 2,415 moments, a lower bound results that is 0.999997066 times as large as 8/33 = 0.242424, a (simpler still) fractional value that had previously been conjectured (J. Phys. A, 40, 14279 (2007)). Furthermore, for alpha = 1/2, employing 3,310 moments, the lower bound is 0.999955 times as large as 29/64 = 0.453125, a rational value previously considered (J. Phys. A, 43, 195302 (2010)).

preprint2012arXiv

Toward Exact 2 x 2 Hilbert-Schmidt Determinantal Probability Distributions via Mellin Transforms and Other Approaches

We attempt to construct the exact univariate probability distributions for 2 x 2 quantum systems that yield the (balanced) univariate Hilbert-Schmidt determinantal moments <(|rho| |rho^{PT}|)^n>, obtained by Slater and Dunkl (J. Phys. A, 45, 095305 [2012]). To begin, we follow--to the extent possible--the Mellin transform-based approach of Penson and Zyczkowski in their study of Fuss-Catalan and Raney distributions (Phys. Rev. E, 83, 061118 [2011]). Further, we approximate the y-intercepts (separability/entanglement boundaries)--at which |rho^{PT}|=0-- of the associated probability distributions based on the (balanced) moments, as well as the previously reported unbalanced determinantal moments <|rho^{PT}|^n>, as a function of the seventy-two values of the Dyson-index-like parameter alpha = 1/2 (rebits), 1 (qubits),...,2 (quaterbits),...35.

preprint2012arXiv

Two-Qubit Separability Probabilities: A Concise Formula

We report a concise answer--in the case of 2 x 2 systems--to the fundamental quantum-information-theoretic question as to "the volume of separable states" posed by Zyczkowski, Horodecki, Sanpera and Lewenstein (Phys. Rev. A, 58, 883 [1998]). We proceed by applying the Mathematica command FindSequenceFunction to a series of conjectured Hilbert-Schmidt generic 2 x 2 (rational-valued) separability probabilities p(a), a = 1, 2,...,32, with a = 1 indexing standard two-qubit systems, and a = 2, two-quater(nionic)bit systems. These 32 inputted values of p(a)--as well as 32 companion non-inputted values for the half-integers, a = 1/2 (two-re[al]bit) systems), 3/2,..., 63/2, are advanced on the basis of high-precision probability-distribution-reconstruction computations, employing 7,501 determinantal moments of partially transposed 4 x 4 density matrices. The function P(a) given by application of the command fully reproduces both of these 32-length sequences, and an equivalent outcome is obtained if the half-integral series is the one inputted. The lengthy expression (containing six hypergeometric functions) obtained for P(a) is, then, impressively condensed (by Qing-Hu Hou and colleagues), using Zeilberger's algorithm. For generic (9-dimensional) two-rebit systems, P(1/2) = 29/64, (15-dimensional) two-qubit systems, P(1) = 8/33, (27-dimensional) two-quaterbit systems, P(2) = 26/323, while for generic classical (3-dimensional) systems, P(0)=1.

preprint2011arXiv

Hilbert-Schmidt Orthogonality of det(rho) and det(rho^{PT}) over the Two-Rebit Systems rho and Further Determinantal Moment Analyses

A complete description of the multitudinous ways in which quantum particles can be entangled requires the use of high-dimensional abstract mathematical spaces. We report here a particularly interesting feature of the nine-dimensional convex set-endowed with Hilbert-Schmidt (Euclidean/flat) measure-composed of two-rebit (4 x 4) density matrices (rho). To each rho is assigned the product of its (nonnegative) determinant |rho| and the determinant of its partial transpose |rho^{PT}| -negative values of which, by the results of Peres and Horodecki, signify the entanglement of rho. Integrating this product, |rho| |rho^{PT}| =|rho rho^{PT}|, over the nine-dimensional space, using the indicated (HS) measure, we obtain the result zero. The two determinants, thus, form a pair of multivariate orthogonal polynomials with respect to HS measure. It is hypothesized that the analogous two determinants for the standard fifteen-dimensional convex set of two-qubit density matrices are similarly orthogonal. However, orthogonality does not hold, we find, if the symmetry of the nine-dimensional two-rebit scenario is broken slightly, nor with the use of non-flat measures, such as the prominent Bures (minimal monotone) measure-nor in the full HS extension to the 15-dimensional convex set of two-qubit density matrices. We discuss relations-involving the HS moments of |ρ^{PT}|-to the long-standing problem of determining the probability that a generic pair of rebits/qubits is separable.

preprint2011arXiv

Joint Hilbert-Schmidt Determinantal Moments of Product Form for the Two-Rebit and Two-Qubit and Higher-Dimensional Quantum Systems

We report formulas for the joint moments of the determinantal products (det{rho})^k (det{rho^PT})^K (k=0, 1, 2,...,N; K = 0, 1, 2, 3, 4) of Hilbert-Schmidt (HS) probability distributions over the two-rebit and (K = 0, 1) two-qubit density matrices rho. Here PT denotes the partial transposition operation of quantum-information-theoretic importance. Each formula is the product of the expression for the moments of (det{rho})^k, k=0,1,2,...,N--already known from work of Cappellini, Sommers and Zyczkowski (Phys. Rev. A 74, 062322 (1996))--and an adjustment factor. The factor is a biproper rational function, with its numerators and denominators both being 3 K-degree polynomials in k. We infer the structure that the denominators follow for arbitrary K in both the two-rebit and two-qubit cases, and the six leading-order coefficients of k of the numerators in the two-rebit scenario.. We also commence an analogous investigation of generic rebit-retrit and qubit-qutrit systems. This research was motivated, in part, by the objective of using the computed moments to well reconstruct the associated HS probabilities (for which we present plots) over the two determinants, and to ascertain--at least to high accuracy--the associated (separability) probabilities of "philosophical, practical and physical" interest that (det{rho^PT})>0.

preprint2011arXiv

Radial and Azimuthal Profiles of Two-Qubit/Rebit Hilbert-Schmidt Separability Probabilities and Related 3-D Visualization Analyses

Firstly, we reduce the long-standing problem of ascertaining the Hilbert-Schmidt probability that a generic pair of qubits is separable to that of determining the specific nature of a one-dimensional (separability) function of the radial coordinate (r) of the unit ball in 15-dimensional Euclidean space, and similarly for a generic pair of rebits, using the 9-dimensional unit ball. Separability probabilities, could, then, be directly obtained by integrating the products of these functions (which we numerically estimate and plot) with jacobian factors of r^m over r in [0,1], with m=17 for the two-rebit case, and m=29 in the two-qubit instance. Secondly, we repeat the analyses, but for the replacement of r as the free variable, by the azimuthal angle phi \in [0,2 pi]--with the associated jacobian factors now being, trivially, unity. So, the separability probabilities, then, become simply the areas under the curves generated. For our analyses, we employ an interesting Cholesky-decomposition parameterization of the 4 x 4 density matrices. Thirdly, in an exploratory investigation, we examine the Hilbert-Schmidt separability probability question using the three-dimensional visualization (cube/tetrahedron/octahedron) of two qubits associated with Avron, Bisker and Kenneth, and Leinaas, Myrheim and Ovrum, among others.

preprint2010arXiv

Advances in delimiting the Hilbert-Schmidt separability probability of real two-qubit systems

We seek to derive the probability--expressed in terms of the Hilbert-Schmidt (Euclidean or flat) metric--that a generic (nine-dimensional) real two-qubit system is separable, by implementing the well-known Peres-Horodecki test on the partial transposes (PT's) of the associated 4 x 4 density matrices). But the full implementation of the test--requiring that the determinant of the PT be nonnegative for separability to hold--appears to be, at least presently, computationally intractable. So, we have previously implemented--using the auxiliary concept of a diagonal-entry-parameterized separability function (DESF)--the weaker implied test of nonnegativity of the six 2 x 2 principal minors of the PT. This yielded an exact upper bound on the separability probability of 1024/{135 pi^2} =0.76854$. Here, we piece together (reflection-symmetric) results obtained by requiring that each of the four 3 x 3 principal minors of the PT, in turn, be nonnegative, giving an improved/reduced upper bound of 22/35 = 0.628571. Then, we conclude that a still further improved upper bound of 1129/2100 = 0.537619 can be found by similarly piecing together the (reflection-symmetric) results of enforcing the simultaneous nonnegativity of certain pairs of the four 3 x 3 principal minors. In deriving our improved upper bounds, we rely repeatedly upon the use of certain integrals over cubes that arise. Finally, we apply an independence assumption to a pair of DESF's that comes close to reproducing our numerical estimate of the true separability function.

preprint2010arXiv

Eigenvalues, Separability and Absolute Separability of Two-Qubit States

Substantial progress has recently been reported in the determination of the Hilbert-Schmidt (HS) separability probabilities for two-qubit and qubit-qutrit (real, complex and quaternionic) systems. An important theoretical concept employed has been that of a separability function. It appears that if one could analogously obtain separability functions parameterized by the eigenvalues of the density matrices in question--rather than the diagonal entries, as originally used--comparable progress could be achieved in obtaining separability probabilities based on the broad, interesting class of monotone metrics (the Bures, being its most prominent [minimal] member). Though large-scale numerical estimations of such eigenvalue-parameterized functions have been undertaken, it seems desirable also to study them in lower-dimensional specialized scenarios in which they can be exactly obtained. In this regard, we employ an Euler-angle parameterization of SO(4) derived by S. Cacciatori (reported in an Appendix)--in the manner of the SU(4)-density matrix parameterization of Tilma, Byrd and Sudarshan. We are, thus, able to find simple exact separability (inverse-sine-like) functions for two real two-qubit (rebit) systems, both having three free eigenvalues and one free Euler angle. We also employ the important Verstraete-Audenaert-de Moor bound to obtain exact HS probabilities that a generic two-qubit state is absolutely separable (that is, can not be entangled by unitary transformations). In this regard, we make copious use of trigonometric identities involving the tetrahedral dihedral angle arccos(1/3).

preprint2010arXiv

Moments of the Hilbert-Schmidt probability distributions over determinants of real two-qubit density matrices and of their partial transposes

The nonnegativity of the determinant of the partial transpose of a two-qubit (4 x 4) density matrix is both a necessary and sufficient condition for its separability. While the determinant is restricted to the interval [0,1/256], the determinant of the partial transpose can range over [-1/16,1/256], with negative values corresponding to entangled states. We report here the exact values of the first nine moments of the probability distribution of the partial transpose over this interval, with respect to the Hilbert-Schmidt (metric volume element) measure on the nine-dimensional convex set of real two-qubit density matrices. Rational functions C_{2 j}(m), yielding the coefficients of the 2j-th power of even polynomials occurring at intermediate steps in our derivation of the m-th moment, emerge. These functions possess poles at finite series of consecutive half-integers (m=-3/2,-1/2,...,(2j-1)/2), and certain (trivial) roots at finite series of consecutive natural numbers (m=0, 1,...). Additionally, the (nontrivial) dominant roots of C_{2 j}(m) approach the same half-integer values (m = (2 j-1)/2, (2 j-3)/2,...), as j increases. The first two moments (mean and variance) found--when employed in the one-sided Chebyshev inequality--give an upper bound of 30397/34749 = 0.874759 on the separability probability of real two-qubit density matrices. We are able to report general formulas for the m-th moment of the Hilbert-Schmidt probability distribution of the density matrix determinant over [0,1/256], in the real, complex and quaternionic two-qubit cases.

preprint2010arXiv

Multiscale Network Reduction Methodologies: Bistochastic and Disparity Filtering of Human Migration Flows between 3,000+ U. S. Counties

To control for multiscale effects in networks, one can transform the matrix of (in general) weighted, directed internodal flows to bistochastic (doubly-stochastic) form, using the iterative proportional fitting (Sinkhorn-Knopp) procedure, which alternatively scales row and column sums to all equal 1. The dominant entries in the bistochasticized table can then be employed for network reduction, using strong component hierarchical clustering. We illustrate various facets of this well-established, widely-applied two-stage algorithm with the 3, 107 x 3, 107 (asymmetric) 1995-2000 intercounty migration flow table for the United States. We compare the results obtained with ones using the disparity filter, for "extracting the "multiscale backbone of complex weighted networks", recently put forth by Serrano, Boguna and Vespignani (SBV) (Proc. Natl. Acad. Sci. 106 [2009], 6483), upon which we have briefly commented (Proc. Natl. Acad. Sci. 106 [2009], E66). The performance of the bistochastic filter appears to be superior-at least in this specific case-in two respects: (1) it requires far fewer links to complete a stongly-connected network backbone; and (2) it "belittles" small flows and nodes less-a principal desideratum of SBV-in the sense that the correlations of the nonzero raw flows are considerably weaker with the corresponding bistochastized links than with the significance levels yielded by the disparity filter. Additional comparative studies--as called for by SBV-of these two filtering procedures, in particular as regards their topological properties, should be of considerable interest. Relatedly, in its many geographic applications, the two-stage procedure has--with rare exceptions-clustered contiguous areas, often reconstructing traditional regions (islands, for example), even though no contiguity constraints, at all, are imposed beforehand.

preprint2009arXiv

A two-stage algorithm for extracting the multiscale backbone of complex weighted networks

The central problem of concern to Serrano, Boguna and Vespignani ("Extracting the multiscale backbone of complex weighted networks", Proc Natl Acad Sci 106:6483-6488 [2009]) can be effectively and elegantly addressed using a well-established two-stage algorithm that has been applied to internal migration flows for numerous nations and several other forms of "transaction flow data".

preprint2009arXiv

Matrix plots of reordered bistochastized transaction flow tables: A United States intercounty migration example

We present a number of variously rearranged matrix plots of the $3, 107 \times 3, 107$ 1995-2000 (asymmetric) intercounty migration table for the United States, principally in its bistochasticized form (all 3,107 row and column sums iteratively proportionally fitted to equal 1). In one set of plots, the counties are seriated on the bases of the subdominant (left and right) eigenvectors of the bistochastic matrix. In another set, we use the ordering of counties in the dendrogram generated by the associated strong component hierarchical clustering. Interesting, diverse features of U. S. intercounty migration emerge--such as a contrast in centralized, hub-like (cosmopolitan/provincial) properties between cosmopolitan "Sunbelt" and provincial "Black Belt" counties. The methodologies employed should also be insightful for the many other diverse forms of interesting transaction flow-type data--interjournal citations being an obvious, much-studied example, where one might expect that the journals Science, Nature and PNAS would display "cosmopolitan" characteristics.

preprint2009arXiv

Ratios of maximal concurrence-parameterized separability functions, and generalized Peres-Horodecki conditions

The probability that a generic real, complex or quaternionic two-qubit state is separable can be considered to be the sum of three contributions. One is from those states that are absolutely separable, that is those (which can not be entangled by unitary transformations) for which the maximal concurrence over spectral orbits (C_{max}) is zero. The other two contributions are from the states for which C_{max} in [0,1/2], and for which C_{max} in [1/2,1]. We have previously (arXiv:0805.0267) found exact formulas for the absolutely separable contributions in terms of the Hilbert-Schmidt metric over the quantum states, and here advance hypotheses as to the exact contributions for C_{max} in [1/2,1]. A crucial element in understanding the other two contributions is the nature of the ratio (R) of the C_{max}-parameterized separability function for the complex states to the square of the comparable function for the real states--both such functions having clearly displayed jump discontinuities at C_{max}=1/2. For C_{max} in [0,1/2], R appears to be of the form 1+ k C_{max}, except near C_{max}=1/2, while for C_{max} in [1/2,1], it seems to be a constant (thus according to the Dyson-index pattern of random matrix theory). Related phenomena also occur for the minimally-degenerate two-qubit states and the qubit-qutrit states. Our results have immediate application to the computation of separability probabilities in terms of other metrics, such as the Bures (minimal monotone) metric. The paper begins with continuous embeddings of the separability probability question in terms of four metrics of interest, using "generalized Peres-Horodecki conditions".

preprint2008arXiv

Discernment of Hubs and Clusters in Socioeconomic Networks

Interest in the analysis of networks has grown rapidly in the new millennium. Consequently, we promote renewed attention to a certain methodological approach introduced in 1974. Over the succeeding decade, this two-stage--double-standardization and hierarchical clustering (single-linkage-like)--procedure was applied to a wide variety of weighted, directed networks of a socioeconomic nature, frequently revealing the presence of ``hubs''. These were, typically--in the numerous instances studied of migration flows between geographic subdivisions within nations--``cosmopolitan/non-provincial'' areas, a prototypical example being the French capital, Paris. Such locations emit and absorb people broadly across their respective nations. Additionally, the two-stage procedure--which ``might very well be the most successful application of cluster analysis'' (R. C. Dubes, 1985)--detected many (physically or socially) isolated, functional groups (regions) of areas, such as the southern islands, Shikoku and Kyushu, of Japan, the Italian islands of Sardinia and Sicily, and the New England region of the United States. Further, we discuss a (complementary) approach developed in 1976, in which the max-flow/min-cut theorem was applied to raw/non-standardized (interindustry, as well as migration) flows.

preprint2008arXiv

Established Clustering Procedures for Network Analysis

In light of the burgeoning interest in network analysis in the new millenium, we bring to the attention of contemporary network theorists, a two-stage double-standarization and hierarchical clustering (single-linkage-like) procedure devised in 1974. In its many applications over the next decade--primarily to the migration flows between geographic subdivisions within nations--the presence was often revealed of ``hubs''. These are, typically, ``cosmopolitan/non-provincial'' areas--such as the French capital, Paris--which send and receive people relatively broadly across their respective nations. Additionally, this two-stage procedure--which ``might very well be the most successful application of cluster analysis'' (R. C. Dubes)--has detected many (physically or socially) isolated groups (regions) of areas, such as those forming the southern islands, Shikoku and Kyushu, of Japan, the Italian islands of Sardinia and Sicily, and the New England region of the United States. Further, we discuss a (complementary) approach developed in 1976, involving the application of the max-flow/min-cut theorem to raw/non-standardized flows.

preprint2008arXiv

Hubs and Clusters in the Evolving U. S. Internal Migration Network

Most nations of the world periodically publish N x N origin-destination tables, recording the number of people who lived in geographic subdivision i at time t and j at t+1. We have developed and widely applied to such national tables and other analogous (weighted, directed) socioeconomic networks, a two-stage--double-standardization and (strong component) hierarchical clustering--procedure. Previous applications of this methodology and related analytical issues are discussed. Its use is illustrated in a large-scale study, employing recorded United States internal migration flows between the 3,000+ county-level units of the nation for the periods 1965-1970 and 1995-2000. Prominent, important features--such as ''cosmopolitan hubs'' and ``functional regions''--are extracted from master dendrograms. The extent to which such characteristics have varied over the intervening thirty years is evaluated.

preprint2006arXiv

Two-Qubit Separability Probabilities and Beta Functions

Due to recent important work of Zyczkowski and Sommers (quant-ph/0302197 and quant-ph/0304041), exact formulas are available (both in terms of the Hilbert-Schmidt and Bures metrics) for the (n^2-1)-dimensional and (n(n-1)/2-1)-dimensional volumes of the complex and real n x n density matrices. However, no comparable formulas are available for the volumes (and, hence, probabilities) of various separable subsets of them. We seek to clarify this situation for the Hilbert-Schmidt metric for the simplest possible case of n=4, that is, the two-qubit systems. Making use of the density matrix (rho) parameterization of Bloore (J. Phys. A 9, 2059 [1976]), we are able to reduce each of the real and complex volume problems to the calculation of a one-dimensional integral, the single relevant variable being a certain ratio of diagonal entries, nu = (rho_{11} rho_{44})/{rho_{22} rho_{33})$. The associated integrand in each case is the product of a known (highly oscillatory near nu=1) jacobian and a certain unknown univariate function, which our extensive numerical (quasi-Monte Carlo) computations indicate is very closely proportional to an (incomplete) beta function B_{nu}(a,b), with a=1/2, b=sqrt{3}in the real case, and a=2 sqrt{6}/5, b =3/sqrt{2} in the complex case. Assuming the full applicability of these specific incomplete beta functions, we undertake separable volume calculations.

preprint2005arXiv

Hilbert-Schmidt Separability Probabilities and Noninformativity of Priors

The Horodecki family employed the Jaynes maximum-entropy principle, fitting the mean (b_{1}) of the Bell-CHSH observable (B). This model was extended by Rajagopal by incorporating the dispersion (σ_{1}^2) of the observable, and by Canosa and Rossignoli, by generalizing the observable (B_α). We further extend the Horodecki one-parameter model in both these manners, obtaining a three-parameter (b_{1},σ_{1}^2,α) two-qubit model, for which we find a highly interesting/intricate continuum (-\infty < α< \infty) of Hilbert-Schmidt (HS) separability probabilities -- in which, the golden ratio is featured. Our model can be contrasted with the three-parameter (b_{q}, σ_{q}^2,q) one of Abe and Rajagopal, which employs a q(Tsallis)-parameter rather than $α$, and has simply q-invariant HS separability probabilities of 1/2. Our results emerge in a study initially focused on embedding certain information metrics over the two-level quantum systems into a q-framework. We find evidence that Srednicki's recently-stated biasedness criterion for noninformative priors yields rankings of priors fully consistent with an information-theoretic test of Clarke, previously applied to quantum systems by Slater.

preprint1999arXiv

Hall Normalization Constants for the Bures Volumes of the n-State Quantum Systems

We report the results of certain integrations of quantum-theoretic interest, relying, in this regard, upon recently developed parameterizations of Boya et al of the n x n density matrices, in terms of squared components of the unit (n-1)-sphere and the n x n unitary matrices. Firstly, we express the normalized volume elements of the Bures (minimal monotone) metric for n = 2 and 3, obtaining thereby "Bures prior probability distributions" over the two- and three-state systems. Then, as an essential first step in extending these results to n > 3, we determine that the "Hall normalization constant" (C_{n}) for the marginal Bures prior probability distribution over the (n-1)-dimensional simplex of the n eigenvalues of the n x n density matrices is, for n = 4, equal to 71680/pi^2. Since we also find that C_{3} = 35/pi, it follows that C_{4} is simply equal to 2^{11} C_{3}/pi. (C_{2} itself is known to equal 2/pi.) The constant C_{5} is also found. It too is associated with a remarkably simple decompositon, involving the product of the eight consecutive prime numbers from 2 to 23. We also preliminarily investigate several cases, n > 5, with the use of quasi-Monte Carlo integration. We hope that the various analyses reported will prove useful in deriving a general formula (which evidence suggests will involve the Bernoulli numbers) for the Hall normalization constant for arbitrary n. This would have diverse applications, including quantum inference and universal quantum coding.

preprint1997arXiv

Asymptotic Redundancies for Universal Quantum Coding

Clarke and Barron have recently shown that the Jeffreys' invariant prior of Bayesian theory yields the common asymptotic (minimax and maximin) redundancy of universal data compression in a parametric setting. We seek a possible analogue of this result for the two-level {\it quantum} systems. We restrict our considerations to prior probability distributions belonging to a certain one-parameter family, $q(u)$, $-\infty < u < 1$. Within this setting, we are able to compute exact redundancy formulas, for which we find the asymptotic limits. We compare our quantum asymptotic redundancy formulas to those derived by naively applying the classical counterparts of Clarke and Barron, and find certain common features. Our results are based on formulas we obtain for the eigenvalues and eigenvectors of $2^n \times 2^n$ (Bayesian density) matrices, $ζ_{n}(u)$. These matrices are the weighted averages (with respect to $q(u)$) of all possible tensor products of $n$ identical $2 \times 2$ density matrices, representing the two-level quantum systems. We propose a form of {\it universal} coding for the situation in which the density matrix describing an ensemble of quantum signal states is unknown. A sequence of $n$ signals would be projected onto the dominant eigenspaces of $\ze_n(u)$.