Tensor LRR and sparse coding-based subspace clustering

Subspace clustering groups a set of samples from a union of several linear subspaces into clusters, so that the samples in the same cluster are drawn from the same linear subspace. In the majority of the existing work on subspace clustering, clusters are built based on feature information, while sample correlations in their original spatial structure are simply ignored. Besides, original high-dimensional feature vector contains noisy/redundant information, and the time complexity grows exponentially with the number of dimensions. To address these issues, we propose a tensor low-rank representation (TLRR) and sparse coding-based (TLRRSC) subspace clustering method by simultaneously considering feature information and spatial structures. TLRR seeks the lowest rank representation over original spatial structures along all spatial directions. Sparse coding learns a dictionary along feature spaces, so that each sample can be represented by a few atoms of the learned dictionary. The affinity matrix used for spectral clustering is built from the joint similarities in both spatial and feature spaces. TLRRSC can well capture the global structure and inherent feature information of data, and provide a robust subspace segmentation from corrupted data. Experimental results on both synthetic and real-world data sets show that TLRRSC outperforms several established state-of-the-art methods.

Heterogeneous tensor decomposition for clustering via manifold optimization

Tensor clustering is an important tool that exploits intrinsically rich structures in real-world multiarray or Tensor datasets. Often in dealing with those datasets, standard practice is to use subspace clustering that is based on vectorizing multiarray data. However, vectorization of tensorial data does not exploit complete structure information. In this paper, we propose a subspace clustering algorithm without adopting any vectorization process. Our approach is based on a novel heterogeneous Tucker decomposition model taking into account cluster membership information. We propose a new clustering algorithm that alternates between different modes of the proposed heterogeneous tensor model. All but the last mode have closed-form updates. Updating the last mode reduces to optimizing over the multinomial manifold for which we investigate second order Riemannian geometry and propose a trust-region algorithm. Numerical experiments show that our proposed algorithm compete effectively with state-of-the-art clustering algorithms that are based on tensor factorization.

CMB and LSS constraints on a single-field model of inflation

A new inflationary scenario whose exponential potential V (Phi) has a quadratic dependence on the field Phi in addition to the standard linear term is confronted with the five-year observations of the Wilkinson-Microwave Anisotropy Probe and the Sloan Digital Sky Survey data. The number of e-folds (N), the ratio of tensor-to-scalar perturbations (r), the spectral scalar index of the primordial power spectrum (n(s)) and its running (dn(s)/d ln k) depend on the dimensionless parameter a multiplying the quadratic term in the potential. In the limit a. 0 all the results of the exponential potential are fully recovered. For values of alpha not equal 0, we find that the model predictions are in good agreement with the current observations of the Cosmic Microwave Background (CMB) anisotropies and Large-Scale Structure (LSS) in the Universe. Copyright (C) EPLA, 2008.

Quantized fields and gravitational particle creation in f (R) expanding universes

The problem of cosmological particle creation for a spatially flat, homogeneous and isotropic universes is discussed in the context of f (R) theories of gravity. Different from cosmological models based on general relativity theory, it is found that a conformal invariant metric does not forbid the creation of massless particles during the early stages (radiation era) of the universe. (C) 2010 Elsevier B.V. All rights reserved.

Integrated marine gravity field in the Brazilian coast from altimeter-derived sea surface gradient and shipborne gravity

Sea surface gradients derived from the Geosat and ERS-1 satellite altimetry geodetic missions were integrated with marine gravity data from the National Geophysical Data Center and Brazilian national surveys. Using the least squares collocation method, models of free-air gravity anomaly and geoid height were calculated for the coast of Brazil with a resolution of 2` x 2`. The integration of satellite and shipborne data showed better statistical results in regions near the coast than using satellite data only, suggesting an improvement when compared to the state-of-the-art global gravity models. Furthermore, these results were obtained with considerably less input information than was used by those reference models. The least squares collocation presented a very low content of high-frequency noise in the predicted gravity anomalies. This may be considered essential to improve the high resolution representation of the gravity field in regions of ocean-continent transition. (C) 2010 Elsevier Ltd. All rights reserved.

