Showing 1–50 of 55 results for author: Christodoulou, D M

Measuring the magnetic dipole moment and magnetospheric fluctuations of accretion-powered pulsars in the Small Magellanic Cloud with an unscented Kalman filter

Authors: Joseph O'Leary, Andrew Melatos, Tom Kimpson, Nicholas J. O'Neill, Patrick M. Meyers, Dimitris M. Christodoulou, Sayantan Bhattacharya, Silas G. T. Laycock

Abstract: Many accretion-powered pulsars rotate in magnetocentrifugal disequilibrium, spinning up or down secularly over multi-year intervals. The magnetic dipole moment $μ$ of such systems cannot be inferred uniquely from the time-averaged aperiodic X-ray flux $\langle L(t) \rangle$ and pulse period $\langle P(t) \rangle$, because the radiative efficiency of the accretion is unknown and degenerate with the… ▽ More Many accretion-powered pulsars rotate in magnetocentrifugal disequilibrium, spinning up or down secularly over multi-year intervals. The magnetic dipole moment $μ$ of such systems cannot be inferred uniquely from the time-averaged aperiodic X-ray flux $\langle L(t) \rangle$ and pulse period $\langle P(t) \rangle$, because the radiative efficiency of the accretion is unknown and degenerate with the mass accretion rate. Here we circumvent the degeneracy by tracking the fluctuations in the unaveraged time series $L(t)$ and $P(t)$ using an unscented Kalman filter, whereupon $μ$ can be estimated uniquely, up to the uncertainties in the mass, radius and distance of the star. The analysis is performed on Rossi X-ray Timing Explorer observations for $24$ X-ray transients in the Small Magellanic Cloud, which have been monitored regularly for $\sim 16$ years. As well as independent estimates of $μ$, the analysis yields time-resolved histories of the mass accretion rate and the Maxwell stress at the disk-magnetosphere boundary for each star, and hence auto- and cross-correlations involving the latter two state variables. The inferred fluctuation statistics convey important information about the complex accretion physics at the disk-magnetosphere boundary. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 42 pages, 10 figures, accepted for publication in The Astrophysical Journal

Measuring the magnetic dipole moment and magnetospheric fluctuations of SXP 18.3 with a Kalman filter

Authors: J. O'Leary, A. Melatos, N. J. O'Neill, P. M. Meyers, D. M. Christodoulou, S. Bhattacharya, S. G. T. Laycock

Abstract: The magnetic dipole moment $μ$ of an accretion-powered pulsar in magnetocentrifugal equilibrium cannot be inferred uniquely from time-averaged pulse period and aperiodic X-ray flux data, because the radiative efficiency $η_0$ of the accretion is unknown, as are the mass, radius, and distance of the star. The degeneracy associated with the radiative efficiency is circumvented, if fluctuations of th… ▽ More The magnetic dipole moment $μ$ of an accretion-powered pulsar in magnetocentrifugal equilibrium cannot be inferred uniquely from time-averaged pulse period and aperiodic X-ray flux data, because the radiative efficiency $η_0$ of the accretion is unknown, as are the mass, radius, and distance of the star. The degeneracy associated with the radiative efficiency is circumvented, if fluctuations of the pulse period and aperiodic X-ray flux are tracked with a Kalman filter, whereupon $μ$ can be measured uniquely up to the uncertainties in the mass, radius, and distance. Here the Kalman filter analysis is demonstrated successfully in practice for the first time on Rossi X-ray Timing Explorer observations of the X-ray transient SXP 18.3 in the Small Magellanic Cloud, which is monitored regularly. The analysis yields $μ= 8.0^{+1.3}_{-1.2} \, \times \, 10^{30} \, {\rm G \, cm^3}$ and $η_0 = 0.04^{+0.02}_{-0.01}$, compared to $μ= 5.0^{+1.0}_{-1.0} \times 10^{30} \, {\rm G \, cm^3}$ as inferred traditionally from time-averaged data assuming $η_0=1$. The analysis also yields time-resolved estimates of two hidden state variables, the mass accretion rate and the Maxwell stress at the disk-magnetosphere boundary. The success of the demonstration confirms that the Kalman filter analysis can be applied in the future to study the magnetic moments and disk-magnetosphere physics of accretion-powered pulsar populations in the Small Magellanic Cloud and elsewhere. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 34 pages, 5 figures, accepted for publication in The Astrophysical Journal

