Nonlinear Gravitational Waves: Their Form and Effects

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The nonlinear essence of gravitational waves

A critical review of gravitational wave theory is made. It is pointed out that the usual linear approach to the gravitational wave theory is neither conceptually consistent nor mathematically justified. Relying upon that analysis it is argued that-analogously to a Yang-Mills propagating field, which must be nonlinear to carry its gauge charge-a gravitational wave must necessarily be nonlinear to transport its own charge-that is, energy-momentum.

Singularities and horizons in the collisions of gravitational waves

This thesis presents a study of the dynamical, nonlinear interaction of colliding gravitational waves, as described by classical general relativity. It is focused mainly on two fundamental questions: First, what is the general structure of the singularities and Killing-Cauchy horizons produced in the collisions of exactly plane-symmetric gravitational waves? Second, under what conditions will the collisions of almost-plane gravitational waves (waves with large but finite transverse sizes) produce singularities?

In the work on the collisions of exactly-plane waves, it is shown that Killing horizons in any plane-symmetric spacetime are unstable against small plane-symmetric perturbations. It is thus concluded that the Killing-Cauchy horizons produced by the collisions of some exactly plane gravitational waves are nongeneric, and that generic initial data for the colliding plane waves always produce "pure" spacetime singularities without such horizons. This conclusion is later proved rigorously (using the full nonlinear theory rather than perturbation theory), in connection with an analysis of the asymptotic singularity structure of a general colliding plane-wave spacetime. This analysis also proves that asymptotically the singularities created by colliding plane waves are of inhomogeneous-Kasner type; the asymptotic Kasner axes and exponents of these singularities in general depend on the spatial coordinate that runs tangentially to the singularity in the non-plane-symmetric direction.

In the work on collisions of almost-plane gravitational waves, first some general properties of single almost-plane gravitational-wave spacetimes are explored. It is shown that, by contrast with an exact plane wave, an almost-plane gravitational wave cannot have a propagation direction that is Killing; i.e., it must diffract and disperse as it propagates. It is also shown that an almost-plane wave cannot be precisely sandwiched between two null wavefronts; i.e., it must leave behind tails in the spacetime region through which it passes. Next, the occurrence of spacetime singularities in the collisions of almost-plane waves is investigated. It is proved that if two colliding, almost-plane gravitational waves are initially exactly plane-symmetric across a central region of sufficiently large but finite transverse dimensions, then their collision produces a spacetime singularity with the same local structure as in the exact-plane-wave collision. Finally, it is shown that a singularity still forms when the central regions are only approximately plane-symmetric initially. Stated more precisely, it is proved that if the colliding almost-plane waves are initially sufficiently close to being exactly plane-symmetric across a bounded central region of sufficiently large transverse dimensions, then their collision necessarily produces spacetime singularities. In this case, nothing is now known about the local and global structures of the singularities.

Evolution and decay of cylindrical and spherical nonlinear acoustic waves generated by a sinusoidal source

The present work gives a comprehensive numerical study of the evolution and decay of cylindrical and spherical nonlinear acoustic waves generated by a sinusoidal source. Using pseudospectral and predictor–corrector implicit finite difference methods, we first reproduced the known analytic results of the plane harmonic problem to a high degree of accuracy. The non-planar harmonic problems, for which the amplitude decay is faster than that for the planar case, are then treated. The results are correlated with the known asymptotic results of Scott (1981) and Enflo (1985). The constant in the old-age formula for the cylindrical canonical problem is found to be 1.85 which is rather close to 2, ‘estimated’ analytically by Enflo. The old-age solutions exhibiting strict symmetry about the maximum are recovered; these provide an excellent analytic check on the numerical solutions. The evolution of the waves for different source geometries is depicted graphically.

Instability Analysis of Nonlinear Surface Waves in a Circular Cylindrical Container Subjected to a Vertical Excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, the fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates damping term and external excitation, was derived from linearized Navier-Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated. (c) 2005 Elsevier SAS. All rights reserved.

Nonlinear Faraday Waves in a Parametrically Excited Circular Cylindrical Container

In the cylindrical coordinate system, a singular perturbation theory of multiple-scale asymptotic expansions was developed to study single standing water wave mode by solving potential equations of water waves in a rigid circular cylinder, which is subjec

An Efficient Front-Tracking Method For Fully Nonlinear Interfacial Waves

A fully nonlinear and dispersive model within the framework of potential theory is developed for interfacial (2-layer) waves. To circumvent the difficulties arisen from the moving boundary problem a viable technique based on the mixed Eulerian and Lagrangian concept is proposed: the computing area is partitioned by a moving mesh system which adjusts its location vertically to conform to the shape of the moving boundaries but keeps frozen in the horizontal direction. Accordingly, a modified dynamic condition is required to properly compute the boundary potentials. To demonstrate the effectiveness of the current method, two important problems for the interfacial wave dynamics, the generation and evolution processes, are investigated. Firstly, analytical solutions for the interfacial wave generations by the interaction between the barotropic tide and topography are derived and compared favorably with the numerical results. Furthermore simulations are performed for the nonlinear interfacial wave evolutions at various water depth ratios and satisfactory agreement is achieved with the existing asymptotical theories. (c) 2008 Elsevier Inc. All rights reserved.

