Discrete Poincaré Lemma

This paper proves a discrete analogue of the Poincar´e lemma in the context of a discrete exterior calculus based on simplicial cochains. The proof requires the construction of a generalized cone operator, p : Ck(K) -> Ck+1(K), as the geometric cone of a simplex cannot, in general, be interpreted as a chain in the simplicial complex. The corresponding cocone operator H : Ck(K) -> Ck−1(K) can be shown to be a homotopy operator, and this yields the discrete Poincar´e lemma. The generalized cone operator is a combinatorial operator that can be constructed for any simplicial complex that can be grown by a process of local augmentation. In particular, regular triangulations and tetrahedralizations of R2 and R3 are presented, for which the discrete Poincar´e lemma is globally valid

Measurement of the Radius of Neutron Stars with High S/N Quiescent Low-mass X-ray Binaries in Globular Clusters

This paper presents the measurement of the neutron star (NS) radius using the thermal spectra from quiescent low-mass X-ray binaries (qLMXBs) inside globular clusters (GCs). Recent observations of NSs have presented evidence that cold ultra dense matter -- present in the core of NSs -- is best described by "normal matter" equations of state (EoSs). Such EoSs predict that the radii of NSs, Rns, are quasi-constant (within measurement errors, of ~10%) for astrophysically relevant masses (Mns > 0.5 Msun). The present work adopts this theoretical prediction as an assumption, and uses it to constrain a single Rns value from five qLMXB targets with available high signal-to-noise X-ray spectroscopic data. Employing a Markov-Chain Monte-Carlo approach, we produce the marginalized posterior distribution for Rns, constrained to be the same value for all five NSs in the sample. An effort was made to include all quantifiable sources of uncertainty into the uncertainty of the quoted radius measurement. These include the uncertainties in the distances to the GCs, the uncertainties due to the Galactic absorption in the direction of the GCs, and the possibility of a hard power-law spectral component for count excesses at high photon energy, which are observed in some qLMXBs in the Galactic plane. Using conservative assumptions,we found that the radius, common to the five qLMXBs and constant for a wide range of masses, lies in the low range of possible NS radii, Rns=9.1(+1.3)(-1.5) km (90%-confidence). Such a value is consistent with low-res equations of state. We compare this result with previous radius measurements of NSs from various analyses of different types of systems. In addition, we compare the spectral analyses of individual qLMXBs to previous works.Comment: Accepted to Apj. 31 pages, 17 figures, 8 table

A Discrete Geometric Optimal Control Framework for Systems with Symmetries

This paper studies the optimal motion control of mechanical systems through a discrete geometric approach. At the core of our formulation is a discrete Lagrange-d’Alembert- Pontryagin variational principle, from which are derived discrete equations of motion that serve as constraints in our optimization framework. We apply this discrete mechanical approach to holonomic systems with symmetries and, as a result, geometric structure and motion invariants are preserved. We illustrate our method by computing optimal trajectories for a simple model of an air vehicle flying through a digital terrain elevation map, and point out some of the numerical benefits that ensue

Non-monotonic density dependence of the diffusion of DNA fragments in low-salt suspensions

The high linear charge density of 20-base-pair oligomers of DNA is shown to lead to a striking non-monotonic dependence of the long-time self-diffusion on the concentration of the DNA in low-salt conditions. This generic non-monotonic behavior results from both the strong coupling between the electrostatic and solvent-mediated hydrodynamic interactions, and from the renormalization of these electrostatic interactions at large separations, and specifically from the dominance of the far-field hydrodynamic interactions caused by the strong repulsion between the DNA fragments.Comment: 4 pages, 2 figures. Physical Review E, accepted on November 24, 200

Geometric, Variational Discretization of Continuum Theories

This study derives geometric, variational discretizations of continuum theories arising in fluid dynamics, magnetohydrodynamics (MHD), and the dynamics of complex fluids. A central role in these discretizations is played by the geometric formulation of fluid dynamics, which views solutions to the governing equations for perfect fluid flow as geodesics on the group of volume-preserving diffeomorphisms of the fluid domain. Inspired by this framework, we construct a finite-dimensional approximation to the diffeomorphism group and its Lie algebra, thereby permitting a variational temporal discretization of geodesics on the spatially discretized diffeomorphism group. The extension to MHD and complex fluid flow is then made through an appeal to the theory of Euler-Poincar\'{e} systems with advection, which provides a generalization of the variational formulation of ideal fluid flow to fluids with one or more advected parameters. Upon deriving a family of structured integrators for these systems, we test their performance via a numerical implementation of the update schemes on a cartesian grid. Among the hallmarks of these new numerical methods are exact preservation of momenta arising from symmetries, automatic satisfaction of solenoidal constraints on vector fields, good long-term energy behavior, robustness with respect to the spatial and temporal resolution of the discretization, and applicability to irregular meshes