Consider a dynamic general relativistic spacetime in which the proper infinitesimal interval along one spatial coordinate direction decreases monotonically with time, while the corresponding intervals increase along other spatial directions. In a system undergoing such complete anisotropic collapse/expansion, we look for the formation of a cosmic double-jet configuration: free test particles in the ambient medium, relative to the collapsing system, gain energy from the gravitational field and asymptotically line up parallel and antiparallel to the direction of collapse such that their Lorentz factors approach infinity. A strong burst of electromagnetic radiation is expected to accompany this event if some of the free test particles carry electric charge. Previous work in this direction involved mainly Ricci-flat spacetimes; hence, we concentrate here on inhomogeneous perfect fluid spacetimes. We briefly explore the possible connection between these theoretical cosmic jets and astrophysical jets. We also discuss other general relativistic scenarios for the formation of cosmic jets.

We present extensions of the treatment contained in our recent paper on nonlocal Newtonian cosmology [C. Chicone and B. Mashhoon, J. Math. Phys. 57, 072501 (2016)]. That is, the implications of the recent nonlocal generalization of Einstein's theory of gravitation are further investigated within the regime of Newtonian cosmology. In particular, we treat the nonlocal problem of structure formation for a spherically symmetric expanding dust model and show numerically that as the central density contrast grows, it tends to decrease slowly with radial distance as the universe expands. The nonlocal violation of Newton's shell theorem provides a physical interpretation of our numerical results.

We explore some of the cosmological implications of the recent classical nonlocal generalization of Einstein's theory of gravitation in which nonlocality is due to the gravitational memory of past events. In the Newtonian regime of this theory, the nonlocal character of gravity simulates dark matter in spiral galaxies and clusters of galaxies. However, dark matter is considered indispensable as well for structure formation in standard models of cosmology. Can nonlocal gravity solve the problem of structure formation without recourse to dark matter? Here we make a beginning in this direction by extending nonlocal gravity in the Newtonian regime to the cosmological domain. The nonlocal analog of the Zel'dovich solution is formulated and the consequences of the resulting nonlocal Zel'dovich model are investigated in detail.

The implications of the recent classical nonlocal generalization of Einstein's theory of gravitation for gravitational physics in the Solar System are investigated. In this theory, the nonlocal character of gravity appears to simulate dark matter. Nonlocal gravity in the Newtonian regime involves a reciprocal kernel with three spatial parameters, of which two have already been determined from the rotation curves of spiral galaxies and the internal dynamics of clusters of galaxies. However, the short-range parameter a_0 remains to be determined. In this connection, the nonlocal contribution to the perihelion precession of a planetary orbit is estimated and a preliminary lower limit on a_0 is determined.

We investigate gravitational radiation in the linear approximation within the framework of the recent nonlocal generalization of Einstein's theory of gravitation. In this theory, nonlocality can simulate dark matter; in fact, in the Newtonian regime, we recover the phenomenological Tohline-Kuhn approach to modified gravity. To account for the observational data regarding the rotation curves of spiral galaxies, nonlocality is associated with a characteristic length scale of order \lambda_0 = 10 kpc. It follows that in nonlocal gravity, the treatment of extremely low-frequency (~ 10^{-12} Hz) gravitational waves with wavelengths of order \lambda_0 would be quite different than in general relativity. However, for radiation of frequency > 10^{-8} Hz, which is the frequency range that is the focus of current observational searches, the corresponding wavelengths are very small compared to \lambda_0. We find that in this frequency regime the nonlocal deviations from general relativity essentially average out and can be safely neglected in practice.

