Regularized Einstein-Gauss-Bonnet (EGB) theory of gravity in four dimensions is a new attempt to include nontrivial contributions of Gauss-Bonnet term. In this paper, we make a detailed analysis on possible constraints of the model parameters of the theory from recent cosmological observations, and some theoretical constraints as well. Our results show that the theory with vanishing bare cosmological constant, $\Lambda_0$, is ruled out by the current observational value of $w_{de}$, and the observations of GW170817 and GRB 170817A as well. For nonvanishing bare cosmological constant, instead, our results show that the current observation of the speed of GWs measured by GW170817 and GRB 170817A would place a constraints on $\tilde{\alpha}$, a dimensionless parameter of the theory, as $-7.78\times10^{-16}\le \tilde{\alpha} \leq 3.33\times10^{-15}$.

We show an inconsistence of the novel 4D Einstein-Gauss-Bonnet gravity by considering a quantum tunneling process of vacua. Using standard semi-classical techniques, we analytically study the vacuum decay rate for all allowed cases in the parameter space. It turns out, without exception, that the theory either encounters a disastrous divergence of vacuum decay rate, or exhibits a confusing complex value of vacuum decay rate, or involves an infinitely instability (a large vacuum mixing). These suggest a strong possibility that the theory, at least the vacuum of the theory, is either unphysical or unstable, or has no well-defined limit as $D\rightarrow 4$.

Geodetic Very Long Baseline Interferometry (VLBI) data are capable of measuring the light deflection caused by the gravitational field of the Sun and large planets with high accuracy. The parameter $\gamma$ of the parametrized Post-Newtonian (PPN) formalism estimated using observations of reference radio sources near the Sun should be equal to unity in the general relativity. We have run several VLBI experiments tracking reference radio sources from 1 to 3 degrees from the Sun. The best formal accuracy of the parameter $\gamma$ achieved in the single-session mode is less than 0.01 percent, or better than the formal accuracy obtained with a global solution included all available observations at arbitrary elongation from the Sun. We are planning more experiments starting from 2020 using better observing conditions near the minimum of the Solar activity cycle.

Context. Radio-loud active galactic nuclei in the early Universe are rare. The quasars J0906+6930 at redshift z=5.47 and J2102+6015 at z=4.57 stand out from the known sample with their compact emission on milliarcsecond (mas) angular scale with high (0.1-Jy level) flux densities measured at GHz radio frequencies. This makes them ideal targets for very long baseline interferometry (VLBI) observations. Aims. By means of VLBI imaging we can reveal the inner radio structure of quasars and model their brightness distribution to better understand the geometry of the jet and the physics of the sources. Methods. We present sensitive high-resolution VLBI images of J0906+6930 and J2102+6015 at two observing frequencies, 2.3 and 8.6 GHz. The data were taken in an astrometric observing programme involving a global five-element radio telescope array. We combined the data from five different epochs from 2017 February to August. Results. For one of the highest redshift blazars known, J0906+6930, we present the first-ever VLBI image obtained at a frequency below 8 GHz. Based on our images at 2.3 and 8.6 GHz, we confirm that this source has a sharply bent helical inner jet structure within ~3 mas from the core. The quasar J2102+6015 shows an elongated radio structure in the east-west direction within the innermost ~2 mas that can be described with a symmetric three-component brightness distribution model at 8.6 GHz. Because of their non-pointlike mas-scale structure, these sources are not ideal as astrometric reference objects. Our results demonstrate that VLBI observing programmes conducted primarily with astrometric or geodetic goals can be utilized for astrophysical purposes as well.

Context. We highlight the capabilities of the geodetic VLBI technique to test General relativity in the classical astrometric style, i.e., measuring the deflection of light in the vicinity of the Sun. Aims. In previous studies, the parameter was estimated by global analyses of thousands of geodetic VLBI sessions. Here we estimate from a single session where the Sun has approached two strong reference radio sources 0229+131 and 0235+164 at an elongation angle of 1-3 degrees. Methods. The AUA020 VLBI session of 1 May 2017 was designed to obtain more than 1000 group delays from the two radio sources. The Solar corona effect was effectively calibrated with the dual-frequency observations even at small elongation from the Sun. Results. We obtained with a precision better than what is obtained through global analyses of thousands of standard geodetic sessions over decades. Current results demonstrate that the modern VLBI technology is capable of establishing new limits on observational test of General Relativity.

