The first neutron star (NS) merger observed by advanced LIGO and Virgo, GW170817, and its fireworks of electromagnetic counterparts across the entire electromagnetic spectrum marked the beginning of multi-messenger astronomy and astrophysics with gravitational waves. The ultraviolet, optical, and near-infrared emission was consistent with being powered by the radioactive decay of nuclei synthesized in the merger ejecta by the rapid neutron capture process (r-process). Starting from an outline of the inferred properties of this 'kilonova' emission, I discuss possible astrophysical sites for r-process nucleosynthesis in NS mergers, arguing that the heaviest r-process elements synthesized in this event most likely originated in outflows from a post-merger accretion disk. I compare the inferred properties of r-process element production in GW170817 to current observational constraints on galactic heavy r-process nucleosynthesis and discuss challenges merger-only models face in explaining the r-process content of our galaxy. Based on the observational properties of GW170817 and recent theoretical progress on r-process nucleosynthesis in collapsars, I then show how GW170817 points to collapsars as the dominant source of r-process enrichment in the Milky Way. These rare core-collapse events arguably better satisfy existing constraints and overcome problems related to r-process enrichment in various environments that NS mergers face. Finally, I comment on the universality of the r-process and on how variations in light r-process elements can be obtained both in NS mergers and collapsars.

Star-to-star dispersion of r-process elements has been observed in a significant number of old, metal-poor globular clusters. We investigate early-time neutron-star mergers as the mechanism for this enrichment. Through both numerical modeling and analytical arguments, we show that neutron-star mergers cannot be induced through dynamical interactions early in the history of the cluster, even when the most liberal assumptions about neutron-star segregation are assumed. Therefore, if neutron-star mergers are the primary mechanism for r-process dispersion in globular clusters, they likely result from the evolution of isolated, primordial binaries in the clusters. Through population modeling, we find that moderate fractions of GCs with enrichment are only possible when a significant number of double neutron-star progenitors proceed through Case BB mass transfer --- under various assumptions for the initial properties of globular clusters, a neutron-star merger with the potential for enrichment will occur in ~15-60% (~30-90%) of globular clusters if this mass transfer proceeds stably (unstably). The strong anti-correlation between the pre-supernova orbital separation and post-supernova systemic velocity due to mass loss in the supernova leads to efficient ejection of most enrichment candidates from their host clusters. Thus, most enrichment events occur shortly after the double neutron stars are born. This requires star-forming gas that can absorb the r-process ejecta to be present in the globular cluster 30-50 Myr after the initial burst of star formation. If scenarios for redistributing gas in globular clusters cannot act on these timescales, the number of neutron-star merger enrichment candidates drops severely, and it is likely that another mechanism, such as r-process enrichment from collapsars, is at play.

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

The binary neutron star merger GW170817 detected by Advanced LIGO/Virgo was followed by a short gamma-ray burst (GRB) and thermal radiation ('kilonova') powered by the radioactive decay of heavy nuclei synthesized in the merger ejecta by the rapid neutron capture process ('r-process'). The large inferred quantity of ejecta is best understood as originating in an outflow from the accretion disk surrounding the newly-formed black hole, the same engine that was likely responsible for the GRB jet. Similar accretion flows accompany the collapse of rotating massive stars ('collapsars'), powering the class of long GRBs and their associated supernovae. Here we show that collapsar accretion disks also produce neutron-rich outflows that synthesize heavy r-process nuclei, despite the comparatively proton-rich composition of the infalling star. Though occurring less frequently than mergers, the much greater accreted mass in collapsars---and their correspondingly larger disk wind ejecta---implicate them as dominant contributors to the Galactic r-process. Collapsars provide the rare r-process source needed to promptly enrich a small fraction of ultra-faint dwarf galaxies early in the Universe which is also compatible with the Galactic chemical enrichment of europium relative to iron over longer timescales. We predict an excess in the late-time near-infrared emission from GRB supernovae, testable by future observations.

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event.