The First TESS Self-Lensing Pulses: Revisiting KIC 12254688

Authors: Nicholas M. Sorabella, Silas G. T. Laycock, Dimitris M. Christodoulou, Sayantan Bhattacharya

Abstract: We report the observations of two self-lensing pulses from KIC 12254688 in Transiting Exoplanet Survey Satellite (TESS) light curves. This system, containing a F2V star and white-dwarf companion, was amongst the first self-lensing binary systems discovered by the Kepler Space Telescope over the past decade. Each observed pulse occurs when the white dwarf transits in front of its companion star, gr… ▽ More We report the observations of two self-lensing pulses from KIC 12254688 in Transiting Exoplanet Survey Satellite (TESS) light curves. This system, containing a F2V star and white-dwarf companion, was amongst the first self-lensing binary systems discovered by the Kepler Space Telescope over the past decade. Each observed pulse occurs when the white dwarf transits in front of its companion star, gravitationally lensing the star's surface, thus making it appear brighter to a distant observer. These two pulses are the very first self-lensing events discovered in TESS observations. We describe the methods by which the data were acquired and detrended, as well as the best-fit binary parameters deduced from our self-lensing+radial velocity model. We highlight the difficulties of finding new self-lensing systems with TESS, and we discuss the types of self-lensing systems that TESS may be more likely to discover in the future. △ Less

Submitted 15 January, 2024; v1 submitted 2 January, 2024; originally announced January 2024.

Comments: Accepted for publication in ApJ Letters, 12 Pages, 6 Figures, and 1 Table

Accreting Black Holes Skewing and Bending the Optical Emission from Massive Wolf-Rayet Companions -- A Case Study of IC10 X-1

Authors: Sayantan Bhattacharya, Dimitris M. Christodoulou, Andre-Nicolas Chene, Silas G. T. Laycock, Breanna A. Binder, Demosthenes Kazanas

Abstract: We present a statistical analysis of the He ii 4686 emission line in the spectra of the black hole and Wolf-Rayet (WR) star of the high-mass X-ray binary IC10 X-1. This line is visibly skewed, and the third moment (skewness) varies with the binary's orbital phase. We describe a new method of extracting such weak/faint features lying barely above a noisy continuum. Using the moments of these featur… ▽ More We present a statistical analysis of the He ii 4686 emission line in the spectra of the black hole and Wolf-Rayet (WR) star of the high-mass X-ray binary IC10 X-1. This line is visibly skewed, and the third moment (skewness) varies with the binary's orbital phase. We describe a new method of extracting such weak/faint features lying barely above a noisy continuum. Using the moments of these features, we have been able to decompose these skewed lines into two symmetric Gaussian profiles as a function of the orbital phase. The astrophysical implications of this decomposition are significant due to the complex nature of wind-accretion stream interactions in such binary systems. Previous studies have already shown a 0.25 phase lag in the radial velocity curve of the star and the X-ray eclipse, which indicates that the He ii emitters might be in the stellar wind, hence not tracing the star's orbital motion. Results from this work further suggest the existence of two separate emitting regions, one in the stellar wind in the shadow of the WR star, and another in the accretion stream that impacts the black hole's outer accretion disk; and the observed skewed He ii lines can be reproduced by superposition of the two corresponding time-dependent Gaussian emission profiles. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: 13 pages, 5 figures, accepted for publication in MNRAS Journal

Probing the Stellar Wind of the Wolf-Rayet Star in IC 10 X-1

Authors: Sayantan Bhattacharya, Silas G. T. Laycock, Andre-Nicolas Chene, Breanna A. Binder, Dimitris M. Christodoulou, Ankur Roy, Nicholas M. Sorabella, Rigel C. Cappallo