Shanghai Astronomical Observatory has upgraded its DiFX cluster to 420 CPU cores and a 432-TB storage system at the end of 2014. An international network connection for the raw data transfer has also been established. The routine operations for IVS sessions including CRF, AOV, and APSG series began in early 2015. In addition to the IVS observations, the correlator is dedicated to astrophysical and astrometric programs with the Chinese VLBI Network and international joint VLBI observations. It also worked with the new-built Tianma 65-m radio telescope and successfully found fringes as high as at X/Ka and Q bands in late 2015. A more powerful platform is planned for the high data rate and massive data correlation tasks in the future.

We present here the results of the first part of the VLBI Ecliptic Plane Survey (VEPS) program. The goal of the program is to find all compact sources within $7.5^\circ$ of the ecliptic plane which are suitable as calibrators for anticipated phase referencing observations of spacecraft and determine their positions with accuracy at the 1.5~nrad level. We run the program in two modes: the search mode and the refining mode. In the search mode, a complete sample of all sources brighter than 50 mJy at 5 GHz listed in the Parkes-MIT-NRAO (PMN) and Green Bank 6~cm (GB6) catalogs, except those previously detected with VLBI, is observed. In the refining mode, the positions of all ecliptic plane sources, including those found in the search mode, are improved. By October 2016, thirteen 24-hr sessions that targeted all sources brighter than 100~mJy have been observed and analyzed. Among 3320 observed target sources, 555 objects have been detected. We also conducted a number of follow-up VLBI experiments in the refining mode and improved the positions of 249 ecliptic plane sources.

We launched the VLBI Ecliptic Plane Survey program in 2015. The goal of this program is to find all compact sources within 7.5 degrees of the ecliptic plane which are suitable as phase calibrators for anticipated phase referencing observations of spacecrafts. We planned to observe a complete sample of the sources brighter than 50 mJy at 5 GHz listed in the PMN and GB6 catalogues that have not yet been observed with VLBI. By April 2016, eight 24-hour sessions have been performed and processed. Among 2227 observed sources, 435 sources were detected in three or more observations. We have also run three 8-hour segments with VLBA for improving positions of 71 ecliptic sources.

In this paper, we study static vacuum solutions of quantum gravity at a fixed Lifshitz point in (2+1) dimensions, and present all the diagonal solutions in closed forms in the infrared limit. The exact solutions represent spacetimes with very rich structures: they can represent generalized BTZ black holes, Lifshitz space-times or Lifshitz solitons, in which the spacetimes are free of any kind of space-time singularities, depending on the choices of the free parameters of the solutions. We also find several classes of exact static non-diagonal solutions, which represent similar space-time structures as those given in the diagonal case. The relevance of these solutions to the non-relativistic Lifshitz-type gauge/gravity duality is discussed.

In this paper, we show that the spin-0 gravitons appearing in Horava-Lifshitz gravity without the projectability condition can be eliminated by extending the gauge symmetries of the foliation-preserving diffeomorphisms to include a local U(1) symmetry. As a result, the problems of stability, ghost, strong coupling, and different speeds in the gravitational sector are automatically resolved. In addition, with the detailed balance condition softly breaking, the number of independent coupling constants can be significantly reduced (from more than 70 down to 15), while the theory is still UV complete and possesses a healthy IR limit, whereby the prediction powers of the theory are considerably improved. The strong coupling problem in the matter sector can be cured by introducing an energy scale $M_{*}$, so that $M_{*} < \Lambda_{\omega}$, where $M_{*}$ denotes the suppression energy of high order derivative terms, and $\Lambda_{\omega}$ the would-be strong coupling energy scale.

In this paper we study the feathering substructures along spiral arms by considering the perturbational gas response to a spiral shock. Feathers are density fluctuations that jut out from the spiral arm to the inter-arm region at pitch angles given by the quantum numbers of the doubly-periodic structure. In a localized asymptotic approximation, related to the shearing sheet except that the inhomogeneities occur in space rather than in time, we derive the linearized perturbation equations for a razor-thin disk with turbulent interstellar gas, frozen-in magnetic field, and gaseous self-gravity. Apart from the modal quantum numbers, the individual normal modes of the system depend on seven dimensionless quantities that characterize the underlying time-independent axisymmetric state plus its steady, nonlinear, two-armed spiral-shock (TASS) response to a hypothesized background density-wave supported by the disk stars of the galaxy. We show that some of these normal modes have positive growth rates. Their over-density contours in the post-shock region are very reminiscent of observed feathering substructures in full magnetohydrodynamic (MHD) simulations. The feathering substructures are parasitic instabilities intrinsic to the system; thus, their study not only provides potential diagnostics for important parameters that characterize the interstellar medium of external galaxies, but also yields a deeper understanding of the basic mechanism that drives the formation of the giant molecular clouds (GMCs) and the OB stars that outline observed grand-design spirals.

