Constraining the Mass Profiles of Stellar Systems: Schwarzschild Modeling of Discrete Velocity Datasets

(ABRIDGED) We present a new Schwarzschild orbit-superposition code designed to model discrete datasets composed of velocities of individual kinematic tracers in a dynamical system. This constitutes an extension of previous implementations that can only address continuous data in the form of (the moments of) velocity distributions, thus avoiding potentially important losses of information due to data binning. Furthermore, the code can handle any combination of available velocity components, i.e., only line-of-sight velocities, only proper motions, or a combination of both. It can also handle a combination of discrete and continuous data. The code finds the distribution function (DF, a function of the three integrals of motion E, Lz, and I3) that best reproduces the available kinematic and photometric observations in a given axisymmetric gravitational potential. The fully numerical approach ensures considerable freedom on the form of the DF f(E,Lz,I3). This allows a very general modeling of the orbital structure, thus avoiding restrictive assumptions about the degree of (an)isotropy of the orbits. We describe the implementation of the discrete code and present a series of tests of its performance based on the modeling of simulated datasets generated from a known DF. We find that the discrete Schwarzschild code recovers the original orbital structure, M/L ratios, and inclination of the input datasets to satisfactory accuracy, as quantified by various statistics. The code will be valuable, e.g., for modeling stellar motions in Galactic globular clusters, and those of individual stars, planetary nebulae, or globular clusters in nearby galaxies. This can shed new light on the total mass distributions of these systems, with central black holes and dark matter halos being of particular interest.Comment: ApJ, in press; 51 pages, 11 figures; manuscript revised following comments by refere

Matching the Hagedorn mass spectrum with Lattice QCD results

Based on recent Lattice QCD (LQCD) results obtained at finite temperature, we discuss modeling of the hadronic phase of QCD in the framework of Hadron Resonance Gas (HRG) with discrete and continuous mass spectra. We focus on fluctuations of conserved charges, and show how a common limiting temperature can be used to constrain the Hagedorn exponential mass spectrum in different sectors of quantum number, through a matching of HRG and LQCD. For strange baryons, the extracted spectra are found to be consistent with all known and expected states listed by the Particle Data Group (PDG). The strange-mesonic sector, however, requires additional states in the intermediate mass range beyond that embodied in the database.Comment: 8 pages, 4 figures, version to appear in phys. rev.

X-ray Constraints on the Intrinsic Shape of the Lenticular Galaxy NGC 1332

We have analyzed ROSAT PSPC X-ray data of the optically elongated S0 galaxy NGC 1332 with the purposes of constraining the intrinsic shape of its underlying mass and presenting a detailed investigation of the uncertainties resulting from the assumptions underlying this type of analysis. The X-ray isophotes are elongated with ellipticity $0.10 - 0.27$ (90% confidence) for semi-major axes 75\arcsec -90\arcsec and have orientations consistent with the optical isophotes (ellipticity $\sim 0.43$). The spectrum is poorly constrained by the PSPC data and cannot rule out sizeable radial temperature gradients or an emission component due to discrete sources equal in magnitude to the hot gas. Using (and clarifying) the "geometric test" for dark matter, we determined that the hypothesis that mass-traces-light is not consistent with the X-ray data at 68% confidence and marginally consistent at 90% confidence independent of the gas temperature profile. Detailed modeling gives constraints on the ellipticity of the underlying mass of \epsilon_{mass} = 0.47 - 0.72 (0.31 - 0.83) at 68% (90%) confidence for isothermal and polytropic models. The total mass of the isothermal models within a=43.6 kpc (D = 20h^{-1}_{80} Mpc) is M_{tot} = (0.38 - 1.7) \times 10^{12}M_{\sun} (90% confidence) corresponding to total blue mass-to-light ratio \Upsilon_B = (31.9 - 143) \Upsilon_{\sun}. Similar results are obtained when the dark matter is fit directly using the known distributions of the stars and gas. When possible rotation of the gas and emission from discrete sources are included flattened mass distributions are still required, although the constraints on \epsilon_{mass}$, but not the total mass, are substantially weakened.Comment: 45 pages (figures missing), PostScript, to appear in ApJ on January 20, 199

Achernar: Rapid Polarization Variability as Evidence of Photospheric and Circumstellar Activity