Abstract: IC 10 X-1 is an eclipsing high mass X-ray binary (HMXB) containing a stellar-mass black hole (BH) and a Wolf-Rayet (WR) donor star with an orbital period of P = 34.9 hr. This binary belongs to a group of systems that can be the progenitors of gravitational wave sources, hence understanding the dynamics of systems such as IC 10 X-1 is of paramount importance. The prominent He II 4686 emission line… ▽ More IC 10 X-1 is an eclipsing high mass X-ray binary (HMXB) containing a stellar-mass black hole (BH) and a Wolf-Rayet (WR) donor star with an orbital period of P = 34.9 hr. This binary belongs to a group of systems that can be the progenitors of gravitational wave sources, hence understanding the dynamics of systems such as IC 10 X-1 is of paramount importance. The prominent He II 4686 emission line (previously used in mass estimates of the BH) is out of phase with the X-ray eclipse, suggesting that this line originates somewhere in the ionized wind of the WR star or in the accretion disk. We obtained 52 spectra from the GEMINI/GMOS archive, observed between 2001 and 2019. We analyzed the spectra both individually, and after binning them by orbital phase to improve the signal-to-noise ratio. The RV curve from the stacked data is similar to historical results, indicating the overall parameters of the binary have remained constant. However, the He II line profile shows a correlation with the X-ray hardness-ratio values, also, we report a pronounced skewness of the line-profile, and the skewness varies with the orbital phase. These results support a paradigm wherein the He II line tracks structures in the stellar wind that are produced by interactions with the BH's ionizing radiation and the accretion flow. We compare the observable signatures of two alternative hypotheses proposed in the literature: wind irradiation plus shadowing, and accretion disk hotspot; and we explore how the line-profile variations fit into each of these models. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Comments: 10 pages, 5 figures, 2 tables, ApJ Accepted

Disentangling the neighbouring pulsars SXP 15.3 and SXP 305

Authors: Itumeleng M. Monageng, Malcolm J. Coe, Lee J. Townsend, Silas G. T. Laycock, Jamie A. Kennea, Ankur Roy, Andrzej Udalski, Sayantan Bhattacharya, Dimitris M. Christodoulou, David A. H. Buckley, Phil A. Evans

Abstract: SXP 15.3 and SXP 305 are two Be X-ray binaries in the Small Magellanic Cloud that are spatially separated by ~7 arcsec. The small separation between these sources has, in the past, resulted in confusion about the origin of the emission from the combined region. We present long-term optical and X-ray monitoring results of both sources, where we study the historic and recent behaviour. In particular… ▽ More SXP 15.3 and SXP 305 are two Be X-ray binaries in the Small Magellanic Cloud that are spatially separated by ~7 arcsec. The small separation between these sources has, in the past, resulted in confusion about the origin of the emission from the combined region. We present long-term optical and X-ray monitoring results of both sources, where we study the historic and recent behaviour. In particular, from data collected as part of the S-CUBED project we see repeating X-ray outbursts from the combined region of the two sources in the recent lightcurve from the Neil Gehrels Swift Observatory, and we investigate the origin of this emission. Using the H-alpha emission line from the Southern African Large Telescope (SALT) and photometric flux from the Optical Gravitational Lensing Experiment (OGLE) to study the changes in the size and structure of the Be disc, we demonstrate that the X-ray emission likely originates from SXP 15.3. Timing analysis reveals unusual behaviour, where the optical outburst profile shows modulation at twice the frequency of the X-ray outbursts. We consider either of these periodicities being the true orbital period in SXP 15.3 and propose models based on the geometric orientations of the Be disc and neutron star to explain the physical origin of the outbursts. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS

The Wolf-Rayet + Black Hole Binary NGC 300 X-1: What is the Mass of the Black Hole?

Authors: Breanna A. Binder, Janelle M. Sy, Michael Eracleous, Dimitris M. Christodoulou, Sayantan Bhattacharya, Rigel Cappallo, Silas Laycock, Paul P. Plucinsky, Benjamin F. Williams

