Density matrix renormalization group algorithm for Bethe lattices of spin-1/2 or spin-1 sites with Heisenberg antiferromagnetic exchange

A density matrix renormalization group (DMRG) algorithm is presented for the Bethe lattice with connectivity Z = 3 and antiferromagnetic exchange between nearest-neighbor spins s = 1/2 or 1 sites in successive generations g. The algorithm is accurate for s = 1 sites. The ground states are magnetic with spin S(g) = 2(g)s, staggered magnetization that persists for large g > 20, and short-range spin correlation functions that decrease exponentially. A finite energy gap to S > S(g) leads to a magnetization plateau in the extended lattice. Closely similar DMRG results for s = 1/2 and 1 are interpreted in terms of an analytical three-site model.

The Ginibre Ensemble and Gaussian Analytic Functions

We show that as n changes, the characteristic polynomial of the n x n random matrix with i.i.d. complex Gaussian entries can be described recursively through a process analogous to Polya's urn scheme. As a result, we get a random analytic function in the limit, which is given by a mixture of Gaussian analytic functions. This suggests another reason why the zeros of Gaussian analytic functions and the Ginibre ensemble exhibit similar local repulsion, but different global behavior. Our approach gives new explicit formulas for the limiting analytic function.

Automated Computation of the Fundamental Matrix for Vision-Based Construction Site Applications

Estimating the fundamental matrix (F), to determine the epipolar geometry between a pair of images or video frames, is a basic step for a wide variety of vision-based functions used in construction operations, such as camera-pair calibration, automatic progress monitoring, and 3D reconstruction. Currently, robust methods (e.g., SIFT + normalized eight-point algorithm + RANSAC) are widely used in the construction community for this purpose. Although they can provide acceptable accuracy, the significant amount of required computational time impedes their adoption in real-time applications, especially video data analysis with many frames per second. Aiming to overcome this limitation, this paper presents and evaluates the accuracy of a solution to find F by combining the use of two speedy and consistent methods: SURF for the selection of a robust set of point correspondences and the normalized eight-point algorithm. This solution is tested extensively on construction site image pairs including changes in viewpoint, scale, illumination, rotation, and moving objects. The results demonstrate that this method can be used for real-time applications (5 image pairs per second with the resolution of 640 × 480) involving scenes of the built environment.

Correlation effects for semiconducting single-wall carbon nanotubes: A density matrix renormalization group study

In this paper we report the applicability of the density matrix renormalization group (DMRG) approach to the cylindrical single wall carbon nanotube (SWCN) for the purpose of its correlation effect. By applying the DMRG approach to the t+U+V model, with t and V being the hopping and Coulomb energies between the nearest neighboring sites, respectively, and U the on-site Coulomb energy, we calculate the phase diagram for the SWCN with chiral numbers (n(1)=3, n(2)=2), which reflects the competition between the correlation energy U and V. Within reasonable parameter ranges, we investigate possible correlated ground states, the lowest excitations, and the corresponding correlation functions in which the connection with the excitonic insulator is particularly addressed.

A transfer matrix approach to conductance in quantum waveguides

A transfer matrix approach is presented for the study of electron conduction in an arbitrarily shaped cavity structure embedded in a quantum wire. Using the boundary conditions for wave functions, the transfer matrix at an interface with a discontinuous potential boundary is obtained for the first time. The total transfer matrix is calculated by multiplication of the transfer matrix for each segment of the structure as well as numerical integration of coupled second-order differential equations. The proposed method is applied to the evaluation of the conductance and the electron probability density in several typical cavity structures. The effect of the geometrical features on the electron transmission is discussed in detail. In the numerical calculations, the method is found to be more efficient than most of the other methods in the literature and the results are found to be in excellent agreement with those obtained by the recursive Green's function method.

Restudy on the wave functions of Gamma(nS) states in the light-front quark model and the radiative decays of Gamma(nS) -> eta(b) + gamma

