Relativistic Singular Isothermal Toroids

We construct self-similar, axisymmetric, time-independent solutions to Einstein's field equations for an isothermal gas with a flat rotation curve in the equatorial plane. The metric scales as $ds^2 \to a^2 ds^2$ under the transformation $r\to a r$ and $t \to a^{1-n} t$, where $n$ is a dimensionless measure of the strength of the gravitational field. The solution space forms a two-parameter family characterized by the ratios of the isothermal sound speed and the equatorial rotation speed to the speed of light. The isodensity surfaces are toroids, empty of matter along the rotation axis. Unlike the Newtonian case, the velocity field is not constant on a cylindrical radius. As the configuration rotates faster, an ergoregion develops in the form of the exterior of a cone centered about the rotation axis. The sequence of solutions terminates when frame dragging becomes infinite and the ergocone closes onto the axis. The fluid velocity of the last solution has finite value in the midplane but reaches the speed of light on the axis.Comment: 11 pages, 8 figure

Computer program documentation for a subcritical wing design code using higher order far-field drag minimization

A subsonic, linearized aerodynamic theory, wing design program for one or two planforms was developed which uses a vortex lattice near field model and a higher order panel method in the far field. The theoretical development of the wake model and its implementation in the vortex lattice design code are summarized and sample results are given. Detailed program usage instructions, sample input and output data, and a program listing are presented in the Appendixes. The far field wake model assumes a wake vortex sheet whose strength varies piecewise linearly in the spanwise direction. From this model analytical expressions for lift coefficient, induced drag coefficient, pitching moment coefficient, and bending moment coefficient were developed. From these relationships a direct optimization scheme is used to determine the optimum wake vorticity distribution for minimum induced drag, subject to constraints on lift, and pitching or bending moment. Integration spanwise yields the bound circulation, which is interpolated in the near field vortex lattice to obtain the design camber surface(s)

Stability of Magnetized Disks and Implications for Planet Formation

This paper considers gravitational perturbations in geometrically thin disks with rotation curves dominated by a central object, but with substantial contributions from magnetic pressure and tension. The treatment is general, but the application is to the circumstellar disks that arise during the gravitational collapse phase of star formation. We find the dispersion relation for spiral density waves in these generalized disks and derive the stability criterion for axisymmetric $(m=0)$ disturbances (the analog of the Toomre parameter $Q_T$) for any radial distribution of the mass-to-flux ratio $\lambda$. The magnetic effects work in two opposing directions: on one hand, magnetic tension and pressure stabilize the disk against gravitational collapse and fragmentation; on the other hand, they also lower the rotation rate making the disk more unstable. For disks around young stars the first effect generally dominates, so that magnetic fields allow disks to be stable for higher surface densities and larger total masses. These results indicate that magnetic fields act to suppress the formation of giant planets through gravitational instability. Finally, even if gravitational instability can form a secondary body, it must lose an enormous amount of magnetic flux in order to become a planet; this latter requirement represents an additional constraint for planet formation via gravitational instability and places a lower limit on the electrical resistivity.Comment: accepted in Ap

Bose-Einstein condensation in linear sigma model at Hartree and large N approximation

The BEC of charged pions is investigated in the framework of O(4) linear sigma model. By using Cornwall-Jackiw-Tomboulis formalism, we have derived the gap equations for the effective masses of the mesons at finite temperature and finite isospin density. The BEC is discussed in chiral limit and non-chiral limit at Hartree approximation and also at large N approximation.Comment: 11 pages, 9 figure

Computing A Glimpse of Randomness

A Chaitin Omega number is the halting probability of a universal Chaitin (self-delimiting Turing) machine. Every Omega number is both computably enumerable (the limit of a computable, increasing, converging sequence of rationals) and random (its binary expansion is an algorithmic random sequence). In particular, every Omega number is strongly non-computable. The aim of this paper is to describe a procedure, which combines Java programming and mathematical proofs, for computing the exact values of the first 64 bits of a Chaitin Omega: 0000001000000100000110001000011010001111110010111011101000010000. Full description of programs and proofs will be given elsewhere.Comment: 16 pages; Experimental Mathematics (accepted