To confront relativity theory with observation, it is necessary to split spacetime into its temporal and spatial components. The (1+3) timelike threading approach involves restrictions on the gravitational potentials $(g_{\mu \nu})$, while the (3+1) spacelike slicing approach involves restrictions on $(g^{\mu \nu})$. These latter coordinate conditions protect chronology within any such coordinate patch. While the threading coordinate conditions can be naturally integrated into the structure of Lorentzian geometry and constitute the standard coordinate conditions in general relativity, this circumstance does not extend to the slicing coordinate conditions. We explore the influence of chronology violation on wave motion. In particular, we consider the propagation of radiation parallel to the rotation axis of stationary G\"odel-type universes characterized by parameters $\eta > 0$ and $\lambda > 0$ such that for $\eta < 1$ ($\eta >1$) chronology is protected (violated). We show that in the WKB approximation such waves can freely propagate only when chronology is protected.

The recent nonlocal generalization of Einstein's theory of gravitation reduces in the Newtonian regime to a nonlocal and nonlinear modification of Poisson's equation of Newtonian gravity. The nonlocally modified Poisson equation implies that nonlocality can simulate dark matter. Observational data regarding dark matter provide limited information about the functional form of the reciprocal kernel, from which the original nonlocal kernel of the theory must be determined. We study this inverse problem of nonlocal gravity in the linear domain, where the applicability of the Fourier transform method is critically examined and the conditions for the existence of the nonlocal kernel are discussed. This approach is illustrated via simple explicit examples for which the kernels are numerically evaluated. We then turn to a general discussion of the modified Poisson equation and present a formal solution of this equation via a successive approximation scheme. The treatment is specialized to the gravitational potential of a point mass, where in the linear regime we recover the Tohline-Kuhn approach to modified gravity.

In dynamic spacetimes in which asymmetric gravitational collapse/expansion is taking place, the timelike geodesic equation appears to exhibit an interesting property: Relative to the collapsing configuration, free test particles undergo gravitational "acceleration" and form a double-jet configuration parallel to the axis of collapse. We illustrate this aspect of peculiar motion in simple spatially homogeneous cosmological models such as the Kasner spacetime. To estimate the effect of spatial inhomogeneities on cosmic jets, timelike geodesics in the Ricci-flat double-Kasner spacetime are studied in detail. While spatial inhomogeneities can significantly modify the structure of cosmic jets, we find that under favorable conditions the double-jet pattern can initially persist over a finite period of time for sufficiently small inhomogeneities.

We present a family of dynamic rotating cylindrically symmetric Ricci-flat gravitational fields whose geodesic motions have the structure of gravitomagnetic jets. These correspond to helical motions of free test particles up and down parallel to the axis of cylindrical symmetry and are reminiscent of the motion of test charges in a magnetic field. The speed of a test particle in a gravitomagnetic jet asymptotically approaches the speed of light. Moreover, numerical evidence suggests that jets are attractors. The possible implications of our results for the role of gravitomagnetism in the formation of astrophysical jets are briefly discussed.

We discuss time-dependent gravitational fields that "accelerate" free test particles to the speed of light resulting in cosmic double-jet configurations. It turns out that complete gravitational collapse along a spatial axis together with corresponding expansion along the other two axes leads to the accelerated motion of free test particles up and down parallel to the collapse axis such that a double-jet pattern is asymptotically formed with respect to the collapsed configuration.

We study a simple class of time-dependent rotating Ricci-flat cylindrically symmetric spacetime manifolds whose geodesics admit gravitomagnetic jets. The helical paths of free test particles in these jets up and down parallel to the rotation axis are analogous to those of charged particles in a magnetic field. The jets are attractors. The jet speed asymptotically approaches the speed of light. In effect, such source-free spacetime regions act as "gravitomagnetic accelerators".

The Newtonian regime of a recent nonlocal extension of general relativity (GR) is investigated. Nonlocality is introduced via a scalar "constitutive" kernel in a special case of the translational gauge theory of gravitation, namely, the teleparallel equivalent of GR. In this theory, the nonlocal aspect of gravity simulates dark matter. A nonlocal and nonlinear generalization of Poisson's equation of Newtonian gravitation is presented. The implications of nonlocality for the gravitational physics in the solar system are briefly studied.