We consider an extended theory of Horava-Lifshitz gravity with the detailed balance condition softly breaking, but without the projectability condition. With the former, the number of independent coupling constants is significantly reduced. With the latter and by extending the original foliation-preserving diffeomorphism symmetry $ {{Diff}}(M, {\cal{F}})$ to include a local U(1) symmetry, the spin-0 gravitons are eliminated. Thus, all the problems related to them disappear, including the instability, strong coupling, and different speeds in the gravitational sector. When the theory couples to a scalar field, we find that the scalar field is not only stable in both the ultraviolet (UV) and infrared (IR), but also free of the strong coupling problem, because of the presence of high-order spatial derivative terms of the scalar field. Furthermore, applying the theory to cosmology, we find that due to the additional U(1) symmetry, the Friedmann-Robertson-Walker (FRW) universe is necessarily flat. We also investigate the scalar, vector, and tensor perturbations of the flat FRW universe, and derive the general linearized field equations for each kind of the perturbations.

Young stars typically have strong magnetic fields, so that the magnetospheres of newly formed close binaries can interact, dissipate energy, and produce synchrotron radiation. The V773 Tau A binary system, a pair of T Tauri stars with a 51 day orbit, displays such a signature, with peak emission taking place near periastron. This paper proposes that the observed emission arises from the change in energy stored in the composite magnetic field of the system. We model the fields using the leading order (dipole) components and show that this picture is consistent with current observations. In this model, the observed radiation accounts for a fraction of the available energy of interaction between the magnetic fields from the two stars. Assuming antisymmetry, we compute the interaction energy $E_{\rm int}$ as a function of the stellar radii, the stellar magnetic field strengths, the binary semi-major axis, and orbital eccentricity, all of which can be measured independently of the synchrotron radiation. The variability in time and energetics of the synchrotron radiation depend on the details of the annihilation of magnetic fields through reconnection events, which generate electric fields that accelerate charged particles, and how those charged particles, especially fast electrons, are removed from the interaction region. However, the major qualitative features are well described by the background changes in the global magnetic configuration driven by the orbital motion. The theory can be tested by observing a collection of pre-main-sequence binary systems.

We apply the holographic principle and the equipartition law of energy to the apparent horizon of a Friedmann-Robertson-Walker universe and derive the Friedmann equation describing the dynamics of the universe. We also show that the equipartition law of energy can be interpreted as the first law of thermodynamics at the apparent horizon.

Previous analyses of magnetospheric accretion and outflow in classical T Tauri stars (CTTSs), within the context of both the X-wind model and other theoretical scenarios, have assumed a dipolar geometry for the stellar magnetic field if it were not perturbed by the presence of an accreting, electrically conducting disk. However, CTTS surveys reveal that accretion hot spots cover a small fraction of the stellar surface, and that the net field polarization on the stellar surface is small. Both facts imply that the magnetic field generated by the star has a complex non-dipolar structure. To address this discrepancy between theory and observations, we re-examine X-wind theory without the dipole constraint. Using simple physical arguments based on the concept of trapped flux, we show that a dipole configuration is in fact not essential. Independent of the precise geometry of the stellar magnetosphere, the requirement for a certain level of trapped flux predicts a definite relationship among various CTTS observables. Moreover, superposition of multipole stellar fields naturally yield small observed hot-spot covering fractions and small net surface polarizations. The generalized X-wind picture remains viable under these conditions, with the outflow from a small annulus near the inner disk edge little affected by the modified geometry, but with inflow highly dependent on the details of how the emergent stellar flux is linked and trapped by the inner disk regions. Our model is consistent with data, including recent spectropolarimetric measurements of the hot spot sizes and field strengths in V2129 Oph and BP Tau.

Strong magnetization makes the disks surrounding young stellar objects rotate at rates that are too sub-Keplerian to enable the thermal launching of disk winds from their surfaces unless the rate of gas diffusion across field lines is dynamically fast. This under-appreciated implication of disk magnetization poses a considerable challenge for disk-wind theory.

