Using recent compilations of detailed X-ray observations and spectral models of exceptional quality, we record the electron cyclotron resonance absorption (ECRA) features that have been detected in 45 pulsating high-mass X-ray binaries (HMXBs) and ultraluminous X-ray (ULX) sources harboring neutron stars, although seven of these detections are still questionable and another 21 are single and/or not independently confirmed. From the comprehensive catalogs of Jaisawal \& Naik and Staubert et al. and from several additional recent observations, we produce two lists of HMXB ECRA sources: a list of 17 sources in which multiple ECRA lines or single very low-energy lines are seen, in which we can reasonably assume that the lowest energy reveals the fundamental cyclotron level for each source; and a `contaminated' list of 38 sources including the 21 detections of single ECRA lines that may (not) be higher-level harmonics. Both lists confirm a previous result that we have obtained independently by modeling the propeller lines of Magellanic HMXB pulsars: the surface dipolar magnetic fields $B_*$ of HMXB neutron stars are segregated around five distinct values with $B_* = 0.28\pm 0.08, 0.55\pm 0.11, 1.3\pm 0.37, 3.0\pm 0.68$, and $7.9\pm 3.1$, in units of TG. An explanation of this phenomenon is currently lacking. We have found no correlation between these $B_*$ values and the corresponding observed spin periods, spin period derivatives, orbital periods, maximum X-ray luminosities, neutron star masses, or companion star masses.

We fit an isothermal oscillatory density model to the disk of HD 163296 in which planets have presumably already formed and they are orbiting at least within the four observed dark gaps. This 156 AU large axisymmetric disk shows various physical properties comparable to those of AS 209, HL Tau, and TW Hya that we have modeled previously; but it compares best to AS 209. The disks of HD 163296 and AS 209 are comparable in size and they share similar values of the power-law index $k\approx 0$ (a radial density profile of the form $\rho(R)\propto R^{-1}$), the rotational parameter $\beta_0$ (to within a factor of 3); a relatively small inner core radius (although this parameter for HD 163296 is exceptionally small, $R_1\simeq 0.15$ AU, presumably due to unresolved planets in the inner 50 AU); the scale length $R_0$ and the Jeans gravitational frequency $\Omega_J$ (to within factors of 1.4); the equation of state ($c_0^2/\rho_0$) and the central density $\rho_0$ (to within factors of 2); and the core angular velocity $\Omega_0$ (to within a factor of 4.5). In the end, we compare all six nebular disks that we have modeled so far.

We fit an isothermal oscillatory density model of Saturn's protoplanetary disk to the present-day major satellites and innermost rings D/C and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Saturnian nebula. This disk does not look like the Jovian and Uranian disks that we modeled previously. Its power-law index is extremely steep ($k=-4.5$) and its radial extent is very narrow ($\Delta R\lesssim 0.9$ Gm), its rotation parameter that measures centrifugal support against self-gravity is somewhat larger ($\beta_0=0.0431$), as is its radial scale length (395 km); but, as was expected, the size of the Saturnian disk, $R_{\rm max}=3.6$ Gm, takes just an intermediate value. On the other hand, the central density of the compact Saturnian core and its angular velocity are both comparable to that of Jupiter's core (density of $\approx 0.3$~g~cm$^{-3}$ in both cases, and rotation period of 5.0 d versus 6.8 d); and significantly less than the corresponding parameters of Uranus' core. As with the other primordial nebulae, this rotation is sufficiently slow to guarantee the disk's long-term stability against self-gravity induced instabilities for millions of years of evolution.

We fit an isothermal oscillatory density model to two ALMA/DSHARP-observed disks, AS 209 and HL Tau, in which planets have presumably already formed and they are orbiting within the observed seven dark gaps in each system. These large disks are roughly similar to our solar nebula, albeit they exhibit milder radial density profiles and they enjoy lower centrifugal support. We find power-law density profiles with index $k=0.0$ (radial densities $\rho(R) \propto R^{-1}$) and centrifugal support against self-gravity so small that it guarantees dynamical stability for millions of years of evolution. The scale lengths of the models differ only by a factor of 1.9, but the inner cores of the disks are very different: HL Tau's core is 8.0 times larger and 3.6 times denser than the core of AS 209. This results in four dark gaps having formed within the core of HL Tau, whereas no dark gap is found in the core of AS 209. On the other hand, the Jeans frequencies and the angular velocities of the cores are comparable to within factors of 1.9 and 1.6, respectively.

Recent ALMA/DSHARP observations reveal that the disk of RU Lup shows presently only four dark gaps at large radii in which planets are presumably already orbiting. This gap arrangement is quite similar to the current structure of the outer planets in our solar system. Although there may be more undetected planets forming in the inner disk of RU Lup, the information obtained from the current observations prompted us to fit preliminary models with an isothermal oscillatory density profile to the outer disk of RU Lup. The best-fit model confirms that several physical properties of the disk of RU Lup are comparable to those of our solar nebula to within factors of 1.0-2.3; but there are some differences as well: central density ($\rho_0$), centrifugal support ($\beta_0$), and core angular velocity ($\Omega_0$) in RU Lup are smaller by factors of 5.5, 8.0, and 18, respectively, whereas the inner core radius ($R_1$) of RU Lup is 25 times larger. The smaller values found for RU Lup are linked to the (apparently gap-empty) inner 10 AU of its disk, a region that is presently inadequately resolved.

We fit an isothermal oscillatory density model to the outer disk of TW Hya in which planets have presumably already formed and they are orbiting within four observed dark gaps. At first sight, this 52 AU small disk does not appear to be similar to our solar nebula; it shows several physical properties comparable to those in HL Tau (size $R_{\rm max}=102$ AU) and very few similarities to AS 209 ($R_{\rm max}=144$ AU). We find a power-law density profile with index $k=-0.2$ (radial densities $\rho(R) \propto R^{-1.2}$) and centrifugal support against self-gravity so small that it virtually guarantees dynamical stability for millions of years of evolution to come. Compared to HL Tau, the scale length $R_0$ and the core size $R_1$ of TW Hya are smaller only by factors of $\sim$2, reflecting the disk's half size. On the opposite end, the Jeans frequency $\Omega_J$ and the angular velocity $\Omega_0$ of the smaller core of TW Hya are larger only by factors of $\sim$2. The only striking difference is that the central density ($\rho_0$) of TW Hya is 5.7 times larger than that of HL Tau, which is understood because the core of TW Hya is only half the size ($R_1$) of HL Tau and about twice as heavy ($\Omega_J$). In the end, we compare the protostellar disks that we have modeled so far.

In view of the many recent observations conducted by ALMA and SPHERE, it is becoming clear that protoplanetary disks form planets in narrow annular gaps at various distances from the central protostars before these protostars are actually fully formed and the gaseous disks have concluded their accretion/dispersal processes. This is in marked contrast to the many multi-planet exoplanetary systems that do not conform to this pristine picture. This major discrepancy calls for an explanation. We provide such an explanation in this work, based on analytical solutions of the cylindrical isothermal Lane-Emden equation with rotation which do not depend on boundary conditions. These ``intrinsic'' solutions of the differential equation attract the solutions of the Cauchy problem and force them to oscillate permanently. The oscillations create density maxima in which dust and planetesimals are trapped and they can form protoplanetary cores during the very early isothermal evolution of such protoplanetary nebulae. We apply this model to our solar nebula that formed in-situ a minimum of eleven protoplanetary cores that have grown to planets which have survived undisturbed to the present day. We are also in the process of applying the same model to the ALMA/DSHARP disks.

We have proposed an alternative model for the formation of our solar system that does not predict any mean-motion resonant interactions, planetary migrations, or self-gravitating instabilities in the very early isothermal solar nebula and before the protosun has formed. Within this context of nonviolent protoplanetary evolution over more than 10 million years, we examine some conundrums and constraints that have been discovered from studies of small bodies in the present-day solar system (Jupiter and Neptune's Trojans and their differences from Kuiper belt objects, the irregular satellites of gaseous giants, the stability of the main asteroid belt, and the Late Heavy Bombardment). These issues that have caused substantial difficulties to models of violent formation do not appear to be problematic for the alternative model, and the reason is the complete lack of violent events during the evolution of protoplanets.

We fit an isothermal oscillatory density model of Jupiter's protoplanetary disk to the present-day Galilean and other nearby satellites and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Jovian nebula. Although the radial density profile of Jupiter's disk was similar to that of the solar nebula, its rotational support against self-gravity was very low, a property that also guaranteed its long-term stability against self-gravity induced instabilities for millions of years.

We fit an isothermal oscillatory density model of Uranus' protoplanetary disk to the present-day major satellites and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Uranian nebula. This disk does not at all look like the Jovian disk that we modeled previously. Its rotation parameter that measures centrifugal support against self-gravity is a lot smaller ($\beta_0=0.00507$), as is the radial scale length (only 27.6 km) and the size of the disk (only 0.60 Gm). On the other hand, the central density of the compact Uranian core is higher by a factor of 180 and its core's angular velocity is about 2.3 times that of Jupiter's core (a rotation period of 3.0 d as opposed to 6.8 d). Yet, the rotation of the disk is sufficiently slow to guarantee its long-term stability against self-gravity induced instabilities for millions of years.

We fit an isothermal oscillatory density model of Neptune's protoplanetary disk to the surviving regular satellites and its innermost ring and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Neptunian nebula. Neptune's regular moons suffered from the retrograde capture of Triton that disrupted the system. Some moons may have been ejected, while others may have survived inside their potential minima. For this reason, the Neptunian nebula does not look like any of the nebulae that we modeled previously. In particular, there must be two density maxima deep inside the core of the nebula where no moons or rings are found nowadays. Even with this strong assumption, the recent discovery of the minor moon N XIV complicates further the modeling effort. With some additional assumptions, the Neptunian nebula still shares many similarities with the Uranian nebula, as was expected from the relative proximity and similar physical conditions of the two systems. For Neptune's primordial disk, we find a steep power-law index ($k=-3.0$), needed to accommodate the arrangement of the outer moons Larissa, N XIV, and Proteus. The rotation parameter that measures centrifugal support against self-gravity is quite small ($\beta_0=0.00808$), as is its radial scale length (13.6 km). The extent of the disk ($R_{\rm max}=0.12$ Gm) is a lot smaller than that of Uranus ($R_{\rm max}=0.60$ Gm) and Triton appears to be responsible for the truncation of the disk. The central density of the compact Neptunian core and its angular velocity are higher than but comparable to those of Uranus' core. In the end, we compare the models of the protoplanetary disks of the four gaseous giants.

We have recently shown that the baryonic Tully-Fisher and Faber-Jackson relations imply that the gravitational "constant" $G$ in the force law varies with acceleration $a$ as $G\propto 1/a$ and vice versa. These results prompt us to reconsider every facet of Newtonian dynamics. Here we show that the integral form of Gauss's law in spherical symmetry remains valid in $G(a)$ gravity, but the differential form depends on the precise distribution of $G(a)M(r)$, where $r$ is the distance from the origin and $M(r)$ is the mass distribution. We derive the differential form of Gauss's law in spherical symmetry, thus the source for Poisson's equation as well. Modified Newtonian dynamics (MOND) and weak-field Weyl gravity are asymptotic limits of $G(a)$ gravity at low and high accelerations, respectively. In these limits, we derive telling approximations to the source in spherical symmetry. It turns out that the source has a strong dependence on surface density $M/r^2$ everywhere in $a$-space except in the deep Newton-Weyl regime of very high accelerations.

We have recently shown that the baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations imply that the gravitational "constant" $G$ in the force law vary with acceleration $a$ as $1/a$. Here we derive the converse from first principles. First we obtain the gravitational potential for all accelerations and we formulate the Lagrangian for the central-force problem. Then action minimization implies the BTF/BFJ relations in the deep MOND limit as well as weak-field Weyl gravity in the Newtonian limit. The results show how we can properly formulate a nonrelativistic conformal theory of modified dynamics that reduces to MOND in its low-acceleration limit and to Weyl gravity in the opposite limit. An unavoidable conclusion is that $a_0$, the transitional acceleration in modified dynamics, does not have a cosmological origin and it may not even be constant among galaxies and galaxy clusters.

The Newtonian gravitational constant $G$ obeys the dimensional relation $[G] [M] [a] = [v]^4$, where $M$, $a$, and $v$ denote mass, acceleration, and speed, respectively. Since the baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations are observed facts, this relation implies that $G\, a = {\rm constant}$. This result cannot be obtained in Newtonian dynamics which cannot explain the origin of the BTF and BFJ relations. An alternative, modified Newtonian dynamics (MOND) assumes that $G=G_0$ is constant in space and derives naturally a characteristic constant acceleration $a=a_0$, as well as the BTF and BFJ relations. This is overkill and it comes with a penalty: MOND cannot explain the origin of $a_0$. A solid physical resolution of this issue is that $G \propto a^{-1}$, which implies that in lower-acceleration environments the gravitational force is boosted relative to its Newtonian value because $G$ increases. This eliminates all problems related to MOND's empirical cutoff $a_0$ and yields a quantitative method for mapping the detailed variations of $G(a)$ across each individual galaxy as well as on larger and smaller scales. On the opposite end, the large accelerations produced by $G(a)$ appear to be linked to the weak-field limit of the fourth-order theory of conformal Weyl gravity.

Three high-mass X-ray binaries have been discovered recently exhibiting enormous spinup rates. Conventional accretion theory predicts extremely high surface dipolar magnetic fields that we believe are unphysical. Instead, we propose quite the opposite scenario: some of these pulsars exhibit weak magnetic fields, so much so that their magnetospheres are crushed by the weight of inflowing matter. The enormous spinup rate is achieved before inflowing matter reaches the pulsar's surface as the penetrating inner disk transfers its excess angular momentum to the receding magnetosphere which, in turn, applies a powerful spinup torque to the pulsar. This mechanism also works in reverse: it spins a pulsar down when the magnetosphere expands beyond corotation and finds itself rotating faster than the accretion disk which then exerts a powerful retarding torque to the magnetic field and to the pulsar itself. The above scenaria cannot be accommodated within the context of neutron-star accretion processes occurring near spin equilibrium, thus they constitute a step toward a new theory of extreme (far from equilibrium) accretion phenomena.

We analyse the most powerful X-ray outbursts from neutron stars in ten Magellanic high-mass X-ray binaries and three pulsating ultraluminous X-ray sources. Most of the outbursts rise to $L_{max}$ which is about the level of the Eddington luminosity, while the rest and more powerful outbursts also appear to recognize that limit when their emissions are assumed to be anisotropic and beamed toward our direction. We use the measurements of pulsar spin periods $P_S$ and their derivatives $\dot{P_S}$ to calculate the X-ray luminosities $L_p$ in their faintest accreting ("propeller") states. In four cases with unknown $\dot{P_S}$, we use the lowest observed X-ray luminosities, which only adds to the heterogeneity of the sample. Then we calculate the ratios $L_p/L_{max}$ and we obtain an outstanding confluence of theory and observations from which we conclude that work done on both fronts is accurate and the results are trustworthy: sources known to reside on the lowest Magellanic propeller line are all located on/near that line, whereas other sources jump higher and reach higher-lying propeller lines. These jumps can be interpreted in only one way, higher-lying pulsars have stronger surface magnetic fields in agreement with empirical results in which $\dot{P_S}$ and $L_p$ values were not used.

In spherical symmetry with radial coordinate $r$, classical Newtonian gravitation supports circular orbits and, for $-1/r$ and $r^2$ potentials only, closed elliptical orbits [1]. Various families of elliptical orbits can be thought of as arising from the action of perturbations on corresponding circular orbits. We show that one elliptical orbit in each family is singled out because its focal length is equal to the radius of the corresponding unperturbed circular orbit. The eccentricity of this special orbit is related to the famous irrational number known as the golden ratio. So inanimate Newtonian gravitation appears to exhibit (but not prefer) the golden ratio which has been previously identified mostly in settings within the animate world.

We present new families of bound, closed, nonelliptical orbits that are supported by various spherical potentials in clear contradiction to Newton's and Bertrand's theorems. We calculate analytically some typical closed orbits of arbitrarily large amplitudes for the spherically symmetric forms of the Newton-Kepler, Hooke, and linear potentials. In all cases, the radial oscillations of the orbits exhibit unequal amplitudes and nonsynchronous frequencies on either side of a circular equilibrium orbit. These parameters are however related through energy conservation and orbit continuity across the equilibrium orbit. The constraints placed on the parameters are all localized to each particular orbit and depend crucially on the chosen initial conditions.

We consider the geometric Titius-Bode rule for the semimajor axes of planetary orbits. We derive an equivalent rule for the midpoints of the segments between consecutive orbits along the radial direction and we interpret it physically in terms of the work done in the gravitational field of the Sun by particles whose orbits are perturbed around each planetary orbit. On such energetic grounds, it is not surprising that some exoplanets in multiple-planet extrasolar systems obey the same relation. But it is surprising that this simple interpretation of the Titius-Bode rule also reveals new properties of the bound closed orbits predicted by Bertrand's theorem and known since 1873.

We have compiled a comprehensive library of all X-ray observations of Magellanic pulsars carried out by {\it XMM-Newton}, {\it Chandra}, and {\it RXTE} in the period 1997-2014. In this work, we use the data from 53 high-mass Be/X-ray binaries in the Small Magellanic Cloud to demonstrate that the distribution of spin-period derivatives vs. spin periods of spinning-down pulsars is not at all different than that of the accreting spinning-up pulsars. The inescapable conclusion is that the up and down samples were drawn from the same continuous parent population, therefore Be/X-ray pulsars that are spinning down over periods spanning 18 years are in fact accreting from retrograde disks. The presence of prograde and retrograde disks in roughly equal numbers supports a new evolutionary scenario for Be/X-ray pulsars in their spin period-period derivative diagram.

The recent discoveries of pulsed X-ray emission from three ultraluminous X-ray (ULX) sources have finally enabled us to recognize a subclass within the ULX class: the great pretenders, neutron stars (NSs) that appear to emit X-ray radiation at isotropic luminosities $L_X = 7\times 10^{39}$~erg~s$^{-1}-1\times 10^{41}$~erg~s$^{-1}$ only because their emissions are strongly beamed toward our direction and our sight lines are offset by only a few degrees from their magnetic-dipole axes. The three known pretenders appear to be stronger emitters than the presumed black holes of the ULX class, such as Holmberg II \& IX X-1, IC10 X-1, and NGC300 X-1. For these three NSs, we have adopted a single reasonable assumption, that their brightest observed outbursts unfold at the Eddington rate, and we have calculated both their propeller states and their surface magnetic-field magnitudes. We find that the results are not at all different from those recently obtained for the Magellanic Be/X-ray pulsars: the three NSs reveal modest magnetic fields of about 0.3-0.4~TG and beamed propeller-line X-ray luminosities of $\sim 10^{36-37}$~erg~s$^{-1}$, substantially below the Eddington limit.

We have collected and analyzed the complete archive of {\itshape XMM-Newton\} (116), {\itshape Chandra\} (151), and {\itshape RXTE\} (952) observations of the Small Magellanic Cloud (SMC), spanning 1997-2014. The resulting observational library provides a comprehensive view of the physical, temporal and statistical properties of the SMC pulsar population across the luminosity range of $L_X= 10^{31.2}$--$10^{38}$~erg~s$^{-1}$. From a sample of 67 pulsars we report $\sim$1654 individual pulsar detections, yielding $\sim$1260 pulse period measurements. Our pipeline generates a suite of products for each pulsar detection: spin period, flux, event list, high time-resolution light-curve, pulse-profile, periodogram, and spectrum. Combining all three satellites, we generated complete histories of the spin periods, pulse amplitudes, pulsed fractions and X-ray luminosities. Some pulsars show variations in pulse period due to the combination of orbital motion and accretion torques. Long-term spin-up/down trends are seen in 12/11 pulsars respectively, pointing to sustained transfer of mass and angular momentum to the neutron star on decadal timescales. Of the sample 30 pulsars have relatively very small spin period derivative and may be close to equilibrium spin. The distributions of pulse-detection and flux as functions of spin-period provide interesting findings: mapping boundaries of accretion-driven X-ray luminosity, and showing that fast pulsars ($P<$10 s) are rarely detected, which yet are more prone to giant outbursts. Accompanying this paper is an initial public release of the library so that it can be used by other researchers. We intend the library to be useful in driving improved models of neutron star magnetospheres and accretion physics.

The 80 high-mass X-ray binary (HMXB) pulsars that are known to reside in the Magellanic Clouds (MCs) have been observed by the XMM-Newton and Chandra X-ray telescopes on a regular basis for 15 years, and the XMM-Newton and Chandra archives contain nearly complete information about the duty cycles of the sources with spin periods P_S < 100 s. We have rerprocessed the archival data from both observatories and we combined the output products with all the published observations of 31 MC pulsars with P_S < 100 s in an attempt to investigate the faintest X-ray emission states of these objects that occur when accretion to the polar caps proceeds at the smallest possible rates. These states determine the so-called propeller lines of the accreting pulsars and yield information about the magnitudes of their surface magnetic fields. We have found that the faintest states of the pulsars segregate into five discrete groups which obey to a high degree of accuracy the theoretical relation between spin period and X-ray luminosity. So the entire population of these pulsars can be described by just five propeller lines and the five corresponding magnetic moments (0.29, 0.53, 1.2, 2.9, and 7.3, in units of 10^30 G cm^3).

We have derived exact axisymmetric solutions of the two-dimensional Lane-Emden equations with rotation. These solutions are intrinsically favored by the differential equations regardless of any adopted boundary conditions and the physical solutions of the Cauchy problem are bound to oscillate about and remain close to these intrinsic solutions. The isothermal solutions are described by power-law density profiles in the radial direction, whereas the polytropic solutions are described by radial density profiles that are powers of the zeroth-order Bessel function of the first kind. Both families of solutions decay exponentially in the vertical direction and both result in increasing or nearly flat radial rotation curves. The results are applicable to gaseous spiral-galaxy disks that exhibit flat rotation curves and to the early stages of protoplanetary disk formation before the central star is formed.

In young starburst galaxies, the X-ray population is expected to be dominated by the relics of the most massive and short-lived stars, black-hole and neutron-star high mass X-ray binaries (XRBs). In the closest such galaxy, IC 10, we have made a multi-wavelength census of these objects. Employing a novel statistical correlation technique, we have matched our list of 110 X-ray point sources, derived from a decade of Chandra observations, against published photometric data. We report an 8 sigma correlation between the celestial coordinates of the two catalogs, with 42 X-ray sources having an optical counterpart. Applying an optical color-magnitude selection to isolate blue supergiant (SG) stars in IC 10, we find 16 matches. Both cases show a statistically significant overabundance versus the expectation value for chance alignments. The blue objects also exhibit systematically higher fx/fv ratios than other stars in the same magnitude range. Blue SG-XRBs include a major class of progenitors of double-degenerate binaries, hence their numbers are an important factor in modeling the rate of gravitational wave sources. We suggest that the anomalous features of the IC 10 stellar population are explained if the age of the IC 10 starburst is close to the time of the peak of interaction for massive binaries.

We have monitored the Cassiopeia dwarf galaxy (IC 10) in a series of 10 Chandra ACIS-S observations to capture its variable and transient X-ray source population, which is expected to be dominated by High Mass X-ray Binaries (HMXBs). We present a sample of 21 X-ray sources that are variable between observations at the 3 sigma level, from a catalog of 110 unique point sources. We find 4 transients (flux variability ratio greater than 10) and a further 8 objects with ratio > 5. The observations span years 2003 - 2010 and reach a limiting luminosity of >10$^{35}$ erg/s, providing sensitivity to X-ray binaries in IC 10 as well as flare stars in the foreground Milky Way. The nature of the variable sources is investigated from light-curves, X-ray spectra, energy quantiles, and optical counterparts. The purpose of this study is to discover the composition of the X-ray binary population in a young starburst environment. IC 10 provides a sharp contrast in stellar population age (<10 My) when compared to the Magellanic Clouds (40-200 My) where most of the known HMXBs reside. We find 10 strong HMXB candidates, 2 probable background Active Galactic Nuclei, 4 foreground flare-stars or active binaries, and 5 not yet classifiable sources. Complete classification of the sample requires optical spectroscopy for radial velocity analysis and deeper X-ray observations to obtain higher S/N spectra and search for pulsations. A catalog and supporting dataset are provided.

{\it NuSTAR} has detected pulsations from the ultraluminous X-ray source X-2 in M82 and archival {\it Chandra} observations have given us a good idea of its duty cycle. The newly discovered pulsar exhibited at least 4 super-Eddington outbursts in the past 15 years but, in its lowest-power state, it radiates just below the Eddington limit and its properties appear to be typical of high-mass X-ray binaries. M82 X-2 has been described as a common neutron star with a 1~TG magnetic field that occasionally accretes above the Eddington rate and as a magnetar-like pulsar with a 10-100~TG magnetic field that reaches above the quantum limit. We argue in favor of the former interpretation. Using standard accretion theory and the available observations, we calculate the stellar magnetic field of this pulsar in two independent ways and we show that it cannot exceed 3~TG in either case. We discuss the implications of our results for other ultraluminous X-ray sources that occasionally exhibit similar powerful outbursts.

The spins of stellar-mass black holes (BHs) and the power outputs of their jets are measurable quantities. Unfortunately, the currently employed methods do not agree and the results are controversial. Two major issues concern the measurements of BH spin and beam (jet) power. The former issue can be resolved by future observations. But the latter issue can be resolved now, if we pay attention to what is expected from theoretical considerations. The question of whether a correlation has been found between the power outputs of few objects and the spins of their BHs is moot because BH beam power does not scale with the square of the spin of the BH. We show that the theoretical BH beam power is a strongly nonlinear function of spin that cannot be approximated by a quadratic relation, as is generally stated when the influence of the magnetic field is not accounted for in the \cite{bla77} model. The BH beam power of ballistic jets should scale a lot more steeply with BH spin irrespective of the magnetic field assumed to thread the horizon and the spin range considered. This behavior may already be visible in the analyses of radio observations by \cite{nar12} and \cite{rus13}. In agreement with previous studies, we also find that the power output that originates in the inner regions of the surrounding accretion disks is higher than that from the BHs and it cannot be ignored in investigations of continuous compact jets from these systems.

Pulsations were recently detected from the ultraluminous X-ray source X-2 in M82. The newly discovered pulsar has been described as a common neutron star with a 1 TG magnetic field that accretes above the Eddington rate and as a magnetar-like pulsar with a 100 TG magnetic field that is above the quantum limit. We show here that this X-ray source is not exotic in any sense. The standard equations of accretion along field lines predict that, for the measured spin period $P_S$ and spinup rate $\dot{P_S}$, the isotropic X-ray luminosity $L_X$ must be near the Eddington limit (i.e., $L_{X}\approx 3.5\times 10^{38}$~erg~s$^{-1}$); and that the surface magnetic field $B$, that does not depend on $P_S$, must be modest (i.e., $B_*\approx 1-10$ TG). The observed higher luminosity can be explained by a moderate amount of geometric beaming that occurs in our direction. Other ultraluminous X-ray sources may also turn out to be common pulsars with similar physical characteristics, but since the emission must occur at a favorable angle to the observer, we expect that very few such pulsars will be discovered in the future.

Helical magnetic fields embedded in the jets of active galactic nuclei (AGNs) are required by the broad range of theoretical models that advocate for electromagnetic launching of the jets. In most models, the direction of the magnetic field is random, but if the axial field is generated by a Cosmic Battery generated by current in the direction of rotation in the accretion disk, there is a correlation between the directions of the spin of the AGN accretion disk and of the axial field, which leads to a specific direction for the axial electric current, azimuthal magnetic field, and the resulting observed transverse Faraday-rotation (FR) gradient across the jet, due to the systematic change in the line-of-sight magnetic field. We consider new observational evidence for the presence of a nested helical magnetic-field structure such as would be brought about by the operation of the Cosmic Battery, and make predictions about the expected behavior of transverse FR gradients observed on decaparsec and kiloparsec scales. The collected results can be understood if the dominant azimuthal field on parsec scales corresponds to an axial electric current flowing inward along the jet, whereas the (weaker) dominant azimuthal field on kiloparsec scales corresponds to a outward-flowing current in the outer sheath of the jet and/or an extended disk wind.

By solving analytically the various types of Lane-Emden equations with rotation, we have discovered two new coupled fundamental properties of rotating, self-gravitating, gaseous disks in equilibrium: Isothermal disks must, on average, exhibit strict power-law density profiles in radius $x$ on their equatorial planes of the form $A x^{k-1}$, where $A$ and $k-1$ are the integration constants, and "flat" rotation curves precisely such as those observed in spiral galaxy disks. Polytropic disks must, on average, exhibit strict density profiles of the form $\left[\ln(A x^k)\right]^n$, where $n$ is the polytropic index, and "flat" rotation curves described by square roots of upper incomplete gamma functions. By "on average," we mean that, irrespective of the chosen boundary conditions, the actual profiles must oscillate around and remain close to the strict mean profiles of the analytic singular equilibrium solutions. We call such singular solutions the "intrinsic" solutions of the differential equations because they are demanded by the second-order equations themselves with no regard to the Cauchy problem. The results are directly applicable to gaseous galaxy disks that have long been known to be isothermal and to protoplanetary disks during the extended isothermal and adiabatic phases of their evolution. In galactic gas dynamics, they have the potential to resolve the dark matter--modified gravity controversy in a sweeping manner, as they render both of these hypotheses unnecessary. In protoplanetary disk research, they provide observers with powerful new probing tool, as they predict a clear and simple connection between the radial density profiles and the rotation curves of self-gravitating disks in their very early (pre-Class 0 and perhaps the youngest Class Young Stellar Objects) phases of evolution.

The relative phasing of the X-ray eclipse ephemeris and optical radial velocity (RV) curve for the X-ray binary IC10 X-1 suggests the He[$\lambda$4686] emission-line originates in a shadowed sector of the stellar wind that avoids ionization by X-rays from the compact object. The line attains maximum blueshift when the wind is directly toward us at mid X-ray eclipse, as is also seen in Cygnus X-3. If the RV curve is unrelated to stellar motion, evidence for a massive black hole evaporates because the mass function of the binary is unknown. The reported X-ray luminosity, spectrum, slow QPO, and broad eclipses caused by absorption/scattering in the WR wind are all consistent with either a low-stellar-mass BH or a NS. For a NS, the centre of mass lies inside the WR envelope whose motion is then far below the observed 370 km/s RV amplitude, while the velocity of the compact object is as high as 600 km/s. The resulting 0.4\% doppler variation of X-ray spectral lines could be confirmed by missions in development. These arguments also apply to other putative BH binaries whose RV and eclipse curves are not yet phase-connected. Theories of BH formation and predicted rates of gravitational wave sources may need revision.

Neutron-star and black-hole X-ray binaries (XRBs) exhibit radio jets, whose properties depend on the X-ray spectral state and history of the source. In particular, black-hole XRBs emit compact, steady radio jets when they are in the so-called hard state, the jets become eruptive as the sources move toward the soft state, disappear in the soft state, and re-appear when the sources return to the hard state. On the other hand, jets from neutron-star X-ray binaries are typically weaker radio emitters than the black-hole ones at the same X-ray luminosity and in some cases radio emission is detected in the soft state. Significant phenomenology has been accumulated so far regarding the spectral states of neutron-star and black-hole XRBs, and there is general agreement about the type of the accretion disk around the compact object in the various spectral states. Our aim is to investigate whether the phenomenology regarding the X-ray emission on one hand and the jet appearance and disappearance on the other can be put together in a consistent physical picture. It has been shown that the so-called Poynting-Robertson Cosmic Battery (PRCB) explains in a natural way the formation of magnetic fields in the disks of AGN and the ejection of jets. We investigate whether the PRCB can also explain the formation, destruction, and variability of jets in XRBs. We find excellent agreement between the conditions under which the PRCB is efficient (i.e., the type of the accretion disk) and the emission or destruction of the radio jet. The disk-jet connection in XRBs is explained in a natural way using the PRCB.

The origin of cosmic magnetic (B) fields remains an open question. It is generally believed that very weak primordial B fields are amplified by dynamo processes, but it appears unlikely that the amplification proceeds fast enough to account for the fields presently observed in galaxies and galaxy clusters. In an alternative scenario, cosmic B fields are generated near the inner edges of accretion disks in Active Galactic Nuclei (AGNs) by azimuthal electric currents due to the difference between the plasma electron and ion velocities that arises when the electrons are retarded by interactions with photons. While dynamo processes show no preference for the polarity of the (presumably random) seed field that they amplify, this alternative mechanism uniquely relates the polarity of the poloidal B field to the angular velocity of the accretion disk, resulting in a unique direction for the toroidal B field induced by disk rotation. Observations of the toroidal fields of 29 AGN jets revealed by parsec-scale Faraday rotation measurements show a clear asymmetry that is consistent with this model, with the probability that this asymmetry came about by chance being less than 1%. This lends support to the hypothesis that the Universe is seeded by B fields that are generated in AGN via this mechanism and subsequently injected into intergalactic space by the jet outflows.

We have derived exact solutions of the isothermal Lane--Emden equation with and without rotation in a cylindrical geometry. The corresponding hydrostatic equilibria are most relevant to the dynamics of the protosolar nebula before and during the stages of planet and satellite formation. The nonrotating solution for the mass density is analytic, nonsingular, monotonically decreasing with radius, and it satisfies easily the usual physical boundary conditions at the center. When differential rotation is added to the Lane--Emden equation, an entire class of exact solutions for the mass density appears. We have determined all of these solutions analytically as well. Within this class, solutions that are power laws or combinations of power laws are not capable of satisfying the associated boundary--value problem, but they are nonetheless of profound importance because they constitute "baselines" to which the actual solutions approach when the central boundary conditions are imposed. Numerical integrations that enforce such physical boundary conditions show that the actual radial equilibrium density profiles emerge from the center close to the nonrotating solution, but once they cross below the corresponding baselines, they cease to be monotonic. The actual solutions are forced to oscillate permanently about the baseline solutions without ever settling onto them because the central boundary conditions strictly prohibit such settling, even in the asymptotic regime of large radii. Based on our results, we expect that quasistatically--evolving protoplanetary disks should develop oscillatory density profiles in their midplanes during their isothermal phase. The peaks in these profiles correspond to local potential minima and their locations are ideal sites for the formation of protoplanets ...

Five planets are presently believed to orbit the primary star of 55 Cnc, but there exists a large 5 AU gap in their distribution between the two outermost planets. This gap has attracted considerable interest because it may contain one or more lower--mass planets whose existence is not contradicted by long-term orbit stability analyses, in fact it is expected according to the "packed planetary systems" hypothesis and an empirical Titius--Bode relation recently proposed for 55 Cnc. Furthermore, the second largest planet is just the second farthest and it orbits very close to the star. Its orbit, the most circular of all, appears to be nearly but not quite commensurable with the orbit of the third planet, casting doubt that any migration or resonant capture of the inner planets has ever occurred and lending support to the idea of "in--situ" giant planet formation by the process of core accretion. All of the above ideas will be tested in the coming years in this natural laboratory as more observations will become available. This opportunity presents itself in conjunction with a physical model that relates the orbits of the observed planets to the structure of the original protoplanetary disk that harbored their formation at the early stages of protostellar collapse. Using only the 5 observed planets of 55 Cnc, this model predicts that the surface density profile of its protoplanetary disk varied with distance $R$ precisely as $\Sigma (R)\propto R^{-3/2}$, as was also found for the minimum--mass solar nebula. Despite this similarity, the disk of 55 Cnc was smaller, heavier, and less rotationally supported than the solar nebula, so this system represents a different mode of multi--planet formation compared to our own solar system.

This paper has been removed by arXiv admin because it was an erroneous duplicate of astro-ph/9411031.

We describe the results of numerical "2.5--dimensional" MHD simulations of an initially unmagnetized disk model orbiting a central point--mass and responding to the continual generation of poloidal magnetic field due to a secular source that emulates the Poynting--Robertson (PR) drag on electrons in the vicinity of a luminous stellar or compact accreting object. The fluid in the disk and in the surrounding hotter atmosphere has finite electrical conductivity and allows for the magnetic field to diffuse freely out of the areas where it is generated, while at the same time, the differential rotation of the disk twists the poloidal field and quickly induces a substantial toroidal--field component. The secular PR term has dual purpose in these simulations as the source of the magnetic field and the trigger of a magnetorotational instability (MRI) in the disk. The MRI is especially mild and does not destroy the disk because a small amount of resistivity dampens the instability efficiently. In simulations with moderate resistivities (diffusion timescales up to $\sim$16 local dynamical times) and after $\sim$100 orbits, the MRI has managed to transfer outward substantial amounts of angular momentum and the inner edge of the disk, along with azimuthal magnetic flux, has flowed toward the central point--mass where a new, magnetized, nuclear disk has formed. The toroidal field in this nuclear disk is amplified by differential rotation and it cannot be contained; when it approaches equipartition, it unwinds vertically and produces episodic jet--like outflows. The poloidal field in the inner region cannot diffuse back out if the characteristic diffusion time is of the order of or larger than the dynamical time; it continues to grow linearly in time undisturbed and without saturation, as the outer sections of many poloidal loops are being drawn radially outward.

We obtain the general forms of the axisymmetric stability criteria in a magnetized compressible Couette flow using an energy variational principle, the so--called interchange or Chandrasekhar's method, which we applied successfully in the incompressible case. This formulation accounts for the simultaneous presence of gravity, rotation, a toroidal magnetic field, a weak axial magnetic field, entropy gradients, and density gradients in the initial equilibrium state. The power of the method lies in its simplicity which allows us to derive extremely compact and physically clear expressions for the relevant stability criteria despite the inclusion of so many physical effects. In the implementation of the method, all the applicable conservation laws are explicitly taken into account during the variations of a quantity with dimensions of energy which we call the "free--energy function." In the limits of plane--parallel atmospheres and homoentropic flows, our formulation easily recovers the stability criteria for suppression of convective and Parker instabilities, as well as some related special cases studied over 40 years ago by Newcomb and Tserkovnikov via laborious variational techniques.

We have searched the IRAS Point Source Catalog and HIRES maps for young stellar objects (YSOs) in the direction of five \HI high-velocity clouds (HVCs). In agreement with optical searches in the halo, no evidence was found for extensive star-forming activity inside the high-latitude HVCs. Specifically, we have found no signs of star formation or YSOs in the direction of the A IV cloud or in the very-high-velocity clouds HVC~110-7-465 and HVC~114-10-440. We have identified only one young star in the direction of the M~I.1 cloud, which shows almost perfect alignment with a knot of \HI emission. Because of the small number of early-type stars observed in the halo, the probability for such a positional coincidence is low; thus, this young star appears to be physically associated with the M~I.1 cloud. We have also identified a good YSO candidate in the \HI shell-like structure observed in the core region of the low-latitude cloud complex H (HVC~131+1-200). This region could be a supernova remnant with several other YSO candidates formed along the shock front produced by the explosion. In agreement with recent theoretical estimates, these results point to a low but significant star-formation rate in intermediate and high Galactic latitude HVCs. For M~I.1 in particular, we estimate that the efficiency of the star-formation process is $M(YSO)/M(\HI)\ga 10^{-4}-10^{-3}$ by mass. Such efficiency is sufficient to account for (a) the existence of the few young blue stars whose ages imply that they were born in the Galactic halo, and (b) the nonprimordial metallicities inferred for some HVCs if their metal content proves to be low.

We use a free-energy minimization approach to describe the secular and dynamical instabilities as well as the bifurcations along equilibrium sequences of rotating, self-gravitating fluid systems. Our approach is fully nonlinear and stems from the Landau-Ginzburg theory of phase transitions. Here we examine higher than 2nd-harmonic disturbances applied to Maclaurin spheroids, the corresponding bifurcating sequences, and their relation to nonlinear fission processes. The triangle and ammonite sequences bifurcate from the two 3rd-harmonic neutral points on the Maclaurin sequence while the square and one-ring sequences bifurcate from two of the three known 4th harmonic neutral points. In the other three cases, secular instability does not set in at the corresponding bifurcation points because the sequences stand and terminate at higher energies relative to the Maclaurin sequence. There is no known bifurcating sequence at the point of 3rd-harmonic dynamical instability. Our nonlinear approach easily identifies resonances between the Maclaurin sequence and various multi-fluid-body sequences that cannot be detected by linear stability analyses. Resonances appear as first-order phase transitions at points where the energies of the two sequences are nearly equal but the lower energy state belongs to one of the multi-fluid-body sequences.

In Papers I and II, we have used a free-energy minimization approach that stems from the Landau-Ginzburg theory of phase transitions to describe in simple and clear physical terms the secular and dynamical instabilities as well as the bifurcations along equilibrium sequences of rotating, self-gravitating fluid systems. Here we investigate the secular and dynamical 3rd-harmonic instabilities that appear first on the Jacobi sequence of incompressible zero- vorticity ellipsoids. Poincare found a bifurcation point on the Jacobi sequence where a 3rd-harmonic mode becomes neutral. A sequence of pear-shaped equilibria branches off at this point but stands at higher energies. Therefore, the Jacobi ellipsoids remain secularly and dynamically stable. Cartan found that dynamical 3rd-harmonic instability also sets in at the Jacobi-pear bifurcation. We find that Cartan's instability leads to differentially rotating objects and not to uniformly rotating pear-shaped equilibria. We demonstrate that the pear-shaped sequence exists at higher energies and at higher rotation relative to the Jacobi sequence. The Jacobi ellipsoid undergoes fission on a secular time scale and a short-period binary is produced. The classical fission hypothesis of binary-star formation of Poincare and Darwin is thus feasible.

We use a free-energy minimization approach to describe the secular and dynamical instabilities as well as the bifurcations along equilibrium sequences of rotating, self-gravitating fluid systems. Our approach is fully nonlinear and stems from the Ginzburg-Landau theory of phase transitions. In this paper, we examine fourth-harmonic axisymmetric disturbances in Maclaurin spheroids and fourth-harmonic nonaxisymmetric disturbances in Jacobi ellipsoids. These two cases are very similar in the framework of phase transitions. Irrespective of whether a nonlinear first-order phase transition occurs between the critical point and the higher turning point or an apparent second-order phase transition occurs beyond the higher turning point, the result is fission (i.e. ``spontaneous breaking'' of the topology) of the original object on a secular time scale: the Maclaurin spheroid becomes a uniformly rotating axisymmetric torus and the Jacobi ellipsoid becomes a binary. The presence of viscosity is crucial since angular momentum needs to be redistributed for uniform rotation to be maintained. The phase transitions of the dynamical systems are briefly discussed in relation to previous numerical simulations of the formation and evolution of protostellar systems.

We have previously introduced the parameter `alpha' as an indicator of stability to m=2 nonaxisymmetric modes in rotating, self-gravitating, axisymmetric, gaseous and stellar systems. This parameter can be written as a function of the total rotational kinetic energy, the total gravitational potential energy, and as a function of the topology/connectedness and the geometric shape of a system. Here we extend the stability criterion to nonaxisymmetric equilibrium systems, such as ellipsoids and elliptical disks and cylinders. We test the validity of this extension by considering predictions for previously published, gaseous and stellar, nonaxisymmetric models. The above formulation and critical values account accurately for the stability properties of m=2 modes in gaseous Riemann S-type ellipsoids (including the Jacobi and Dedekind ellipsoids) and elliptical Riemann disks as well as in stellar elliptical Freeman disks and cylinders: all these systems are dynamically stable except for the stellar elliptical Freeman disks that exhibit a relatively small region of m=2 dynamical instability.