We study the linear post-Newtonian approximation to general relativity known as gravitoelectromagnetism (GEM); in particular, we examine the similarities and differences between GEM and electrodynamics. Notwithstanding some significant differences between them, we find that a special nonstationary metric in GEM can be employed to show {\it explicitly} that it is possible to introduce gravitational induction within GEM in close analogy with Faraday's law of induction and Lenz's law in electrodynamics. Some of the physical implications of gravitational induction are briefly discussed.

We consider bound geodesic orbits of test masses in the exterior gravitational field of a rotating astronomical source whose proper angular momentum varies linearly with time. The linear perturbation approach of Lense and Thirring is herein extended to the nonstationary case. In particular, we investigate the instability of Lense-Thirring precessing orbits due to the slow temporal variation of the gravitomagnetic field of the source.

Fermi coordinates are directly constructed in de Sitter and Goedel spacetimes and the corresponding exact coordinate transformations are given explicitly. The quasi-inertial Fermi coordinates are then employed to discuss the dynamics of a free test particle in these spacetimes and the results are compared to the corresponding generalized Jacobi equations that contain only the lowest-order tidal terms. The domain of validity of the generalized Jacobi equation is thus examined in these cases. Furthermore, the difficulty of constructing explicit Fermi coordinates in black-hole spacetimes is demonstrated.

The motion of free nearby test particles relative to a stable equatorial circular geodesic orbit about a Kerr source is investigated. It is shown that the nonlinear generalized Jacobi equation can be transformed in this case to an autonomous form. Tidal dynamics beyond the critical speed c/sqrt(2) is studied. We show, in particular, that a free test particle vertically launched from the circular orbit parallel or antiparallel to the Kerr rotation axis is tidally accelerated if its initial relative speed exceeds c/sqrt(2). Possible applications of our results to high-energy astrophysics are briefly mentioned.

We investigate the motion of free relativistic particles relative to the ambient medium around a gravitationally collapsed system. If the relative speed exceeds a critical value given by c/sqrt(2), the gravitational tidal effects exhibit novel features that are contrary to Newtonian expectations. In particular, ultrarelativistic jet clumps moving freely outward along the rotation axis strongly decelerate with respect to the ambient medium, while ultrarelativistic particles strongly accelerate in directions normal to the jet axis. The implications of these direct consequences of general relativity for jets in microquasars and the origin of the high-energy cosmic rays are briefly mentioned.

The relative motion of a classical relativistic spinning test particle is studied with respect to a nearby free test particle in the gravitational field of a rotating source. The effects of the spin-curvature coupling force are elucidated and the implications of the results for the motion of rotating plasma clumps in astrophysical jets are discussed.

Imagine a swarm of free particles near a point P outside a gravitating mass M and a free reference particle at P that is on a radial escape trajectory away from M. Relative to this reference particle and in a Fermi normal coordinate system constructed along its worldline, the particles in the swarm that move along the radial direction and are ultrarelativistic (that is, they have speeds above c/sqrt(2)) decelerate toward this terminal speed. On the other hand, the swarm particles that are ultrarelativistic and move in directions normal to the radial (jet) directions accelerate to almost the speed of light by the gravitational tidal force of the mass M. The implications of these effects as well as the influence of the higher-order terms on the tidal acceleration mechanism are investigated. The observational evidence in support of these general relativistic effects is briefly discussed.

We point out novel consequences of general relativity involving tidal dynamics of ultrarelativistic relative motion. Specifically, we use the generalized Jacobi equation and its extension to study the force-free dynamics of relativistic flows near a massive rotating source. We show that along the rotation axis of the gravitational source, relativistic tidal effects strongly decelerate an initially ultrarelativistic flow with respect to the ambient medium, contrary to Newtonian expectations. Moreover, an initially ultrarelativistic flow perpendicular to the axis of rotation is strongly accelerated by the relativistic tidal forces. The astrophysical implications of these results for jets and ultrahigh energy cosmic rays are briefly mentioned.