We present the results of a high accuracy ($\sigma \approx 0.005$) polarization monitoring of the Be Star Achernar that was carried out between July 7th and November 5th, 2006. Our results indicate that, after a near quiescent phase from 1998 to 2002, Achernar is presently in an active phase and has built a circumstellar disk. We detect variations both in the polarization level and position angle in timescales as short as one hour and as long as several weeks. Detailed modeling of the observed polarization strongly suggests that the short-term variations originate from discrete mass ejection events which produce transient inhomogeneities in the inner disk. Long-term variations, on the other hand, can be explained by the formation of an inner ring following one or several mass ejection events.Comment: 16 pages, 5 figures, Accepted to Ap

Analysis and design of a capsule landing system and surface vehicle control system for Mars exploration

Problems related to the design and control of a mobile planetary vehicle to implement a systematic plan for the exploration of Mars are reported. Problem areas include: vehicle configuration, control, dynamics, systems and propulsion; systems analysis, terrain modeling and path selection; and chemical analysis of specimens. These tasks are summarized: vehicle model design, mathematical model of vehicle dynamics, experimental vehicle dynamics, obstacle negotiation, electrochemical controls, remote control, collapsibility and deployment, construction of a wheel tester, wheel analysis, payload design, system design optimization, effect of design assumptions, accessory optimal design, on-board computer subsystem, laser range measurement, discrete obstacle detection, obstacle detection systems, terrain modeling, path selection system simulation and evaluation, gas chromatograph/mass spectrometer system concepts, and chromatograph model evaluation and improvement

A discrete chemo-dynamical model of the dwarf spheroidal galaxy Sculptor: mass profile, velocity anisotropy and internal rotation

We present a new discrete chemo-dynamical axisymmetric modeling technique, which we apply to the dwarf spheroidal galaxy Sculptor. The major improvement over previous Jeans models is that realistic chemical distributions are included directly in the dynamical modelling of the discrete data. This avoids loss of information due to spatial binning and eliminates the need for hard cuts to remove contaminants and to separate stars based on their chemical properties. Using a combined likelihood in position, metallicity and kinematics, we find that our models naturally separate Sculptor stars into a metal-rich and a metal-poor population. Allowing for non-spherical symmetry, our approach provides a central slope of the dark matter density of $\gamma = 0.5 \pm 0.3$. The metal-rich population is nearly isotropic (with $\beta_r^{red} = 0.0\pm0.1$) while the metal-poor population is tangentially anisotropic (with $\beta_r^{blue} = -0.2\pm0.1$) around the half light radius of $0.26$ kpc. A weak internal rotation of the metal-rich population is revealed with $v_{max}/\sigma_0 = 0.15 \pm 0.15$. We run tests using mock data to show that a discrete dataset with $\sim 6000$ stars is required to distinguish between a core ($\gamma = 0$) and cusp ($\gamma = 1$), and to constrain the possible internal rotation to better than $1\,\sigma$ confidence with our model. We conclude that our discrete chemo-dynamical modelling technique provides a flexible and powerful tool to robustly constrain the internal dynamics of multiple populations, and the total mass distribution in a stellar system.Comment: Accepted by MNRA

The role of fragmentation in interaction of meteoroids with the Earth's atmosphere

As a rule, when analyzing the mechanism of quasi continuous fragmentation (QCF) it is assumed that fragments separated from the parent meteoroid (PM) are of equal initial mass. In reality, this may not be so. A major difficulty is the lack of observational data on the function of the fragment initial mass distribution and so one must resort to theoretical modeling. A discrete distribution is considered which excluded to a certain extent some mathematical difficulties. The calculation results are given and discussed

Detailed modeling of hydrodynamics mass transfer and chemical reactions in a bubble column using a discrete bubble model

A 3D discrete bubble model is adopted to investigate complex behavior involving hydrodynamics, mass transfer and chemical reactions in a gas¿liquid bubble column reactor. In this model a continuum description is adopted for the liquid phase and additionally each individual bubble is tracked in a Lagrangian framework, while accounting for bubble¿bubble and bubble¿wall interactions via an encounter model. The mass transfer rate is calculated for each individual bubble using a surface renewal model accounting for the instantaneous and local properties of the liquid phase in its vicinity. The distributions in space of chemical species residing in the liquid phase are computed from the coupled species balances considering the mass transfer from bubbles and reactions between the species. The model has been applied to simulate chemisorption of CO2 bubbles in NaOH solutions. Our results show that apart from hydrodynamics behavior, the model is able to predict the bubble size distribution as well as temporal and spatial variations of each chemical species involved