We analyze previous results on the stability of uniformly and differentially rotating, self-gravitating, gaseous and stellar, axisymmetric systems to derive a new stability criterion for the appearance of toroidal, m=2 Intermediate (I) and bar modes. In the process, we demonstrate that the bar modes in stellar systems and the m=2 I-modes in gaseous systems have many common physical characteristics and only one substantial difference: because of the anisotropy of the stress tensor, dynamical instability sets in at lower rotation in stellar systems. This difference is reflected also in the new stability criterion. The new stability parameter "alpha" is formulated first for uniformly rotating systems and is based on the angular momentum content rather than on the energy content of a system. For stability of stellar systems "alpha" = 0.254-0.258, while it is in the range of 0.341-0.354 for gaseous systems. For uniform rotation, one can write "alpha" as a function of the total (rotational) kinetic and gravitational energies, and of a parameter which is characteristic of the topology/connectedness and the geometric shape of a system. "Alpha" can be extended to and calculated for a variety of differentially rotating, gaseous and stellar, axisymmetric disk and spheroidal models whose equilibrium structures and stability characteristics are known. We also estimate "alpha" for gaseous toroidal models and for stellar disk systems embedded in an inert or responsive "halo". We find that the new stability criterion holds equally well for all these previously published axisymmetric models.