Abstract: We present new X-ray and UV observations of the Wolf-Rayet + black hole binary system NGC 300 X-1 with the Chandra X-ray Observatory and the Hubble Space Telescope Cosmic Origins Spectrograph. When combined with archival X-ray observations, our X-ray and UV observations sample the entire binary orbit, providing clues to the system geometry and interaction between the black hole accretion disk and… ▽ More We present new X-ray and UV observations of the Wolf-Rayet + black hole binary system NGC 300 X-1 with the Chandra X-ray Observatory and the Hubble Space Telescope Cosmic Origins Spectrograph. When combined with archival X-ray observations, our X-ray and UV observations sample the entire binary orbit, providing clues to the system geometry and interaction between the black hole accretion disk and the donor star wind. We measure a binary orbital period of 32.7921$\pm$0.0003 hr, in agreement with previous studies, and perform phase-resolved spectroscopy using the X-ray data. The X-ray light curve reveals a deep eclipse, consistent with inclination angles of $i=60-75^{\circ}$, and a pre-eclipse excess consistent with an accretion stream impacting the disk edge. We further measure radial velocity variations for several prominent FUV spectral lines, most notably He II $λ$1640 and C IV $λ$1550. We find that the He II emission lines systematically lag the expected Wolf-Rayet star orbital motion by a phase difference $Δφ\sim0.3$, while C IV $λ$1550 matches the phase of the anticipated radial velocity curve of the Wolf-Rayet donor. We assume the C IV $λ$1550 emission line follows a sinusoidal radial velocity curve (semi-amplitude = 250 km s$^{-1}$) and infer a BH mass of 17$\pm$4 M$_{\odot}$. Our observations are consistent with the presence of a wind-Roche lobe overflow accretion disk, where an accretion stream forms from gravitationally focused wind material and impacts the edge of the black hole accretion disk. △ Less

Submitted 13 February, 2021; originally announced February 2021.

Comments: Accepted to ApJ. 22 pages, 16 figures, 7 tables

On the observed clustering of major bodies in solar and extrasolar subsystems

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: Major (exo)planetary and satellite bodies seem to concentrate at intermediate areas of the radial distributions of all the objects present in each (sub)system. We prove rigorously that the secular evolution of (exo)planets and satellites necessarily results in the observed intermediate accumulation of the massive objects in all such subsystems. We quantify a "middle" as the mean of mean motions (o… ▽ More Major (exo)planetary and satellite bodies seem to concentrate at intermediate areas of the radial distributions of all the objects present in each (sub)system. We prove rigorously that the secular evolution of (exo)planets and satellites necessarily results in the observed intermediate accumulation of the massive objects in all such subsystems. We quantify a "middle" as the mean of mean motions (orbital angular velocities) of three or more massive objects involved. Orbital evolution is expected to be halted or severely diminished when the survivors settle near mean-motion resonances and substantial angular-momentum transfer between bodies ceases to occur (gravitational Landau damping). The dynamics is opposite in direction to what has been theorized for viscous and magnetized accretion disks in which gas spreads out and away from either side of any conceivable intermediate area. The results are bound to change the way we think about planet and moon formation and evolution. △ Less

Submitted 20 January, 2021; v1 submitted 6 January, 2021; originally announced January 2021.

Comments: Submitted to RAA

Meta-analysis of Electron Cyclotron Resonance Absorption Features Detected in High-Mass X-ray Binaries

Authors: Dimitris M. Christodoulou, Silas G. T. Laycock, Demosthenes Kazanas

Abstract: 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 no… ▽ More 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. △ Less

Submitted 16 May, 2019; v1 submitted 13 May, 2019; originally announced May 2019.

Comments: Submitted to RAA. This version incorporates the referee's comments and matches the published version

Theoretical model of HD 163296 presently forming in-situ planets and comparison with the models of AS 209, HL Tau, and TW Hya

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 A… ▽ More 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 $ρ(R)\propto R^{-1}$), the rotational parameter $β_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 $Ω_J$ (to within factors of 1.4); the equation of state ($c_0^2/ρ_0$) and the central density $ρ_0$ (to within factors of 2); and the core angular velocity $Ω_0$ (to within a factor of 4.5). In the end, we compare all six nebular disks that we have modeled so far. △ Less

Submitted 3 March, 2019; originally announced March 2019.

Comments: Part 10 and the final part (HD 163296)