The interaction of accretion disks with the magnetospheres of young stars can produce X-winds and funnel flows. With the assumption of axial symmetry and steady state flow, the problem can be formulated in terms of quantities that are conserved along streamlines, such as the Bernoulli integral (BI), plus a partial differential equation (PDE), called the Grad-Shafranov equation (GSE), that governs the distribution of streamlines in the meridional plane. The GSE plus BI yields a PDE of mixed type, elliptic before critical surfaces where the flow speed equals certain characteristic wave speeds are crossed and hyperbolic afterward. The computational difficulties are exacerbated by the locations of the critical surfaces not being known in advance. To overcome these obstacles, we consider a variational principle by which the GSE can be attacked by extremizing an action integral, with all other conserved quantities of the problem explicitly included as part of the overall formulation. To simplify actual applications we adopt the cold limit of a negligibly small ratio of the sound speed to the speed of Keplerian rotation in the disk where the X-wind is launched. We also ignore the obstructing effects of any magnetic fields that might thread a disk approximated to be infinitesimally thin. We then introduce trial functions with adjustable coefficients to minimize the variations that give the GSE. We tabulate the resulting coefficients so that other workers can have analytic forms to reconstruct X-wind solutions for various astronomical, cosmochemical, and meteoritical applications.

This paper re-examines the problem of ambipolar diffusion as a mechanism for the production and runaway evolution of centrally condensed molecular cloud cores, a process that has been termed the gravomagneto catastrophe. Our calculation applies in the geometric limit of a highly flattened core and allows for a semi-analytic treatment of the full problem, although physical fixes are required to resolve a poor representation of the central region. A noteworthy feature of the overall formulation is that the solutions for the ambipolar diffusion portion of the evolution for negative times ($t < 0$) match smoothly onto the collapse solutions for positive times ($t > 0$). The treatment shows that the resulting cores display non-zero, but sub-magnetosonic, inward velocities at the end of the diffusion epoch, in agreement with current observations. Another important result is the derivation of an analytic relationship between the dimensionless mass to flux ratio $\lambda_0\equiv f_0^{-1}$ of the central regions produced by runaway core condensation and the dimensionless measure of the rate of ambipolar diffusion $\epsilon$. In conjunction with previous work showing that ambipolar diffusion takes place more quickly in the presence of turbulent fluctuations, i.e., that the effective value of $\epsilon$ can be enhanced by turbulence, the resultant theory provides a viable working hypothesis for the formation of isolated molecular-cloud cores and their subsequent collapse to form stars and planetary systems.

We consider the formation of OB associations from two perspectives: (a) the fractional gas consumption in star formation,epsilon, per dynamical time scale t_dyn in a galaxy, and (b) the origin of the so-called Kennicutt-Schmidt law that the rate of star formation per unit area is proportional to a power, alpha, of the surface density in HI and H_2 gas when certain thresholds are crossed. The empirical findings that epsilon is approximately 0.01 and alpha is approximately 1.4 or 1.5 have simple explanations if the rate of star formation is magnetically regulated. An empirical test of the ideas resides in an analysis of why giant OB associations are ``strung out like pearls along the arms" of spiral galaxies.

We consider the accretion process in a disk with magnetic fields that are dragged in from the interstellar medium by gravitational collapse. Two diffusive processes are at work in the system: (1) "viscous" torques exerted by turbulent and magnetic stresses, and (2) "resistive" redistribution of mass with respect to the magnetic flux arising from the imperfect conduction of current. In steady state, self-consistency between the two rates of drift requires that a relationship exists between the coefficients of turbulent viscosity and turbulent resistivity. Ignoring any interactions with a stellar magnetosphere, we solve the steady-state equations for a magnetized disk under the gravitational attraction of a mass point and threaded by an amount of magnetic flux consistent with calculations of magnetized gravitational collapse in star formation. Our model mean-field equations have an exact analytical solution that corresponds to magnetically diluted Keplerian rotation about the central mass point. The solution yields the strength of the magnetic field and the surface density as functions of radial position in the disk and their connection with the departure from pure Keplerian rotation in representative cases. We compare the predictions of the theory with the available observations concerning T Tauri stars, FU Orionis stars, and low- and high-mass protostars. Finally, we speculate on the physical causes for high and low states of the accretion disks that surround young stellar objects. One of the more important results of this study is the physical derivation of analytic expressions for the turbulent viscosity and turbulent resistivity.

