Least square projection: A fast high-precision multidimensional projection technique and its application to document mapping

The problem of projecting multidimensional data into lower dimensions has been pursued by many researchers due to its potential application to data analyses of various kinds. This paper presents a novel multidimensional projection technique based on least square approximations. The approximations compute the coordinates of a set of projected points based on the coordinates of a reduced number of control points with defined geometry. We name the technique Least Square Projections ( LSP). From an initial projection of the control points, LSP defines the positioning of their neighboring points through a numerical solution that aims at preserving a similarity relationship between the points given by a metric in mD. In order to perform the projection, a small number of distance calculations are necessary, and no repositioning of the points is required to obtain a final solution with satisfactory precision. The results show the capability of the technique to form groups of points by degree of similarity in 2D. We illustrate that capability through its application to mapping collections of textual documents from varied sources, a strategic yet difficult application. LSP is faster and more accurate than other existing high-quality methods, particularly where it was mostly tested, that is, for mapping text sets.

3-manifolds in Euclidean space from a contact viewpoint

We study the geometry of 3-manifolds generically embedded in R(n) by means of the analysis of the singularities of the distance-squared and height functions on them. We describe the local structure of the discriminant (associated to the distribution of asymptotic directions), the ridges and the flat ridges.

On the lower bound on the exchange-correlation energy in two dimensions

We study the properties of the lower bound on the exchange-correlation energy in two dimensions. First we review the derivation of the bound and show how it can be written in a simple density-functional form. This form allows an explicit determination of the prefactor of the bound and testing its tightness. Next we focus on finite two-dimensional systems and examine how their distance from the bound depends on the system geometry. The results for the high-density limit suggest that a finite system that comes as close as possible to the ultimate bound on the exchange-correlation energy has circular geometry and a weak confining potential with a negative curvature. (c) 2009 Elsevier B.V. All rights reserved.

Shot Noise in a Spin-Diode Geometry

We apply the master equation technique to calculate shot noise in a system composed of single level quantum dot attached to a normal metal lead and to a ferromagnetic lead (NM-QD-FM). It is known that this system operates as a spin-diode, giving unpolarized currents for forward bias and polarized current for reverse bias. This effect is observed when only one electron can tunnel at a time through the dot, due to the strong intradot Coulomb interaction. We find that the shot noise also presents a signature of this spin-diode effect, with a super-Poissonian shot noise for forward and a sub-Poissonian shot noise for reverse bias voltages. The shot noise thus can provide further experimental evidence of the spin-rectification in the NM-QD-FM geometry.

SRTM resample with short distance-low nugget kriging

The shuttle radar topography mission (SRTM), was flow on the space shuttle Endeavour in February 2000, with the objective of acquiring a digital elevation model of all land between 60 degrees north latitude and 56 degrees south latitude, using interferometric synthetic aperture radar (InSAR) techniques. The SRTM data are distributed at horizontal resolution of 1 arc-second (similar to 30m) for areas within the USA and at 3 arc-second (similar to 90m) resolution for the rest of the world. A resolution of 90m can be considered suitable for the small or medium-scale analysis, but it is too coarse for more detailed purposes. One alternative is to interpolate the SRTM data at a finer resolution; it will not increase the level of detail of the original digital elevation model (DEM), but it will lead to a surface where there is the coherence of angular properties (i.e. slope, aspect) between neighbouring pixels, which is an important characteristic when dealing with terrain analysis. This work intents to show how the proper adjustment of variogram and kriging parameters, namely the nugget effect and the maximum distance within which values are used in interpolation, can be set to achieve quality results on resampling SRTM data from 3"" to 1"". We present for a test area in western USA, which includes different adjustment schemes (changes in nugget effect value and in the interpolation radius) and comparisons with the original 1"" model of the area, with the national elevation dataset (NED) DEMs, and with other interpolation methods (splines and inverse distance weighted (IDW)). The basic concepts for using kriging to resample terrain data are: (i) working only with the immediate neighbourhood of the predicted point, due to the high spatial correlation of the topographic surface and omnidirectional behaviour of variogram in short distances; (ii) adding a very small random variation to the coordinates of the points prior to interpolation, to avoid punctual artifacts generated by predicted points with the same location than original data points and; (iii) using a small value of nugget effect, to avoid smoothing that can obliterate terrain features. Drainages derived from the surfaces interpolated by kriging and by splines have a good agreement with streams derived from the 1"" NED, with correct identification of watersheds, even though a few differences occur in the positions of some rivers in flat areas. Although the 1"" surfaces resampled by kriging and splines are very similar, we consider the results produced by kriging as superior, since the spline-interpolated surface still presented some noise and linear artifacts, which were removed by kriging.