Theoretical model of the outer disk of TW Hya presently forming in-situ planets and comparison with models of AS 209 and HL Tau

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 (… ▽ More 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 $ρ(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 $Ω_J$ and the angular velocity $Ω_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 ($ρ_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 ($Ω_J$). In the end, we compare the protostellar disks that we have modeled so far. △ Less

Submitted 3 March, 2019; v1 submitted 9 February, 2019; originally announced February 2019.

Comments: Updated version, Part 9 (TW Hya). arXiv admin note: text overlap with arXiv:1901.10642

Preliminary model of the outer disk of RU Lup presently showing only four dark gaps

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 observ… ▽ More 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 ($ρ_0$), centrifugal support ($β_0$), and core angular velocity ($Ω_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. △ Less

Submitted 3 March, 2019; v1 submitted 31 January, 2019; originally announced February 2019.

Comments: Updated version, Part 8 (RU Lup). arXiv admin note: text overlap with arXiv:1901.10642

Theoretical models of the protostellar disks of AS 209 and HL Tau presently forming in-situ planets

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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… ▽ More 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 $ρ(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. △ Less

Submitted 3 March, 2019; v1 submitted 29 January, 2019; originally announced January 2019.

Comments: Updated version, Part 7 (AS 209 and HL Tau)

Models of Saturn's protoplanetary disk forming in-situ its regular satellites and innermost rings before the planet is formed

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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. It… ▽ More 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 ($ΔR\lesssim 0.9$ Gm), its rotation parameter that measures centrifugal support against self-gravity is somewhat larger ($β_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. △ Less

Submitted 3 March, 2019; v1 submitted 22 January, 2019; originally announced January 2019.

Comments: Updated version, Part 5 (Saturn). arXiv admin note: substantial text overlap with arXiv:1901.06448; text overlap with arXiv:1901.05131

Model of Neptune's protoplanetary disk forming in-situ its surviving regular satellites after Triton's capture and comparison of the protoplanetary disks of the four gaseous giants

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 th… ▽ More 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 ($β_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. △ Less

Submitted 3 March, 2019; v1 submitted 25 January, 2019; originally announced January 2019.

Comments: Updated version, Part 6 (Neptune). arXiv admin note: text overlap with arXiv:1901.06448

Models of a protoplanetary disk forming in-situ the major Uranian satellites before the planet is formed

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 measure… ▽ More 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 ($β_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. △ Less

Submitted 3 March, 2019; v1 submitted 18 January, 2019; originally announced January 2019.

Comments: Updated version, Part 4 (Uranus). arXiv admin note: text overlap with arXiv:1901.05131

Models of a protoplanetary disk forming in-situ the Galilean and smaller nearby satellites before Jupiter is formed

Authors: D. M. Christodoulou, D. Kazanas

Abstract: 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,… ▽ More 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. △ Less

Submitted 3 March, 2019; v1 submitted 15 January, 2019; originally announced January 2019.

Comments: Updated version, Part 3 (Jupiter)

Conundrums and constraints concerning the formation of our solar system -- An alternative view

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 constr… ▽ More 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. △ Less

Submitted 3 March, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: Updated version, Part 2 (Small Bodies in the Solar System)

On the formation of our solar system and many other protoplanetary systems observed by ALMA and SPHERE

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 t… ▽ More 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. △ Less

Submitted 3 March, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: Updated version, Part 1 (Solar Nebula). arXiv admin note: text overlap with arXiv:0706.3205

Gauss's Law and the Source for Poisson's Equation in Modified Gravity with Varying G

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 differe… ▽ More 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. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: MNRAS, submitted

Gravitational Potential and Nonrelativistic Lagrangian in Modified Gravity with Varying G

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 imp… ▽ More 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. △ Less

Submitted 21 November, 2018; originally announced November 2018.

Comments: To appear in MNRASL

Interposing a Varying Gravitational Constant Between Modified Newtonian Dynamics and Weak Weyl Gravity

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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… ▽ More 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. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: An original letter to appear in MNRAS

Not an Oxymoron: Some X-ray Binary Pulsars with Enormous Spinup Rates Reveal Weak Magnetic Fields

Authors: D. M. Christodoulou, S. G. T. Laycock, D. Kazanas

Abstract: 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. Th… ▽ More 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. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: Submitted to MNRAS

A Striking Confluence Between Theory and Observations of High-Mass X-ray Binary Pulsars

Authors: D. M. Christodoulou, S. G. T. Laycock, D. Kazanas

Abstract: 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 ou… ▽ More 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. △ Less

Submitted 9 January, 2018; v1 submitted 12 December, 2017; originally announced December 2017.

Comments: Added LMC X-4 and commented on the cyclotron absorption line of SMC X-2. 4 pages, 1 figure, 2 tables, submitted to MNRASL

Golden Elliptical Orbits in Newtonian Gravitation

Authors: Dimitris M. Christodoulou

Abstract: 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 leng… ▽ More 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. △ Less

Submitted 27 November, 2017; originally announced November 2017.

Comments: Submitted to Forum Geometricorum, 2 figures for the two golden ellipses of Newtonian dynamics. Based on the results of arXiv:1705.09356

New Bound Closed Orbits in Spherical Potentials

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 e… ▽ More 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. △ Less

Submitted 28 September, 2017; v1 submitted 16 July, 2017; originally announced July 2017.

Comments: There are issues to be resolved with the dynamical equations that determine the orbits

A Physical Interpretation of the Titius-Bode Rule and its Connection to the Closed Orbits of Bertrand's Theorem

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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… ▽ More 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. △ Less

Submitted 11 September, 2017; v1 submitted 25 May, 2017; originally announced May 2017.

Comments: Final version, issues resolved, submitted to RAA

Retrograde Accretion Disks in High-Mass Be/X-ray Binaries

Authors: D. M. Christodoulou, S. G. T. Laycock, D. Kazanas

Abstract: 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 differe… ▽ More 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. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: To appear in MNRAS Letters

The Great Pretenders Among the ULX Class

Authors: D. M. Christodoulou, S. G. T. Laycock, D. Kazanas, R. Cappallo, I. Contopoulos

Abstract: 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 directi… ▽ More 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. △ Less

Submitted 24 March, 2017; originally announced March 2017.

Comments: To appear in Research in Astronomy and Astrophysics

A Comprehensive Library of X-ray Pulsars in the Small Magellanic Cloud: Time Evolution of their Luminosities and Spin Periods

Authors: J. Yang, S. G. T. Laycock, D. M. Christodoulou, S. Fingerman, M. J. Coe, J. J. Drake

Abstract: 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… ▽ More 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. △ Less

Submitted 15 March, 2017; originally announced March 2017.

Comments: 17 pages, 11 + 58 (appendix) figures. To appear in the Astrophysical Journal Supplement

Determination of the Magnetic Fields of Magellanic X-Ray Pulsars

Authors: Dimitris M. Christodoulou, Silas G. T. Laycock, Jun Yang, Samuel Fingerman

Abstract: 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 observa… ▽ More 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). △ Less

Submitted 10 March, 2017; originally announced March 2017.

Comments: To appear in Reserch in Astronomy and Astrophysics

Exact Axisymmetric Solutions of the 2-D Lane-Emden Equations with Rotation

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 d… ▽ More 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. △ Less

Submitted 14 January, 2017; originally announced January 2017.

Journal ref: Journal of Modern Physics, 2016, 7, 2177-2187

Blue Supergiant X-ray Binaries in the Nearby Dwarf Galaxy IC 10

Authors: Silas G. T. Laycock, Dimitris M. Christodoulou, Benjamin F. Williams, Breanna Binder, Andrea Prestwich

Abstract: 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, derive… ▽ More 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. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: Accepted to ApJ, 7 pages, 5 figures

The X-ray Binary Population of the Nearby Dwarf Starburst Galaxy IC 10: Variable and Transient X-ray Sources

Authors: Silas G. T. Laycock, Rigel Cappallo, Benjamin F. Williams, Andrea Prestwich, Breanna Binder, Dimitris M. Christodoulou

Abstract: 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… ▽ More 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. △ Less

Submitted 25 November, 2016; originally announced November 2016.

Comments: 17 pages, 8 figures. Accepted to ApJ

The X-ray Pulsar M82 X-2 on its Propeller Line

Authors: D. M. Christodoulou, D. Kazanas, S. G. T. Laycock

Abstract: {\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 b… ▽ More {\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. △ Less

Submitted 22 June, 2016; originally announced June 2016.

Comments: 7 pages, submitted to MNRAS. This is a new calculation of the magnetic field of M82 X-2 based on the latest observations from the Chandra archive

On the Theoretical Framework of Magnetized Outflows from Stellar-Mass Black Holes and Related Observations

Authors: D. M. Christodoulou, I. Contopoulos, D. Kazanas, J. F. Steiner, D. B. Papadopoulos, S. G. T. Laycock

Abstract: 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 wh… ▽ More 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. △ Less

Submitted 21 June, 2016; originally announced June 2016.

Comments: 9 pages, 8 figures, 2 tables. Accepted for publication to MNRAS

Dominance of outflowing electric currents on decaparsec to kiloparsec scales in extragalactic jets

Authors: Dimitris M. Christodoulou, Denise C. Gabuzda, Sebastian Knuettel, Ioannis Contopoulos, Demosthenes Kazanas, Colm P. Coughlan

Abstract: 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 correlat… ▽ More 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. △ Less

Submitted 29 May, 2016; originally announced May 2016.

Comments: 12 pages, 4 figures, 1 table, accepted for publication in Astronomy & Astrophysics

The Case Against Dark Matter and Modified Gravity: Flat Rotation Curves Are a Rigorous Requirement in Rotating Self-Gravitating Newtonian Gaseous Disks

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 "f… ▽ More 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. △ Less

Submitted 24 October, 2015; v1 submitted 19 October, 2015; originally announced October 2015.

Comments: Appendix was added to address vigorous objections against using cylindrical coordinates to model the equatorial planes of vertically thin disks

Journal ref: Journal of Modern Physics, 7, 680-698 (2016)

Revisiting the Dynamical Case for a Massive Black Hole in IC10 X-1

Authors: Silas G. T. Laycock, Thomas J. Maccarone, Dimitris M. Christodoulou

Abstract: 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[$λ$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… ▽ More 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[$λ$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. △ Less

Submitted 11 June, 2015; originally announced June 2015.

Comments: Accepted for Publication in MNRAS

The Magnetic Field of the Ultraluminous X-ray Pulsar M82 X-2

Authors: Dimitris M. Christodoulou, Silas G. T. Laycock, Demosthenes Kazanas

Abstract: 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 o… ▽ More 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. △ Less

Submitted 21 June, 2016; v1 submitted 19 November, 2014; originally announced November 2014.

Comments: This paper has been withdrawn by the authors. New observations have led to a different calculation that will be presented in another paper

Formation and destruction of jets in X-ray binaries

Authors: N. D. Kylafis, I. Contopoulos, D. Kazanas, D. M. Christodoulou

Abstract: 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 har… ▽ More 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. △ Less

Submitted 20 December, 2011; v1 submitted 8 April, 2011; originally announced April 2011.

Comments: Accepted by A&A

The invariant Twist of Magnetic Fields in the Relativistic Jets of Active Galactic Nuclei

Authors: Ioannis Contopoulos, Dimitris M. Christodoulou, Demosthenes Kazanas, Denise C. Gabuzda

Abstract: 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 accretio… ▽ More 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. △ Less

Submitted 21 July, 2009; originally announced July 2009.

Comments: Accepted for publication in the Astrophysical Journal Letters

55 Cancri: A Laboratory for Testing Numerous Conjectures about Planet Formation

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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 s… ▽ More 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 $Σ(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. △ Less

Submitted 5 November, 2008; originally announced November 2008.

Comments: The physical model applied to the solar system in 0706.3205 (updated version 2) is now applied to the 5 planets of 55 Cancri. 5 tables, 2 figures

Exact Solutions of the Isothermal Lane--Emden Equation with Rotation and Implications for the Formation of Planets and Satellites

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: 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,… ▽ More 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 ... △ Less

Submitted 5 November, 2008; v1 submitted 21 June, 2007; originally announced June 2007.

Comments: Modified version with longer historical introduction, discussion of model stability, and updated discussion of multi-planet extrasolar systems

Simulations of the Poynting--Robertson Cosmic Battery in Resistive Accretion Disks

Authors: Dimitris M. Christodoulou, Ioannis Contopoulos, Demosthenes Kazanas

Abstract: 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… ▽ More 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. △ Less

Submitted 13 June, 2007; originally announced June 2007.

Comments: 27 pages, 55 figures

The Cosmic Battery Revisited

Authors: Ioannis Contopoulos Demosthenes Kazanas Dimitris M. Christodoulou

Abstract: We reinvestigate the generation and accumulation of magnetic flux in optically thin accretion flows around active gravitating objects. The source of the magnetic field is the azimuthal electric current associated with the Poynting-Robertson drag on the electrons of the accreting plasma. This current generates magnetic field loops which open up because of the differential rotation of the flow. We… ▽ More We reinvestigate the generation and accumulation of magnetic flux in optically thin accretion flows around active gravitating objects. The source of the magnetic field is the azimuthal electric current associated with the Poynting-Robertson drag on the electrons of the accreting plasma. This current generates magnetic field loops which open up because of the differential rotation of the flow. We show through simple numerical simulations that what regulates the generation and accumulation of magnetic flux near the center is the value of the plasma conductivity. Although the conductivity is usually considered to be effectively infinite for the fully ionized plasmas expected near the inner edge of accretion disks, the turbulence of those plasmas may actually render them much less conducting due to the presence of anomalous resistivity. We have discovered that if the resistivity is sufficiently high throughout the turbulent disk while it is suppressed interior to its inner edge, an interesting steady-state process is established: accretion carries and accumulates magnetic flux of one polarity inside the inner edge of the disk, whereas magnetic diffusion releases magnetic flux of the opposite polarity to large distances. In this scenario, magnetic flux of one polarity grows and accumulates at a steady rate in the region inside the inner edge and up to the point of equipartition when it becomes dynamically important. We argue that this inward growth and outward expulsion of oppositely-directed magnetic fields that we propose may account for the ~30 min cyclic variability observed in the galactic microquasar GRS1915+105. △ Less

Submitted 31 August, 2006; originally announced August 2006.

Journal ref: Astrophys.J.652:1451-1456,2006

Interchange Method in Compressible Magnetized Couette Flow: Magnetorotational and Magnetoconvective Instabilities

Authors: Dimitris M. Christodoulou, John Contopoulos, Demosthenes Kazanas

Abstract: 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, entr… ▽ More 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. △ Less

Submitted 22 November, 2002; originally announced November 2002.

Comments: 24 pages, AASTEX, accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J. 586 (2003) 372-383

Infrared Search for Young Stars in HI High-velocity Clouds

Authors: Z. Ivezic, D. M. Christodoulou

Abstract: 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… ▽ More 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. △ Less

Submitted 25 May, 1997; originally announced May 1997.

Comments: 9 pages, 4 JPEG figures. PostScript figures available from authors

Phase-Transition theory of Instabilities. IV. Critical Points on the Maclaurin Sequence and Nonlinear Fission Processes

Authors: D. M. Christodoulou, D. Kazanas, I. Shlosman, J. E. Tohline

Abstract: 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 bi… ▽ More 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. △ Less

Submitted 22 May, 1995; originally announced May 1995.

Comments: 23 pages, postscript, compressed, uuencoded. Figs. (6) available by anonymous ftp from ftp://asta.pa.uky.edu/shlosman/paper4/ , get *.ps.Z). To appear in ApJ

Phase-Transition Theory of Instabilities. III. The Third-Harmonic Bifurcation on the Jacobi Sequence and the Fission Problem

Authors: D. M. Christodoulou, D. Kazanas, I. Shlosman, J. E. Tohline

Abstract: 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… ▽ More 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. △ Less

Submitted 2 May, 1995; originally announced May 1995.

Comments: 22 pages (including 2 figures), postscript, compressed, uuencoded. Also available by anonymous ftp at ftp://asta.pa.uky.edu/shlosman/paper3 , mget *.ps.Z). To appear in ApJ