The light-front quark model has been applied to calculate the transition matrix elements of heavy hadron decays. However, it is noted that using the traditional wave functions of the light-front quark model given in the literature, the theoretically determined decay constants of the Gamma(nS) obviously contradict the data. This implies that the wave functions must be modified. Keeping the orthogonality among the nS states and fitting their decay constants, we obtain a series of the wave functions for Gamma(nS). Based on these wave functions and by analogy with the hydrogen atom, we suggest a modified analytical form for the Gamma(nS) wave functions. Using the modified wave functions, the obtained decay constants are close to the experimental data. Then we calculate the rates of radiative decays of Gamma(nS) -> eta(b) + gamma. Our predictions are consistent with the experimental data on decays Gamma(3S) -> eta(b) + gamma within the theoretical and experimental errors.

Financial modelling with 2-EPT probability density functions

The class of all Exponential-Polynomial-Trigonometric (EPT) functions is classical and equal to the Euler-d’Alembert class of solutions of linear differential equations with constant coefficients. The class of non-negative EPT functions defined on [0;1) was discussed in Hanzon and Holland (2010) of which EPT probability density functions are an important subclass. EPT functions can be represented as ceAxb, where A is a square matrix, b a column vector and c a row vector where the triple (A; b; c) is the minimal realization of the EPT function. The minimal triple is only unique up to a basis transformation. Here the class of 2-EPT probability density functions on R is defined and shown to be closed under a variety of operations. The class is also generalised to include mixtures with the pointmass at zero. This class coincides with the class of probability density functions with rational characteristic functions. It is illustrated that the Variance Gamma density is a 2-EPT density under a parameter restriction. A discrete 2-EPT process is a process which has stochastically independent 2-EPT random variables as increments. It is shown that the distribution of the minimum and maximum of such a process is an EPT density mixed with a pointmass at zero. The Laplace Transform of these distributions correspond to the discrete time Wiener-Hopf factors of the discrete time 2-EPT process. A distribution of daily log-returns, observed over the period 1931-2011 from a prominent US index, is approximated with a 2-EPT density function. Without the non-negativity condition, it is illustrated how this problem is transformed into a discrete time rational approximation problem. The rational approximation software RARL2 is used to carry out this approximation. The non-negativity constraint is then imposed via a convex optimisation procedure after the unconstrained approximation. Sufficient and necessary conditions are derived to characterise infinitely divisible EPT and 2-EPT functions. Infinitely divisible 2-EPT density functions generate 2-EPT Lévy processes. An assets log returns can be modelled as a 2-EPT Lévy process. Closed form pricing formulae are then derived for European Options with specific times to maturity. Formulae for discretely monitored Lookback Options and 2-Period Bermudan Options are also provided. Certain Greeks, including Delta and Gamma, of these options are also computed analytically. MATLAB scripts are provided for calculations involving 2-EPT functions. Numerical option pricing examples illustrate the effectiveness of the 2-EPT approach to financial modelling.

Stochastic comparability and dual q-functions

In this paper we discuss the relationship and characterization of stochastic comparability, duality, and Feller–Reuter–Riley transition functions which are closely linked with each other for continuous time Markov chains. A necessary and sufficient condition for two Feller minimal transition functions to be stochastically comparable is given in terms of their density q-matrices only. Moreover, a necessary and sufficient condition under which a transition function is a dual for some stochastically monotone q-function is given in terms of, again, its density q-matrix. Finally, for a class of q-matrices, the necessary and sufficient condition for a transition function to be a Feller–Reuter–Riley transition function is also given.

Applications of Feller-Reuter-Riley transition functions