Topics in linear and nonlinear dispersive waves

The various singularities and instabilities which arise in the modulation theory of dispersive wavetrains are studied. Primary interest is in the theory of nonlinear waves, but a study of associated questions in linear theory provides background information and is of independent interest.

The full modulation theory is developed in general terms. In the first approximation for slow modulations, the modulation equations are solved. In both the linear and nonlinear theories, singularities and regions of multivalued modulations are predicted. Higher order effects are considered to evaluate this first order theory. An improved approximation is presented which gives the true behavior in the singular regions. For the linear case, the end result can be interpreted as the overlap of elementary wavetrains. In the nonlinear case, it is found that a sufficiently strong nonlinearity prevents this overlap. Transition zones with a predictable structure replace the singular regions.

For linear problems, exact solutions are found by Fourier integrals and other superposition techniques. These show the true behavior when breaking modulations are predicted.

A numerical study is made for the anharmonic lattice to assess the nonlinear theory. This confirms the theoretical predictions of nonlinear group velocities, group splitting, and wavetrain instability, as well as higher order effects in the singular regions.

Nonlinear dispersive waves in nonlinear optics

A study is made of solutions of the macroscopic Maxwell equations in nonlinear media. Both nonlinear and dispersive terms are responsible for effects that are not taken into account in the geometrical optics approximation. The nonlinear terms can, depending on the nature of the nonlinearity, cause plane waves to focus when the amplitude varies across the wavefront. The dispersive terms prevent the singularities that nonlinearity alone would produce. Solutions are found which de scribe periodic plane waves in fully nonlinear media. Equations describing the evolution of the amplitude, frequency and wave number are generated by means of averaged Lagrangian techniques. The equations are solved for near linear media to produce the form of focusing waves which develop a singularity at the focal point. When higher dispersion is included nonlinear and dispersive effects can balance and one finds amplitude profiles that propagate with straight rays.

Part I. On the breaking of nonlinear dispersive waves. Part II. Variational principles in continuum mechanics

A model equation for water waves has been suggested by Whitham to study, qualitatively at least, the different kinds of breaking. This is an integro-differential equation which combines a typical nonlinear convection term with an integral for the dispersive effects and is of independent mathematical interest. For an approximate kernel of the form e^(-b|x|) it is shown first that solitary waves have a maximum height with sharp crests and secondly that waves which are sufficiently asymmetric break into "bores." The second part applies to a wide class of bounded kernels, but the kernel giving the correct dispersion effects of water waves has a square root singularity and the present argument does not go through. Nevertheless the possibility of the two kinds of breaking in such integro-differential equations is demonstrated.

Difficulties arise in finding variational principles for continuum mechanics problems in the Eulerian (field) description. The reason is found to be that continuum equations in the original field variables lack a mathematical "self-adjointness" property which is necessary for Euler equations. This is a feature of the Eulerian description and occurs in non-dissipative problems which have variational principles for their Lagrangian description. To overcome this difficulty a "potential representation" approach is used which consists of transforming to new (Eulerian) variables whose equations are self-adjoint. The transformations to the velocity potential or stream function in fluids or the scaler and vector potentials in electromagnetism often lead to variational principles in this way. As yet no general procedure is available for finding suitable transformations. Existing variational principles for the inviscid fluid equations in the Eulerian description are reviewed and some ideas on the form of the appropriate transformations and Lagrangians for fluid problems are obtained. These ideas are developed in a series of examples which include finding variational principles for Rossby waves and for the internal waves of a stratified fluid.

I. Singular perturbation problems involving singular points and turning points. II. On the averaged Lagrangian technique for nonlinear dispersive waves

In Part I a class of linear boundary value problems is considered which is a simple model of boundary layer theory. The effect of zeros and singularities of the coefficients of the equations at the point where the boundary layer occurs is considered. The usual boundary layer techniques are still applicable in some cases and are used to derive uniform asymptotic expansions. In other cases it is shown that the inner and outer expansions do not overlap due to the presence of a turning point outside the boundary layer. The region near the turning point is described by a two-variable expansion. In these cases a related initial value problem is solved and then used to show formally that for the boundary value problem either a solution exists, except for a discrete set of eigenvalues, whose asymptotic behaviour is found, or the solution is non-unique. A proof is given of the validity of the two-variable expansion; in a special case this proof also demonstrates the validity of the inner and outer expansions.