We formulate the problem of magnetic field dissipation during the accretion phase of low-mass star formation, and we carry out the first step of an iterative solution procedure by assuming that the gas is in free-fall along radial field lines. This so-called ``kinematic approximation'' ignores the back reaction of the Lorentz force on the accretion flow. In quasi steady-state, and assuming the resistivity coefficient to be spatially uniform, the problem is analytically soluble in terms of Legendre's polynomials and confluent hypergeometric functions. The dissipation of the magnetic field occurs inside a region of radius inversely proportional to the mass of the central star (the ``Ohm radius''), where the magnetic field becomes asymptotically straight and uniform. In our solution, the magnetic flux problem of star formation is avoided because the magnetic flux dragged in the accreting protostar is always zero. Our results imply that the effective resistivity of the infalling gas must be higher by several orders of magnitude than the microscopic electric resistivity, to avoid conflict with measurements of paleomagnetism in meteorites and with the observed luminosity of regions of low-mass star formation.

We study the self-similar collapse of an isothermal magnetized rotating cloud in the ideal magnetohydrodynamic (MHD) regime. In the limit of small distance from the accreting protostar we find an analytic solution that corresponds to free-fall onto a central mass point. The density distribution is not spherically symmetric but depends on the mass loading of magnetic field lines, which can be obtained by matching our inner solution to an outer collapse solution previously computed by Allen, Shu & Li. The concentration of magnetic field trapped by the central mass point under field-freezing, independent on the details of the starting state, creates a split monopole configuration where the magnetic field strength increases as the inverse square of the distance from the center. Under such conditions, the inflow eventually becomes subalfvenic and the outward transfer of angular momentum by magnetic braking very efficient, thus preventing the formation of a centrifugally supported disk. Instead, the azimuthal velocity of the infalling gas decreases to zero at the center, and the gas spirals into the star. Therefore, the dissipation of dynamically important levels of magnetic field is a fundamental requisite for the formation of protoplanetary disks around young stars.

Two explanations exist for the short-lived radionuclides present in the solar system when the CAIs first formed. They originated either from the ejecta of a supernova or by the in situ irradiation of nebular dust by energetic particles. With a half-life of only 53 days, Beryllium-7 is then the key discriminant, since it can be made only by irradiation. We calculate the yield of Be-7. Within model uncertainties associated mainly with nuclear cross sections, we obtain agreement with the experimental value. Moreover, if Be-7 and Be-10 have the same origin, the irradiation time must be short. The x-wind model provides a natural astrophysical setting that gives the requisite conditions. The decoupling of the Al-26 and Be-10 observed in some rare CAIs receives a quantitative explanation when rare gradual events are considered. Finally, we show that the presence of supernova-produced Fe-60 in the solar accretion disk does not necessarily mean that other short-lived radionuclides have a stellar origin.

We perform a linear stability analysis of the axisymmetric, relativistic, self-similar isothermal disk against non-axisymmetric perturbations. Two sets of neutral modes are discovered. The first set corresponds to marginally unstable perturbations driven by gravitational radiation, and the other signals the onset of bifurcation to non-axisymmetric equilibrium solutions to the Einstein equations.

We investigate whether magnetic tension can define the masses of forming stars by holding up the subcritical envelope of a molecular cloud that suffers gravitational collapse of its supercritical core. We perform an equilibrium analysis of the initial and final states assuming perfect field freezing, no rotation, isothermality, and a completely flattened configuration. The sheet geometry allows us to separate the magnetic tension into a levitation associated with the split monopole formed by the trapped flux of the central star and a suspension associated with curved field lines that thread the static pseudodisk and envelope of material external to the star. We find solutions where the eigenvalue for the stellar mass is a fixed multiple of the initial core mass of the cloud. We verify the analytically derived result by an explicit numerical simulation of a closely related 3-D axisymmetric system. However, with field freezing, the implied surface magnetic fields much exceed measured values for young stars. If the pinch by the central split monopole were to be eliminated by magnetic reconnection, then magnetic suspension alone cannot keep the subcritical envelope (i.e., the entire model cloud) from falling onto the star. We argue that this answer has general validity, even if the initial state lacked any kind of symmetry, possessed rotation, and had a substantial level of turbulence. These findings strongly support a picture for the halt of infall that invokes dynamic levitation by YSO winds and jets, but the breakdown of ideal magnetohydrodynamics is required to allow the appearance in the problem of a rapidly rotating, centrifugally supported disk. We use these results to calculate the initial mass function and star formation efficiency for the distributed and clustered modes of star formation.