A free-energy minimization approach is used to address the secular & dynamical instabilities & the bifurcations along sequences of rotating, self-gravitating fluid and stellar systems. Our approach stems from the Landau-Ginzburg theory of phase transitions. We focus on the Maclaurin sequence of oblate spheroids & on the effects of second-harmonic disturbances. Second-order phase transitions appear on the Maclaurin sequence also at the points of dynamical instability & of bifurcation of the Dedekind sequence. Distinguishing characteristic of each second-order phase transition is the (non)conservation of an integral of motion (e.g. circulation) which determines the appearance of the transition. Circulation is not conserved in stellar systems because the stress-tensor gradient terms that appear in the Jeans equations of motion include viscosity-like off-diagonal terms of the same order of magnitude as the conventional pressure gradient terms. This explains why the Jacobi bifurcation is a point of dynamical instability in stellar systems but only a point of secular instability in viscous fluids. The second-order phase transitions are discussed in relation to the dynamical instability of stellar systems, the lambda-transition of liquid He-4, the second-order phase transition in superconductivity & the mechanism of spontaneous symmetry breaking.

We construct models of highly-inclined (``polar'') rings in an external potential including both self--gravity and dissipation due to a drag force. We adopt the oblate spheroidal scale--free logarithmic potential with axis ratio $q=0.85$ and an initial inclination of 80$^\circ$ for the self--gravitating rings. Stellar (dissipationless) rings suffer from mass loss during their evolution which drives a secular change of the mean inclination toward the poles of the potential. As much as half of the ring mass escapes in the process and forms an inner and an outer shell of precessing orbits. If the remaining mass is more than $\sim$ 0.02 of the enclosed galaxy mass, rings remain bound and and are not destroyed by differential precession. The rings precess at a constant rate for more than a precession period $\tau_p$ finding the configuration predicted by Sparke in 1986 which warps at larger radii toward the poles of the potential. We model shear viscosity with a velocity-dependent drag force and find that nuclear inflow dominates over self--gravity if the characteristic viscous inflow time scale $\tau_{vi}$ is shorter than $\sim 25 \tau_p$. Rings with $\tau_{vi}/\tau_p \ \sles \ 25$ collapse toward the nucleus of the potential within one precession period independent of the amount of self--gravity. Our results imply that stars and gas in real polar rings exhibit markedly different dynamical evolutions.