Merging binaries consisting of two neutron stars (NSs) or an NS and a stellar-mass black hole typically form a massive accretion torus around the remnant black hole or long-lived NS. Outflows from these neutrino-cooled accretion disks represent an important site for $r$-process nucleosynthesis and the generation of kilonovae. We present the first three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulations including weak interactions and a realistic equation of state of such accretion disks over viscous timescales ($380\,\mathrm{ms}$). We witness the emergence of steady-state MHD turbulence, a magnetic dynamo with an $\sim\!20\,\mathrm{ms}$ cycle, and the generation of a `hot' disk corona that launches powerful thermal outflows aided by the energy released as free nucleons recombine into $\alpha$-particles. We identify a self-regulation mechanism that keeps the midplane electron fraction low ($Y_e\sim0.1$) over viscous timescales. This neutron-rich reservoir, in turn, feeds outflows that retain a sufficiently low value of $Y_e\approx 0.2$ to robustly synthesize third-peak $r$-process elements. The quasi-spherical outflows are projected to unbind $40\%$ of the initial disk mass with typical asymptotic escape velocities of $0.1c$, and may thus represent the dominant mass ejection mechanism in NS-NS mergers. Including neutrino absorption, our findings agree with previous hydrodynamical $\alpha-$disk simulations that the entire range of $r$-process nuclei from the first to the third $r$-process peak can be synthesized in the outflows, in good agreement with observed solar system abundances. The asymptotic escape velocities and the quantity of ejecta, when extrapolated to moderately higher disk masses, are consistent with those needed to explain the red kilonova emission following the NS merger GW170817.

General-relativistic magnetohydrodynamic (GRMHD) simulations are an important tool to study a variety of astrophysical systems such as neutron star mergers, core-collapse supernovae, and accretion onto compact objects. A conservative GRMHD scheme numerically evolves a set of conservation equations for 'conserved' quantities and requires the computation of certain primitive variables at every time step. This recovery procedure constitutes a core part of any conservative GRMHD scheme and it is closely tied to the equation of state (EOS) of the fluid. In the quest to include nuclear physics, weak interactions, and neutrino physics, state-of-the-art GRMHD simulations employ finite-temperature, composition-dependent EOSs. While different schemes have individually been proposed, the recovery problem still remains a major source of error, failure, and inefficiency in GRMHD simulations with advanced microphysics. The strengths and weaknesses of the different schemes when compared to each other remain unclear. Here we present the first systematic comparison of various recovery schemes used in different dynamical spacetime GRMHD codes for both analytic and tabulated microphysical EOSs. We assess the schemes in terms of (i) speed, (ii) accuracy, and (iii) robustness. We find large variations among the different schemes and that there is not a single ideal scheme. While the computationally most efficient schemes are less robust, the most robust schemes are computationally less efficient. More robust schemes may require an order of magnitude more calls to the EOS, which are computationally expensive. We propose an optimal strategy of an efficient three-dimensional Newton-Raphson scheme and a slower but more robust one-dimensional scheme as a fall-back.

This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core collapse supernovae are making rapid progress. Detection of both, electron neutrinos and antineutrinos from the next galactic supernova will constrain the composition of neutrino-driven winds and provide unique nucleosynthesis information. Finally FRIB and other rare-isotope beam facilities will soon have dramatic new capabilities to synthesize many neutron-rich nuclei that are involved in the r-process. The new capabilities can significantly improve our understanding of the r-process and likely resolve one of the main outstanding problems in classical nuclear astrophysics. However, to make best use of the new experimental capabilities and to fully interpret the results, a great deal of infrastructure is needed in many related areas of astrophysics, astronomy, and nuclear theory. We will place these experiments in context by discussing astrophysical simulations and observations of r-process sites, observations of stellar abundances, galactic chemical evolution, and nuclear theory for the structure and reactions of very neutron-rich nuclei. This review paper was initiated at a three-week International Collaborations in Nuclear Theory program in June 2016 where we explored promising r-process experiments and discussed their likely impact, and their astrophysical, astronomical, and nuclear theory context.

The merger of binary neutron stars, or of a neutron star and a stellar-mass black hole, can result in the formation of a massive rotating torus around a spinning black hole. In addition to providing collimating media for gamma-ray burst jets, unbound outflows from these disks are an important source of mass ejection and rapid neutron capture (r-process) nucleosynthesis. We present the first three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations of neutrino-cooled accretion disks in neutron star mergers, including a realistic equation of state valid at low densities and temperatures, self-consistent evolution of the electron fraction, and neutrino cooling through an approximate leakage scheme. After initial magnetic field amplification by magnetic winding, we witness the vigorous onset of turbulence driven by the magneto-rotational instability (MRI). The disk quickly reaches a balance between heating from MRI-driven turbulence and neutrino cooling, which regulates the midplane electron fraction to a low equilibrium value $Y_\text{e} \approx 0.1$. Over the 380 ms duration of the simulation, we find that a fraction $\approx 20\%$ of the initial torus mass is unbound in powerful outflows with velocities $v \approx 0.03-0.1\,c$ and electron fractions $Y_\text{e} \approx 0.1-0.25$. Post-processing the outflows through a nuclear reaction network shows the production of a robust second and third peak r-process. Though broadly consistent with the results of previous axisymmetric hydrodynamical simulations, extrapolation of our results to late times suggests that the total ejecta mass from GRMHD disks is significantly higher. Our results provide strong evidence that post-merger disk outflows are an important site for the r-process.

Merging binary neutron stars (BNSs) represent the ultimate targets for multimessenger astronomy, being among the most promising sources of gravitational waves (GWs), and, at the same time, likely accompanied by a variety of electromagnetic counterparts across the entire spectrum, possibly including short gamma-ray bursts (SGRBs) and kilonova/macronova transients. Numerical relativity simulations play a central role in the study of these events. In particular, given the importance of magnetic fields, various aspects of this investigation require general relativistic magnetohydrodynamics (GRMHD). So far, most GRMHD simulations focused the attention on BNS mergers leading to the formation of a hypermassive NS, which, in turn, collapses within few tens of ms into a black hole surrounded by an accretion disk. However, recent observations suggest that a significant fraction of these systems could form a long-lived NS remnant, which will either collapse on much longer timescales or remain indefinitely stable. Despite the profound implications for the evolution and the emission properties of the system, a detailed investigation of this alternative evolution channel is still missing. Here, we follow this direction and present a first detailed GRMHD study of BNS mergers forming a long-lived NS. We consider magnetized binaries with different mass ratios and equations of state and analyze the structure of the NS remnants, the rotation profiles, the accretion disks, the evolution and amplification of magnetic fields, and the ejection of matter. Moreover, we discuss the connection with the central engine of SGRBs and provide order-of-magnitude estimates for the kilonova/macronova signal. Finally, we study the GW emission, with particular attention to the post-merger phase.

Binary neutron star (BNS) mergers are the leading model to explain the phenomenology of short gamma-ray bursts (SGRBs), which are among the most luminous explosions in the universe. Recent observations of long-lasting X-ray afterglows of SGRBs challenge standard paradigms and indicate that in a large fraction of events a long-lived neutron star (NS) may be formed rather than a black hole. Understanding the mechanisms underlying these afterglows is necessary in order to address the open questions concerning the nature of SGRB central engines. However, recent theoretical progress has been hampered by the fact that the timescales of interest for the afterglow emission are inaccessible to numerical relativity simulations. Here we present a detailed model to bridge the gap between numerical simulations of the merger process and the relevant timescales for the afterglows, assuming that the merger results in a long-lived NS. This model is formulated in terms of a set of coupled differential equations that follow the evolution of the post-merger system and predict its electromagnetic (EM) emission in a self-consistent way, starting from initial data that can be extracted from BNS merger simulations and taking into account the most relevant radiative processes. Moreover, the model can accomodate the collapse of the remnant NS at any time during the evolution as well as different scenarios for the prompt SGRB emission. A second major reason of interest for BNS mergers is that they are considered the most promising source of gravitational waves (GWs) for detection with the advanced ground-based detector network LIGO/Virgo coming online this year. Multimessenger astronomy with joint EM and GW observations of the merger and post-merger phase can greatly enhance the scientific output of either type of observation. However, the actual benefit depends on ...

Recent observations indicate that in a large fraction of binary neutron star (BNS) mergers a long-lived neutron star (NS) may be formed rather than a black hole. Unambiguous electromagnetic (EM) signatures of such a scenario would strongly impact our knowledge on how short gamma-ray bursts (SGRBs) and their afterglow radiation are generated. Furthermore, such EM signals would have profound implications for multimessenger astronomy with joint EM and gravitational-wave (GW) observations of BNS mergers, which will soon become reality with the ground-based advanced LIGO/Virgo GW detector network starting its first science run this year. Here we explore such EM signatures based on the model presented in a companion paper, which provides a self-consistent evolution of the post-merger system and its EM emission starting from an early baryonic wind phase and resulting in a final pulsar wind nebula that is confined by the previously ejected material. Lightcurves and spectra are computed for a wide range of post-merger physical properties and particular attention is paid to the emission in the X-ray band. In the context of SGRB afterglow modeling, we present X-ray lightcurves corresponding to the 'standard' and the recently proposed 'time-reversal' scenario (SGRB prompt emission produced at the time of merger or at the time of collapse of the long-lived NS). The resulting afterglow lightcurve morphologies include, in particular, single and two-plateau features with timescales and luminosities that are in good agreement with the observations by the Swift satellite. Furthermore, we compute the X-ray signal that should precede the SGRB in the time-reversal scenario. If found, such a signal would represent smoking-gun evidence for this scenario. Finally, we find a bright, highly isotropic EM transient signal peaking in the X-ray band ...

Recent observations by the Swift satellite have revealed long-lasting ($\sim 10^2-10^5\,\mathrm{s}$), "plateau-like" X-ray afterglows in the vast majority of short gamma-ray bursts events. This has put forward the idea of a long-lived millisecond magnetar central engine being generated in a binary neutron star (BNS) merger and being responsible for the sustained energy injection over these timescales ("magnetar model"). We elaborate here on recent simulations that investigate the early evolution of such a merger remnant in general-relativistic magnetohydrodynamics. These simulations reveal very different conditions than those usually assumed for dipole spin-down emission in the magnetar model. In particular, the surrounding of the newly formed NS is polluted by baryons due to a dense, highly magnetized and isotropic wind from the stellar surface that is induced by magnetic field amplification in the interior of the star. The timescales and luminosities of this wind are compatible with early X-ray afterglows, such as the "extended emission". These isotropic winds are a generic feature of BNS merger remnants and thus represent an attractive alternative to current models of early X-ray afterglows. Further implications to BNS mergers and short gamma-ray bursts are discussed.

After decades of observations the physical mechanisms that generate short gamma-ray bursts (SGRBs) still remain unclear. Observational evidence provides support to the idea that SGRBs originate from the merger of compact binaries, consisting of two neutron stars (NSs) or a NS and a black hole (BH). Theoretical models and numerical simulations seem to converge to an explanation in which the central engine of SGRBs is given by a spinning BH surrounded by a hot accretion torus. Such a BH-torus system can be formed in compact binary mergers and is able to launch a relativistic jet, which can then produce the SGRB. This basic scenario, however, has recently been challenged by Swift satellite observations, which have revealed long-lasting X-ray afterglows in association with a large fraction of SGRB events. The long durations of these afterglows (from minutes to several hours) cannot be explained by the $\sim\text{s}$ accretion timescale of the torus onto the BH, and, instead, suggest a long-lived NS as the persistent source of radiation. Yet, if the merger results in a massive NS the conditions to generate a relativistic jet and thus the prompt SGRB emission are hardly met. Here we consider an alternative scenario that can reconcile the two aspects and account for both the prompt and the X-ray afterglow emission. Implications for future observations, multi-messenger astronomy and for constraining NS properties are discussed, as well as potential challenges for the model.

Short gamma-ray bursts (SGRBs) are among the most luminous explosions in the Universe and their origin still remains uncertain. Observational evidence favors the association with binary neutron star or neutron star-black hole (NS-BH) binary mergers. Leading models relate SGRBs to a relativistic jet launched by the BH-torus system resulting from the merger. However, recent observations have revealed a large fraction of SGRB events accompanied by X-ray afterglows with durations $\sim10^2\!-\!10^5~\mathrm{s}$, suggesting continuous energy injection from a long-lived central engine, which is incompatible with the short ($\lesssim1~\mathrm{s}$) accretion timescale of a BH-torus system. The formation of a supramassive NS, resisting the collapse on much longer spin-down timescales, can explain these afterglow durations, but leaves serious doubts on whether a relativistic jet can be launched at merger. Here we present a novel scenario accommodating both aspects, where the SGRB is produced after the collapse of a supramassive NS. Early differential rotation and subsequent spin-down emission generate an optically thick environment around the NS consisting of a photon-pair nebula and an outer shell of baryon-loaded ejecta. While the jet easily drills through this environment, spin-down radiation diffuses outwards on much longer timescales and accumulates a delay that allows the SGRB to be observed before (part of) the long-lasting X-ray signal. By analyzing diffusion timescales for a wide range of physical parameters, we find delays that can generally reach $\sim10^5~\mathrm{s}$, compatible with observations. The success of this fundamental test makes this "time-reversal" scenario an attractive alternative to current SGRB models.

Besides being among the most promising sources of gravitational waves, merging neutron star binaries also represent a leading scenario to explain the phenomenology of short gamma-ray bursts (SGRBs). Recent observations have revealed a large subclass of SGRBs with roughly constant luminosity in their X-ray afterglows, lasting $10\!-\!10^4$ s. These features are generally taken as evidence of a long-lived central engine powered by the magnetic spin-down of a uniformly rotating, magnetized object. We propose a different scenario in which the central engine powering the X-ray emission is a differentially rotating hypermassive neutron star (HMNS) that launches a quasi-isotropic and baryon-loaded wind driven by the magnetic field, which is built-up through differential rotation. Our model is supported by long-term, three-dimensional, general-relativistic, and ideal magnetohydrodynamic simulations, showing that this isotropic emission is a very robust feature. For a given HMNS, the presence of a collimated component depends sensitively on the initial magnetic field geometry, while the stationary electromagnetic luminosity depends only on the magnetic energy initially stored in the system. We show that our model is compatible with the observed timescales and luminosities and express the latter in terms of a simple scaling relation.

The universe is expected to be permeated by a stochastic background of gravitational radiation of astrophysical and cosmological origin. This background is capable of exciting oscillations in solar-like stars. Here we show that solar-like oscillators can be employed as giant hydrodynamical detectors for such a background in the muHz to mHz frequency range, which has remained essentially unexplored until today. We demonstrate this approach by using high-precision radial velocity data for the Sun to constrain the normalized energy density of the stochastic gravitational-wave background around 0.11 mHz. These results open up the possibility for asteroseismic missions like CoRoT and Kepler to probe fundamental physics.

We present a hydrodynamic model that describes excitation of linear stellar oscillations by a stochastic background of gravitational waves (SBGW) of astrophysical and cosmological origin. We find that this excitation mechanism is capable of generating solar g-mode amplitudes close to or comparable with values expected from excitation by turbulent convection, which is considered to be the main driving force for solar-like oscillations. A method is presented that places direct upper bounds on the SBGW in a frequency range, in which the SBGW is expected to contain rich astrophysical information. Employing estimates for solar g-mode amplitudes, the proposed method is demonstrated to have the potential to compete with sensitivities reached by gravitational wave experiments in other frequency ranges.

Mergers of binary neutron stars likely lead to the formation of a hypermassive neutron star (HMNS), which is metastable and eventually collapses to a black hole. This merger scenario is thought to explain the phenomenology of short gamma-ray bursts (SGRBs). The very high energies observed in SGRBs have been suggested to stem from neutrino-antineutrino annihilation and/or from very strong magnetic fields created during or after the merger by mechanisms like the magnetorotational instability (MRI). Here, we report on results that show for the first time the development of the MRI in HMNSs in three-dimensional, fully general-relativistic magnetohydrodynamic simulations. This instability amplifies magnetic fields exponentially and could be a vital ingredient in solving the SGRB puzzle.

A differentially rotating hypermassive neutron star (HMNS) is a metastable object which can be formed in the merger of neutron-star binaries. The eventual collapse of the HMNS into a black hole is a key element in generating the physical conditions expected to accompany the launch of a short gamma-ray burst. We investigate the influence of magnetic fields on HMNSs by performing three-dimensional simulations in general-relativistic magnetohydrodynamics. In particular, we provide direct evidence for the occurrence of the magnetorotational instability (MRI) in HMNS interiors. For the first time in simulations of these systems, rapidly-growing and spatially-periodic structures are observed to form with features like those of the channel flows produced by the MRI in other systems. Moreover, the growth time and wavelength of the fastest-growing mode are extracted and compared successfully with analytical predictions. The MRI emerges as an important mechanism to amplify magnetic fields over the lifetime of the HMNS, whose collapse to a black hole is accelerated. The evidence provided here that the MRI can actually develop in HMNSs could have a profound impact on the outcome of the merger of neutron-star binaries and on its connection to short gamma-ray bursts.

Starting from a general relativistic framework a hydrodynamic formalism is derived that yields the mean-square amplitudes and rms surface velocities of normal modes of non-relativistic stars excited by arbitrary gravitational wave (GW) radiation. In particular, stationary GW fields are considered and the resulting formulae are evaluated for two general types of GW radiation: radiation from a particular astrophysical source (e.g., a binary system) and a stochastic background of gravitational waves (SBGW). Expected sources and signal strengths for both types of GW radiation are reviewed and discussed. Numerical results for the Sun show that low-order quadrupolar g modes are excited more strongly than p modes by orders of magnitude. Maximal rms surface velocities in the case of excitation by astrophysical sources are found to be v {\le} 10^(-8) mm/s, assuming GW strain amplitudes of h {\le} 10^(-20). It is shown that current models for an SBGW produced by cosmic strings, with Omega_GW ~ 10^(-8)-10^(-5) in the frequency range of solar g modes, are able to produce maximal solar g-mode rms surface velocities of 10^(-5)-10^(-3) mm/s. This result lies close to or within the amplitude range of 10^(-3)-1 mm/s expected from excitation by turbulent convection, which is currently considered to be responsible for stellar g-mode excitation. It is concluded that studying g-mode observations of stars other than the Sun, in which excitation by GWs could be even more effective due to different stellar structures, might provide a new method to either detect GWs or to deduce a significant direct upper limit on an SBGW at intermediate frequencies between the pulsar bound and the bounds from interferometric detectors on Earth.

The excitation of solar and solar-like g modes in non-relativistic stars by arbitrary external gravitational wave fields is studied starting from the full field equations of general relativity. We develop a formalism that yields the mean-square amplitudes and surface velocities of global normal modes excited in such a way. The isotropic elastic sphere model of a star is adopted to demonstrate this formalism and for calculative simplicity. It is shown that gravitational waves solely couple to quadrupolar spheroidal eigenmodes and that normal modes are only sensitive to the spherical component of the gravitational waves having the same azimuthal order. The mean-square amplitudes in case of stationary external gravitational waves are given by a simple expression, a product of a factor depending on the resonant properties of the star and the power spectral density of the gravitational waves' spherical accelerations. Both mean-square amplitudes and surface velocities show a characteristic R^8-dependence (effective R^2-dependence) on the radius of the star. This finding increases the relevance of this excitation mechanism in case of stars larger than the Sun.