We study numerically the collapse of rotating, magnetized molecular cloud cores, focusing on rotation and magnetic braking during the main accretion phase of isolated star formation. Motivated by previous numerical work and analytic considerations, we idealize the pre-collapse core as a magnetized singular isothermal toroid, with a constant rotational speed everywhere. The collapse starts from the center, and propagates outwards in an inside-out fashion, satisfying exact self-similarity in space and time. For rotation rates and field strengths typical of dense low-mass cores, the main feature remains the flattening of the mass distribution along field lines -- the formation of a pseudodisk, as in the nonrotating cases. The density distribution of the pseudodisk is little affected by rotation. On the other hand, the rotation rate is strongly modified by pseudodisk formation. Most of the centrally accreted material reaches the vicinity of the protostar through the pseudodisk. The specific angular momentum can be greatly reduced on the way, by an order of magnitude or more, even when the pre-collapse field strength is substantially below the critical value for dominant cloud support. The efficient magnetic braking is due to the pinched geometry of the magnetic field in the pseudodisk, which strengthens the magnetic field and lengthens the level arm for braking. Both effects enhance the magnetic transport of angular momentum from inside to outside. The excess angular momentum is carried away in a low-speed outflow that has, despite claims made by other workers, little in common with observed bipolar molecular outflows. We discuss the implications of our calculations for the formation of true disks that are supported against gravity by rotation.

We study numerically the collapse of non-rotating, self-gravitating, magnetized, singular isothermal toroids characterized by sound speed, $a$, and level of magnetic to thermal support, $H_0$. In qualitative agreement with previous treatments by Galli & Shu and other workers, we find that the infalling material is deflected by the field lines towards the equatorial plane, creating a high-density, flattened structure -- a pseudodisk. The pseudodisk contracts dynamically in the radial direction, dragging the field lines threading it into a highly pinched configuration that resembles a split monopole. The oppositely directed field lines across the midplane and the large implied stresses may play a role in how magnetic flux is lost in the actual situation in the presence of finite resistivity or ambipolar diffusion. The infall rate into the central regions is given to 5% uncertainty by the formula, $\dot M = (1+H_0)a^3/G$, where $G$ is the universal gravitational constant, anticipated by semi-analytical studies of the self-similar gravitational collapses of the singular isothermal sphere and isopedically magnetized disks. The introduction of finite initial rotation results in a complex interplay between pseudodisk and true (Keplerian) disk formation that is examined in a companion paper.

We use hydrodynamical simulations to investigate the response of geometrically thin, self-gravitating, singular isothermal disks of gas to imposed rigidly rotating spiral potentials. By minimizing reflection-induced feedback from boundaries, and by restricting our attention to models where the swing parameter $X \sim 10$, we minimize the swing amplification of global normal modes even in models where Toomre's $Q_g \sim 1-2$ in the gas disk. We perform two classes of simulations: short-term ones over a few galactic revolutions where the background spiral forcing is large, and long-term ones over many galactic revolutions where the spiral forcing is considerably smaller. In both classes of simulations, the initial response of the gas disk is smooth and mimics the driving spiral field. At late times, many of the models evince substructure akin to the so-called branches, spurs, and feathers observed in real spiral galaxies. We comment on the parts played respectively by ultraharmonic resonances, reflection off internal features produced by nonlinear dredging, and local, transient, gravitational instabilites within spiral arms in the generation of such features. Our simulations reinforce the idea that spiral structure in the gaseous component becomes increasingly flocculent and disordered with the passage of time, even when the background population of old disk stars is a grand-design spiral. We speculate that truly chaotic behavior arises when many overlapping ultraharmonic resonances develop in reaction to an imposed spiral forcing that has itself a nonlinear, yet smooth, wave profile.

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

(Abridged) We consider the expansion of an initially self-gravitating, static, singular, isothermal cloud core. For t>0, the gas is ionized and heated to a higher uniform temperature by the formation of a luminous, but massless, star in its center. The approximation that the mass and gravity of the central star is negligible for the subsequent motion of the HII region holds for distances much greater than about 100 AU and for the massive cloud cores that give rise to high-mass stars. If the initial ionization and heating is approximated to occur instantaneously at t=0, then the subsequent flow (for r >> 100 AU) caused by the resulting imbalance between self-gravity and thermal pressure is self-similar. Because of the steep density profile, pressure gradients produce a shock front that travels into the cloud, accelerating the gas to supersonic velocities in what has been called the ``champagne phase.'' We then study the self-similar solutions of the expansion of HII regions embedded in molecular clouds characterized by more general power-law density distributions with exponent 3/2 < n < 3. In these cases, the shock velocity is an increasing function of the exponent n, and diverges as n tends to 3. We show that this happens because the model includes an origin, where the pressure driving the shock diverges because the enclosed heated mass is infinite. Our results imply that the continued photoevaporation of massive reservoirs of neutral gas (e.g., surrounding disks and/or globules) nearby to the embedded ionizing source is required in order to maintain over a significant timescale the emission measure observed in champagne flows.