Fermi coordinates are the natural generalization of inertial Cartesian coordinates to accelerated systems and gravitational fields. We study the motion of ultrarelativistic particles and light rays in Fermi coordinates and investigate inertial and tidal effects beyond the critical speed c/sqrt(2). In particular, we discuss the black-hole tidal acceleration mechanism for ultrarelativistic particles in connection with a possible origin for high-energy cosmic rays.

In the description of \emph{relative} motion in accelerated systems and gravitational fields, inertial and tidal accelerations must be taken into account, respectively. These involve a critical speed that in the first approximation can be simply illustrated in the case of motion in one dimension. For one-dimensional motion, such first-order accelerations are multiplied by $(1-V^2/V_c^2)$, where $V_c=c/\sqrt{2}$ is the critical speed. If the speed of relative motion exceeds $V_c$, there is a sign reversal with consequences that are contrary to Newtonian expectations.

Dynamics of relativistic outflows along the rotation axis of a Kerr black hole is investigated using a simple model that takes into account the relativistic tidal force of the central source as well as the Lorentz force due to the large-scale electromagnetic field which is assumed to be present in the ambient medium. The evolution of the speed of the flow relative to the ambient medium is studied. In the force-free case, the resulting equation of motion predicts rapid deceleration of the initial flow and an asymptotic relative speed with a Lorentz factor of 2^1/2. In the presence of the Lorentz force, the long-term relative speed of the clump tends to the ambient electrical drift speed.

The Jacobi equation in pseudo-Riemannian geometry determines the linearized geodesic flow. The linearization ignores the relative velocity of the geodesics. The generalized Jacobi equation takes the relative velocity into account; that is, when the geodesics are neighboring but their relative velocity is arbitrary the corresponding geodesic deviation equation is the generalized Jacobi equation. The Hamiltonian structure of this nonlinear equation is analyzed in this paper. The tidal accelerations for test particles in the field of a plane gravitational wave and the exterior field of a rotating mass are investigated. In the latter case, the existence of an attractor of uniform relative radial motion with speed $2^{-1/2}c\approx 0.7 c$ is pointed out. The astrophysical implications of this result for the terminal speed of a relativistic jet is briefly explored.

We define the general Hill system and briefly analyze its dynamical behavior. A particular Hill system representing the interaction of a Keplerian binary system with a normally incident circularly polarized gravitational wave is discussed in detail. In this case, we compute the Poincar\'e-Melnikov function explicitly and determine its zeros. Moreover, we provide numerical evidence in favor of chaos in this system. The partially averaged equations for the Hill system are used to predict the regular behavior of the Keplerian orbit at resonance with the external radiation.

The general phenomena associated with sustained resonance are studied in this paper in connection with relativistic binary pulsars. We represent such a system by two point masses in a Keplerian binary system that evolves via gravitational radiation damping as well as an external tidal perturbation. For further simplification, we assume that the external tidal perturbation is caused by a normally incident circularly polarized monochromatic gravitational wave. In this case, the second-order partially averaged equations are studied and a theorem of C. Robinson is employed to prove that for certain values of the physical parameters resonance capture followed by sustained resonance is possible in the averaged system. We conjecture that sustained resonance can occur in the physical system when the perturbing influences nearly balance each other.

The long-term dynamical evolution of a Keplerian binary orbit due to the emission and absorption of gravitational radiation is investigated. This work extends our previous results on transient chaos in the planar case to the three dimensional Kepler system. Specifically, we consider the nonlinear evolution of the relative orbit due to gravitational radiation damping as well as external gravitational radiation that is obliquely incident on the initial orbital plane. The variation of orbital inclination, especially during resonance capture, turns out to be very sensitive to the initial conditions. Moreover, we discuss the novel phenomenon of chaotic transition.