A Feller–Reuter–Riley function is a Markov transition function whose corresponding semigroup maps the set of the real-valued continuous functions vanishing at infinity into itself. The aim of this paper is to investigate applications of such functions in the dual problem, Markov branching processes, and the Williams-matrix. The remarkable property of a Feller–Reuter–Riley function is that it is a Feller minimal transition function with a stable q-matrix. By using this property we are able to prove that, in the theory of branching processes, the branching property is equivalent to the requirement that the corresponding transition function satisfies the Kolmogorov forward equations associated with a stable q-matrix. It follows that the probabilistic definition and the analytic definition for Markov branching processes are actually equivalent. Also, by using this property, together with the Resolvent Decomposition Theorem, a simple analytical proof of the Williams' existence theorem with respect to the Williams-matrix is obtained. The close link between the dual problem and the Feller–Reuter–Riley transition functions is revealed. It enables us to prove that a dual transition function must satisfy the Kolmogorov forward equations. A necessary and sufficient condition for a dual transition function satisfying the Kolmogorov backward equations is also provided.

B-spline methods in R-matrix theory for scattering in two-electron systems

The use of B-spline basis sets in R-matrix theory for scattering processes has been investigated. In the present approach a B-spline basis is used for the description of the inner region, which is matched to the physical outgoing wavefunctions by the R-matrix. Using B-splines, continuum basis functions can be determined easily, while pseudostates can be included naturally. The accuracy for low-energy scattering processes is demonstrated by calculating inelastic scattering cross sections for e colliding on H. Very good agreement with other calculations has been obtained. Further extensions of the codes to quasi two-electron systems and general atoms are discussed as well as the application to (multi) photoionization.

Giving structure to the biofilm matrix: an overview of individual strategies and emerging common themes

Biofilms are communities of microbial cells that underpin diverse processes including sewage bioremediation, plant growth promotion, chronic infections and industrial biofouling. The cells resident in the biofilm are encased within a self-produced exopolymeric matrix that commonly comprises lipids, proteins that frequently exhibit amyloid-like properties, eDNA and exopolysaccharides. This matrix fulfils a variety of functions for the community, from providing structural rigidity and protection from the external environment to controlling gene regulation and nutrient adsorption. Critical to the development of novel strategies to control biofilm infections, or the capability to capitalize on the power of biofilm formation for industrial and biotechnological uses, is an in-depth knowledge of the biofilm matrix. This is with respect to the structure of the individual components, the nature of the interactions between the molecules and the three-dimensional spatial organization. We highlight recent advances in the understanding of the structural and functional role that carbohydrates and proteins play within the biofilm matrix to provide three-dimensional architectural integrity and functionality to the biofilm community. We highlight, where relevant, experimental techniques that are allowing the boundaries of our understanding of the biofilm matrix to be extended using Escherichia coli, Staphylococcus aureus, Vibrio cholerae, and Bacillus subtilis as exemplars.

Photoionization cross sections for O-like S IX: a Breit-Pauli R-matrix calculation

In this paper we present photoionization cross sections for the lowest five states of O-like S IX (1s(2)2s(2)2p(4) P-3(0,1,2), D-1(2), S-1(0)). The relativistic Breit-Pauli R-matrix codes were utilized including all terms of the 2s(2)2p(3), 2s2p(4), 2p(5), 2s(2)2p(2)3s, 3p, 3d and 2s2p(3)3s, 3p, 3d configurations in the expansion of the collision wavefunction for S X. It was also found that to achieve convergence of the low-lying energy separations of the target levels, an additional 21 configuration functions needed to be included in the configuration interaction expansion, incorporating two-electron excitations from the 2s and 2p shells to the 3s, 3p and 3d shells. The present work thus constitutes the most sophisticated photoionization evaluation for ground and metastable levels of the S IX ion. Direct comparisons have been made with the only available data found on the OPEN-ADAS database between level resolved contributions of the spectrum. This comparison for the background cross section exhibits excellent agreement at all photon energies for each partial photoionization cross section contribution investigated. Finally, the autoionizing bound states arising from numerous open channels have also been investigated and identified using the QB approach, a procedure for analyzing resonances in atomic and molecular collision theory which exploits the analytic properties of R-matrix theory. Major Rydberg resonance series are also presented and tabulated for the dominant linewidths considered.