We investigate the formation of planetesimals via the gravitational instability of solids that have settled to the midplane of a circumstellar disk. Vertical shear between the gas and a subdisk of solids induces turbulent mixing which inhibits gravitational instability. Working in the limit of small, well-coupled particles, we find that the mixing becomes ineffective when the surface density ratio of solids to gas exceeds a critical value. Solids in excess of this precipitation limit can undergo midplane gravitational instability and form planetesimals. However, this saturation effect typically requires increasing the local ratio of solid to gaseous surface density by factors of two to ten times cosmic abundances, depending on the exact properties of the gas disk. We discuss existing astrophysical mechanisms for augmenting the ratio of solids to gas in protoplanetary disks by such factors, and investigate a particular process that depends on the radial variations of orbital drift speeds induced by gas drag. This mechanism can concentrate millimeter sized chondrules to the supercritical surface density in < few x 10^6 years, a suggestive timescale for the disappearance of dusty disks in T Tauri stars. We discuss the relevance of our results to some outstanding puzzles in planet formation theory: the size of the observed solar system, and the rapid type I migration of Earth mass bodies.

Recent observations show that chromospheric activity in late-M and L dwarfs is much lower than in the earlier M types, in spite of comparatively rapid rotation. We investigate the possibility that this drop-off in activity results from the very high electrical resistivities in the dense, cool and predominantly neutral atmospheres of late-M and L dwarfs. We calculate magnetic field diffusivities in the atmospheres of objects with effective temperatures in the range 3000-1500 (mid-M to L), using the atmospheric structure models of Allard and Hauschildt. We find that the combination of very low ionization fraction and high density in these atmospheres results in very large resistivities due to neutral-charged particle collisions, and efficient field diffusion. The resistivities are found to increase with both decreasing optical depth, and decreasing effective temperature. As a result, any existing magnetic fields are increasingly decoupled from atmospheric motions as one moves from mid-M to L; we quantify this through a simple Reynolds number calculation. This, coupled with the difficulty in transporting magnetic stresses through the highly resistive atmosphere, can account for the observed drop in activity from mid-M to L, assuming activity in these objects is magnetically driven. We also examine the issue of acoustic heating, and find that this appears inadequate to explain the observed H-alpha fluxes in mid-M to L dwarfs. Consequently, magnetic heating does seem to be the most viable mechanism for generating activity in these objects. Finally, we speculate on a possible flare mechanism in these cool dwarfs.

We formulate and solve by semi-analytic means the axisymmetric equilibria of relativistic self-similar disks of infinitesimal vertical thickness. These disks are supported in the horizontal directions against their self-gravity by a combination of isothermal (two-dimensional) pressure and a flat rotation curve. The dragging of inertial frames restricts possible solutions to rotation speeds that are always less than 0.438 times the speed of light, a result first obtained by Lynden-Bell and Pineault in 1978 for a cold disk. We show that prograde circular orbits of massive test particles exist and are stable for all of our model disks, but retrograde circular orbits cannot be maintained with particle velocities less than the speed of light once the disk develops an ergoregion. We also compute photon trajectories, planar and non-planar, in the resulting spacetime, for disks with and without ergoregions. We find that all photon orbits, except for a set of measure zero, tend to be focused by the gravity of the flattened mass-energy distribution toward the plane of the disk. This result suggests that strongly relativistic, rapidly rotating, compact objects may have difficulty ejecting collimated beams of matter or light along the rotation axes until the flows get well beyond the flattened parts of the relativistic mass distribution (which cannot happen in the self-similar models considered in this paper).

In order to compare the x-wind with observations, one needs to be able to calculate its thermal and ionization properties. We formulate the physical basis for the streamline-by-streamline integration of the ionization and heat equations of the steady x-wind. In addition to the well-known processes associated with the interaction of stellar and accretion-funnel hot-spot radiation with the wind, we include X-ray heating and ionization, mechanical heating, and a revised calculation of ambipolar diffusion heating. The mechanical heating arises from fluctuations produced by star-disk interactions of the time dependent x-wind that are carried by the wind to large distances where they are dissipated in shocks, MHD waves, and turbulent cascades. We model the time-averaged heating by the scale-free volumetric heating rate, $\Gamma_{\rm mech} = \alpha \rho v^3 s^{-1}$, where $\rho$ and $v$ are the local mass density and wind speed, respectively, $s$ is the distance from the origin, and $\alpha$ is a phenomenological constant. When we consider a partially-revealed but active young stellar object, we find that choosing $\alpha \sim 10^{-3}$ in our numerical calculations produces temperatures and electron fractions that are high enough for the x-wind jet to radiate in the optical forbidden lines at the level and on the spatial scales that are observed. We also discuss a variety of applications of our thermal-chemical calculations that can lead to further observational checks of x-wind theory.