Magnetotelluric deep soundings, gravity and geoid in the south Sao Francisco craton: Geophysical indicators of cratonic lithosphere rejuvenation and crustal underplating

In the south Sao Francisco craton a circular and 8-m amplitude geoid anomaly coincides with the outcropping terrain of an Archean-Paleoproterozoic basement. Broadband magnetotelluric (MT) data inversions of two radial profiles within the positive geoid and Bouguer gravity anomaly yield geo-electrical crustal sections, whereby the lower crust is locally more conductive (10 to 100 Omega m) in spatial coincidence with a denser lower crust modeled by the gravity data. This anomalous lower crust may have resulted from magmatic underplating, associated with Mesoarchean and Proterozoic episodes of tholeiitic dike intrusion. Long-period MT soundings reveal a low electrical resistivity mantle (20 to 200 Omega m) from depths beyond 120 km. Forward geoid modeling, using the scope of the low electrical resistivity region within the mantle as a constraint, entails a density increase (40 to 50 kg/m(3)) possibly due to Fe enrichment of mantle minerals. However, this factor alone does not explain the observed resistivity. A supplemented presence of small amounts of percolated carbonatite melting (similar to 0.005 vol.%), dissolved water and enhanced oxygen fugacity within the peridotitic mantle are viable agents that could explain the less resistive upper mantle. We propose that metasomatic processes confined in the sub-continental lithospheric mantle foster the conditions for a low degree melting with variable CO(2), H(2)O and Fe content. Even though the precise age of this metasomatism is unknown it might be older than the Early Cretaceous based on the evidence that a high-degree of melting in a lithospheric mantle impregnated with carbonatites originated the tholeiitic dike intrusions dispersed from the southeastern border of the Sao Francisco craton, during the onset of the lithosphere extension and break-up of the western Gondwana. The proxies are the NE Parana and Espinhaco (130 Ma, Ar/Ar ages) tholeiitic dikes, which contain (similar to 3%) carbonatites in their composition. The occurrence of a positive geoid anomaly (+ 10 m) and pre-tholeiites (age > 138 Ma), carbonatites and kimberlites along the west African continental margin (Angola and Namibia) reinforces the presumed age of the Sao Francisco-Congo craton rejuvenation to be prior to its fragmentation in the Lower Cretaceous. (C) 2010 Elsevier B.V. All rights reserved.

Reappraisal of the effective elastic thickness for the sub-Andes using 3-D finite element flexural modelling, gravity and geological constraints

P>Estimates of effective elastic thickness (T(e)) for the western portion of the South American Plate using, independently, forward flexural modelling and coherence analysis, suggest different thermomechanical properties for the same continental lithosphere. We present a review of these T(e) estimates and carry out a critical reappraisal using a common methodology of 3-D finite element method to solve a differential equation for the bending of a thin elastic plate. The finite element flexural model incorporates lateral variations of T(e) and the Andes topography as the load. Three T(e) maps for the entire Andes were analysed: Stewart & Watts (1997), Tassara et al. (2007) and Perez-Gussinye et al. (2007). The predicted flexural deformation obtained for each T(e) map was compared with the depth to the base of the foreland basin sequence. Likewise, the gravity effect of flexurally induced crust-mantle deformation was compared with the observed Bouguer gravity. T(e) estimates using forward flexural modelling by Stewart & Watts (1997) better predict the geological and gravity data for most of the Andean system, particularly in the Central Andes, where T(e) ranges from greater than 70 km in the sub-Andes to less than 15 km under the Andes Cordillera. The misfit between the calculated and observed foreland basin subsidence and the gravity anomaly for the Maranon basin in Peru and the Bermejo basin in Argentina, regardless of the assumed T(e) map, may be due to a dynamic topography component associated with the shallow subduction of the Nazca Plate beneath the Andes at these latitudes.