Evolution of matrix and bone gamma-carboxyglutamic acid proteins in vertebrates

The evolution of calcified tissues is a defining feature in vertebrate evolution. Investigating the evolution of proteins involved in tissue calcification should help elucidate how calcified tissues have evolved. The purpose of this study was to collect and compare sequences of matrix and bone γ-carboxyglutamic acid proteins (MGP and BGP, respectively) to identify common features and determine the evolutionary relationship between MGP and BGP. Thirteen cDNAs and genes were cloned using standard methods or reconstructed through the use of comparative genomics and data mining. These sequences were compared with available annotated sequences (a total of 48 complete or nearly complete sequences, 28 BGPs and 20 MGPs) have been identified across 32 different species (representing most classes of vertebrates), and evolutionarily conserved features in both MGP and BGP were analyzed using bioinformatic tools and the Tree-Puzzle software. We propose that: 1) MGP and BGP genes originated from two genome duplications that occurred around 500 and 400 million years ago before jawless and jawed fish evolved, respectively; 2) MGP appeared first concomitantly with the emergence of cartilaginous structures, and BGP appeared thereafter along with bony structures; and 3) BGP derives from MGP. We also propose a highly specific pattern definition for the Gla domain of BGP and MGP. Previous Section Next Section BGP1 (bone Gla protein or osteocalcin) and MGP (matrix Gla protein) belong to the growing family of vitamin K-dependent (VKD) proteins, the members of which are involved in a broad range of biological functions such as skeletogenesis and bone maintenance (BGP and MGP), hemostasis (prothrombin, clotting factors VII, IX, and X, and proteins C, S, and Z), growth control (gas6), and potentially signal transduction (proline-rich Gla proteins 1 and 2). VKD proteins are characterized by the presence of several Gla residues resulting from the post-translational vitamin K-dependent γ-carboxylation of specific glutamates, through which they can bind to calcium-containing mineral such as hydroxyapatite. To date, VKD proteins have only been clearly identified in vertebrates (1) although the presence of a γ-glutamyl carboxylase has been reported in the fruit fly Drosophila melanogaster (2) and in marine snails belonging to the genus Conus (3). Gla residues have also been found in neuropeptides from Conus venoms (4), suggesting a wider prevalence of γ-carboxylation.

Non-negative matrix decomposition approaches to frequency domain analysis of music audio signals

On étudie l’application des algorithmes de décomposition matricielles tel que la Factorisation Matricielle Non-négative (FMN), aux représentations fréquentielles de signaux audio musicaux. Ces algorithmes, dirigés par une fonction d’erreur de reconstruction, apprennent un ensemble de fonctions de base et un ensemble de coef- ficients correspondants qui approximent le signal d’entrée. On compare l’utilisation de trois fonctions d’erreur de reconstruction quand la FMN est appliquée à des gammes monophoniques et harmonisées: moindre carré, divergence Kullback-Leibler, et une mesure de divergence dépendente de la phase, introduite récemment. Des nouvelles méthodes pour interpréter les décompositions résultantes sont présentées et sont comparées aux méthodes utilisées précédemment qui nécessitent des connaissances du domaine acoustique. Finalement, on analyse la capacité de généralisation des fonctions de bases apprises par rapport à trois paramètres musicaux: l’amplitude, la durée et le type d’instrument. Pour ce faire, on introduit deux algorithmes d’étiquetage des fonctions de bases qui performent mieux que l’approche précédente dans la majorité de nos tests, la tâche d’instrument avec audio monophonique étant la seule exception importante.

Negative Fukui functions: new insights based on electronegativity equalization

The occurrence of negative values for Fukui functions was studied through the electronegativity equalization method. Using algebraic relations between Fukui functions and different other conceptual DFT quantities on the one hand and the hardness matrix on the other hand, expressions were obtained for Fukui functions for several archetypical small molecules. Based on EEM calculations for large molecular sets, no negative Fukui functions were found