We review the difficulties of the classical fission and fragmentation hypotheses for the formation of binary and multiple stars. A crucial missing ingredient in previous theoretical studies is the inclusion of dynamically important levels of magnetic fields. As a minimal model for a candidate presursor to the formation of binary and multiple stars, we therefore formulate and solve the problem of the equilibria of isopedically magnetized, singular isothermal disks, without the assumption of axial symmetry. Considerable analytical progress can be made if we restrict our attention to models that are scale-free, i.e., that have surface densities that vary inversely with distance from the rotation axis of the system. In agreement with earlier analysis by Syer and Tremaine, we find that lopsided (M=1) configurations exist at any dimensionless rotation rate, including zero. Multiple-lobed (M = 2, 3, 4, ...) configurations bifurcate from an underlying axisymmetric sequence at progressively higher dimensionless rates of rotation, but such nonaxisymmetric sequences always terminate in shockwaves before they have a chance to fission into M=2, 3, 4, ... separate bodies. On the basis of our experience in this paper, and the preceding Paper I, we advance the hypothesis that binary and multiple star-formation from smooth (i.e., not highly turbulent) starting states that are supercritical but in unstable mechanical balance requires the rapid (i.e., dynamical) loss of magnetic flux at some stage of the ensuing gravitational collapse.

A crucial missing ingredient in previous theoretical studies of fragmentation is the inclusion of dynamically important levels of magnetic fields. As a minimal model for a candidate presursor to the formation of binary and multiple stars, we therefore consider the equilibrium configuration of isopedically magnetized, scale-free, singular isothermal disks, without the assumption of axial symmetry. We find that lopsided (M=1) configurations exist at any dimensionless rotation rate, including zero. Multiple-lobed (M = 2, 3, 4, ...) configurations bifurcate from an underlying axisymmetric sequence at progressively higher dimensionless rates of rotation, but such nonaxisymmetric sequences always terminate in shockwaves before they have a chance to fission into separate bodies. We advance the hypothesis that binary and multiple star-formation from smooth (i.e., not highly turbulent) starting states that are supercritical but in unstable mechanical balance requires the rapid (i.e., dynamical) loss of magnetic flux at some stage of the ensuing gravitational collapse.

We investigate the equilibrium properties of self-gravitating magnetized clouds with polytropic equations of state with negative index n. In particular, we consider scale-free isopedic configurations that have constant dimensionless spherical mass-to-flux ratio lambda_r and that may constitute ``pivotal'' states for subsequent dynamical collapse to form groups or clusters of stars. For given Gamma = 1 + 1/n, equilibria with smaller values of lambda_r are more flattened, ranging from spherical configurations with lambda_r = infinity to completely flattened states for lambda_r=1. For a given amount of support provided by the magnetic field as measured by the dimensionless parameter H_0, equilibria with smaller values of Gamma are more flattened. However, logatropic (defined by Gamma = 0) disks do not exist. The only possible scale-free isopedic equilibria with logatropic equation of state are spherical uniformly magnetized clouds.

We report high spectral resolution observations of the CO vibrational overtone emission from the young stellar object 1548C27; our observations include both the CO(2-0) and CO(5-3) bandhead regions. These data and similar observations of the young stellar object WL16, reported in a previous contribution to this journal, provide some of the most compelling evidence to date for the existence of inner disks around young stars. We describe the simple procedure that we use to synthesize bandhead emission from disks including the effect of thermal dissociation of CO and non-LTE excitation of the vibrational levels. Using this spectral synthesis procedure to extract the kinematics and physical properties of the emitting gas from the overtone data, we show how these high signal-to-noise data are also powerful probes of the stellar and inner disk properties of these systems. Our modeling is consistent with the identification of WL16 and 1548C27 as Herbig AeBe stars with stellar masses of approximately 2 and 4 solar masses, respectively. Thus, the kinematic signature of rotating disks in the overtone spectra of these sources provides strong support for the role of accretion disks in the formation of intermediate mass stars. For both WL16 and 1548C27, we interpret our modeling results as indicating that the overtone emission arises from a temperature inversion region in the inner disk atmosphere. We also find evidence for suprathermally broadened lines and are able to place useful constraints on the radial temperature and column density distributions of the CO line-formation region of the disk atmosphere.