PACKMOL: A Package for Building Initial Configurations for Molecular Dynamics Simulations

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system`s overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules Must keep safe pairwise distances. Obtaining Such a molecular arrangement can be considered it packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules Must be greater than some specified tolerance. We have developed a code able to pack millions of atoms. grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straight forward. In addition. different atoms belonging to the same molecule may also be restricted to different spatial regions, in Such a way that more ordered molecular arrangements call be built, as micelles. lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies front a few seconds to a few Minutes in a personal Computer. The input files are simple and Currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and call be downloaded front http://www.ime.unicamp.br/similar to martinez/packmol/. (C) 2009 Wiley Periodicals. Inc. J Comput Chem 30: 2157-2164, 2009

Perfect Simulation of a Coupling Achieving the (d)over-bar-distance Between Ordered Pairs of Binary Chains of Infinite Order

We explicitly construct a stationary coupling attaining Ornstein`s (d) over bar -distance between ordered pairs of binary chains of infinite order. Our main tool is a representation of the transition probabilities of the coupled bivariate chain of infinite order as a countable mixture of Markov transition probabilities of increasing order. Under suitable conditions on the loss of memory of the chains, this representation implies that the coupled chain can be represented as a concatenation of i.i.d. sequences of bivariate finite random strings of symbols. The perfect simulation algorithm is based on the fact that we can identify the first regeneration point to the left of the origin almost surely.

Decay of geometry for Fibonacci critical covering maps of the circle

We study the growth of Df `` (f(c)) when f is a Fibonacci critical covering map of the circle with negative Schwarzian derivative, degree d >= 2 and critical point c of order l > 1. As an application we prove that f exhibits exponential decay of geometry if and only if l <= 2, and in this case it has an absolutely continuous invariant probability measure, although not satisfying the so-called Collet-Eckmann condition. (C) 2009 Elsevier Masson SAS. All rights reserved.

Periodic geodesics and geometry of compact Lorentzian manifolds with a Killing vector field

We study the geometry and the periodic geodesics of a compact Lorentzian manifold that has a Killing vector field which is timelike somewhere. Using a compactness argument for subgroups of the isometry group, we prove the existence of one timelike non self-intersecting periodic geodesic. If the Killing vector field is nowhere vanishing, then there are at least two distinct periodic geodesics; as a special case, compact stationary manifolds have at least two periodic timelike geodesics. We also discuss some properties of the topology of such manifolds. In particular, we show that a compact manifold M admits a Lorentzian metric with a nowhere vanishing Killing vector field which is timelike somewhere if and only if M admits a smooth circle action without fixed points.

COMPARISON RESULTS FOR CONJUGATE AND FOCAL POINTS IN SEMI-RIEMANNIAN GEOMETRY VIA MASLOV INDEX

We prove an estimate on the difference of Maslov indices relative to the choice of two distinct reference Lagrangians of a continuous path in the Lagrangian Grassmannian of a symplectic space. We discuss some applications to the study of conjugate and focal points along a geodesic in a semi-Riemannian manifold.

Understanding the Effect of Adsorption Geometry over Substrate Selectivity in the Surface-Enhanced Raman Scattering Spectra of Simazine and Atrazine

This paper studies the selectivity of Well-defined Au and Ag nanostructures as substrates for the SERS, (surface-enhanced Raman scattering) detection of simazine (6-chloro-N,N`-diethyl-1,3,5-triazine-2,4-diamine) and atrazine (6-chloro-N-ethyl-N`-isopropyl-1,3,5-triazine-2,4-diamine). Our data showed that simazine and atrazine displayed similar SERS spectra when the Au was employed as substrate. Conversely, distinct SERS signatures were obtained upon the utilization of Ag substrates. Density functional theory (DFT) calculations and vibrational assignments suggested that, while simazine and atrazine adsorbed on Au via the N3 position of the triazine ring, simazine adsorbed on Ag via N3 and atrazine via N5. The results presented herein demonstrated that the adsorption geometry of analyte molecules can play a central role over substrate selectivity in SERS, which is particularly important in applications involving ultrasensitive analysis of mixtures containing structurally similar molecules.