Nonlinear dispersive wave equations which are governed by variational principles are considered in Part II. It is shown that the averaged Lagrangian variational principle is in fact exact. This result is used to construct perturbation schemes to enable higher order terms in the equations for the slowly varying quantities to be calculated. A simple scheme applicable to linear or near-linear equations is first derived. The specific form of the first order correction terms is derived for several examples. The stability of constant solutions to these equations is considered and it is shown that the correction terms lead to the instability cut-off found by Benjamin. A general stability criterion is given which explicitly demonstrates the conditions under which this cut-off occurs. The corrected set of equations are nonlinear dispersive equations and their stationary solutions are investigated. A more sophisticated scheme is developed for fully nonlinear equations by using an extension of the Hamiltonian formalism recently introduced by Whitham. Finally the averaged Lagrangian technique is extended to treat slowly varying multiply-periodic solutions. The adiabatic invariants for a separable mechanical system are derived by this method.

The importance of spin for observing gravitational waves from coalescing compact binaries with LIGO and Virgo

General Relativity predicts the existence of gravitational waves, which carry information about the physical and dynamical properties of their source. One of the many promising sources of gravitational waves observable by ground-based instruments, such as in LIGO and Virgo, is the coalescence of two compact objects (neutron star or black hole). Black holes and neutron stars sometimes form binaries with short orbital periods, radiating so strongly in gravitational waves that they coalesce on astrophysically short timescales. General Relativity gives precise predictions for the form of the signal emitted by these systems. The most recent searches for theses events used waveform models that neglected the effects of black hole and neutron star spin. However, real astrophysical compact objects, especially black holes, are expected to have large spins. We demonstrate here a data analysis infrastructure which achieves an improved sensitivity to spinning compact binaries by the inclusion of spin effects in the template waveforms. This infrastructure is designed for scalable, low-latency data analysis, ideal for rapid electromagnetic followup of gravitational wave events.

The search for gravitational waves from the coalescence of black hole binary systems in data from the LIGO and Virgo detectors. Or: A dark walk through a random forest

The LIGO and Virgo gravitational-wave observatories are complex and extremely sensitive strain detectors that can be used to search for a wide variety of gravitational waves from astrophysical and cosmological sources. In this thesis, I motivate the search for the gravitational wave signals from coalescing black hole binary systems with total mass between 25 and 100 solar masses. The mechanisms for formation of such systems are not well-understood, and we do not have many observational constraints on the parameters that guide the formation scenarios. Detection of gravitational waves from such systems — or, in the absence of detection, the tightening of upper limits on the rate of such coalescences — will provide valuable information that can inform the astrophysics of the formation of these systems. I review the search for these systems and place upper limits on the rate of black hole binary coalescences with total mass between 25 and 100 solar masses. I then show how the sensitivity of this search can be improved by up to 40% by the the application of the multivariate statistical classifier known as a random forest of bagged decision trees to more effectively discriminate between signal and non-Gaussian instrumental noise. I also discuss the use of this classifier in the search for the ringdown signal from the merger of two black holes with total mass between 50 and 450 solar masses and present upper limits. I also apply multivariate statistical classifiers to the problem of quantifying the non-Gaussianity of LIGO data. Despite these improvements, no gravitational-wave signals have been detected in LIGO data so far. However, the use of multivariate statistical classification can significantly improve the sensitivity of the Advanced LIGO detectors to such signals.

Gauge Invariant Spectral Cauchy Characteristic Extraction of Gravitational Waves in Computational General Relativity

We present a complete system for Spectral Cauchy characteristic extraction (Spectral CCE). Implemented in C++ within the Spectral Einstein Code (SpEC), the method employs numerous innovative algorithms to efficiently calculate the Bondi strain, news, and flux.

Spectral CCE was envisioned to ensure physically accurate gravitational wave-forms computed for the Laser Interferometer Gravitational wave Observatory (LIGO) and similar experiments, while working toward a template bank with more than a thousand waveforms to span the binary black hole (BBH) problem’s seven-dimensional parameter space.

The Bondi strain, news, and flux are physical quantities central to efforts to understand and detect astrophysical gravitational wave sources within the Simulations of eXtreme Spacetime (SXS) collaboration, with the ultimate aim of providing the first strong field probe of the Einstein field equation.

In a series of included papers, we demonstrate stability, convergence, and gauge invariance. We also demonstrate agreement between Spectral CCE and the legacy Pitt null code, while achieving a factor of 200 improvement in computational efficiency.

Spectral CCE represents a significant computational advance. It is the foundation upon which further capability will be built, specifically enabling the complete calculation of junk-free, gauge-free, and physically valid waveform data on the fly within SpEC.

Quantum Emulation of Gravitational Waves

Gravitational waves, as predicted by Einstein's general relativity theory, appear as ripples in the fabric of spacetime traveling at the speed of light. We prove that the propagation of small amplitude gravitational waves in a curved spacetime is equivalent to the propagation of a subspace of electromagnetic states. We use this result to propose the use of entangled photons to emulate the evolution of gravitational waves in curved spacetimes by means of experimental electromagnetic setups featuring metamaterials.