Programs to compute magnetization to density ratio and the magnetization inclination from 3-D gravity and magnetic anomalies

Vector field formulation based on the Poisson theorem allows an automatic determination of rock physical properties (magnetization to density ratio-MDR-and the magnetization inclination-MI) from combined processing of gravity and magnetic geophysical data. The basic assumptions (i.e., Poisson conditions) are: that gravity and magnetic fields share common sources, and that these sources have a uniform magnetization direction and MDR. In addition, the previously existing formulation was restricted to profile data, and assumed sufficiently elongated (2-D) sources. For sources that violate Poisson conditions or have a 3-D geometry, the apparent values of MDR and MI that are generated in this way have an unclear relationship to the actual properties in the subsurface. We present Fortran programs that estimate MDR and MI values for 3-D sources through processing of gridded gravity and magnetic data. Tests with simple geophysical models indicate that magnetization polarity can be successfully recovered by MDR-MI processing, even in cases where juxtaposed bodies cannot be clearly distinguished on the basis of anomaly data. These results may be useful in crustal studies, especially in mapping magnetization polarity from marine-based gravity and magnetic data. (c) 2007 Elsevier Ltd. All rights reserved.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

Instabilities in thermal gravity with a cosmological constant

It is shown that in quantum gravity at finite temperature, the effective potential evaluated in the tadpole approximation can have a local minimum below a certain critical temperature. However, when the leading higher order thermal loop corrections are included, one finds that no static solution exists at high temperature. (C) 2008 Elsevier B.V. All rights reserved.

Spacetime noncommutativity with a bifermionic parameter

We consider a Moyal plane and propose to make the noncommutativity parameter Theta(mu nu) bifermionic, i.e. composed of two fermionic (Grassmann odd) parameters. The Moyal product then contains a finite number of derivatives, which avoid the difficulties of the standard approach. As an example, we construct a two-dimensional noncommutative field theory model based on the Moyal product with a bifermionic parameter and show that it has a locally conserved energy-momentum tensor. The model has no problem with the canonical quantization and appears to be renormalizable.

Pseudoclassical description of scalar particle in non-abelian background and path-integral representations

Path-integral representations for a scalar particle propagator in non-Abelian external backgrounds are derived. To this aim, we generalize the procedure proposed by Gitman and Schvartsman of path-integral construction to any representation of SU(N) given in terms of antisymmetric generators. And for arbitrary representations of SU(N), we present an alternative construction by means of fermionic coherent states. From the path-integral representations we derive pseudoclassical actions for a scalar particle placed in non-Abelian backgrounds. These actions are classically analyzed and then quantized to prove their consistency.

Field on Poincare Group and Quantum Description of Orientable Objects

We propose an approach to the quantum-mechanical description of relativistic orientable objects. It generalizes Wigner`s ideas concerning the treatment of nonrelativistic orientable objects (in particular, a nonrelativistic rotator) with the help of two reference frames (space-fixed and body-fixed). A technical realization of this generalization (for instance, in 3+1 dimensions) amounts to introducing wave functions that depend on elements of the Poincar, group G. A complete set of transformations that test the symmetries of an orientable object and of the embedding space belongs to the group I =GxG. All such transformations can be studied by considering a generalized regular representation of G in the space of scalar functions on the group, f(x,z), that depend on the Minkowski space points xaG/Spin(3,1) as well as on the orientation variables given by the elements z of a matrix ZaSpin(3,1). In particular, the field f(x,z) is a generating function of the usual spin-tensor multi-component fields. In the theory under consideration, there are four different types of spinors, and an orientable object is characterized by ten quantum numbers. We study the corresponding relativistic wave equations and their symmetry properties.

Scalar form factors and nuclear interactions

The scalar-isoscalar term in the two-pion exchange NN potential is abnormally large and does not respect the hierarchy of effects predicted by chiral perturbation theory. We argue that this anomaly is associated with non-perturbative effects, which are also present in the pi N scalar form factor.