Showing 1–47 of 47 results for author: Novak, J

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Addressing the problem of the LIGO-Virgo-KAGRA visibility in the scientific literature

Authors: Pablo Barneo, Giuseppe Cabras, Pierre-Francois Cohadon, Livia Conti, Davide Guerra, Edoardo Milotti, Jerome Novak, Agata Trovato, Andrea Virtuoso

Abstract: As members of the Virgo Collaboration -- one of the large scientific collaborations that explore the universe of gravitational waves together with the LIGO Scientific Collaboration and the KAGRA Collaboration -- we became aware of biased citation practices that exclude Virgo, as well as KAGRA, from achievements that collectively belong to the wider LIGO/Virgo/KAGRA Collaboration. Here, we frame th… ▽ More As members of the Virgo Collaboration -- one of the large scientific collaborations that explore the universe of gravitational waves together with the LIGO Scientific Collaboration and the KAGRA Collaboration -- we became aware of biased citation practices that exclude Virgo, as well as KAGRA, from achievements that collectively belong to the wider LIGO/Virgo/KAGRA Collaboration. Here, we frame these practices in the context of Merton's Matthew effect, extending the reach of this well studied cognitive bias to include large international scientific collaborations. We provide qualitative evidence of its occurrence, displaying the network of links among published papers in the scientific literature related to Gravitational Wave science. We note how the keyword LIGO is linked to a much larger number of papers and variety of subjects than the keyword Virgo. We support these qualitative observations with a quantitative study based on a year-long monitoring of the relevant literature, where we scan all new preprints appearing in the arXiv electronic preprint database. Over the course of one year we identified all preprints failing to assign due credits to Virgo. As a further step, we undertook positive actions by asking the authors of problematic papers to correct them. We also report on a more in-depth investigation which we performed on problematic preprints that appeared in the first three months of the period under consideration, checking how frequently their authors reacted positively to our request and corrected their papers. Finally, we measure the global impact of papers classified as problematic and observe that, thanks to the changes implemented in response to our requests, the global impact (measured as the number of citations of papers which still contain Virgo visibility issues) was halved. We conclude the paper with general considerations for future work in a wider perspective. △ Less

Submitted 29 January, 2024; originally announced February 2024.

Comments: 17 pages, 9 figures, supplementary material (data and software) in linked public repository

Journal ref: The European Physical Journal H, (2024) 49: 2

New lower background and higher rate technique for anti-neutrino detection using Tungsten 183 Isotope

Authors: Jarred Novak, Nickolas Solomey, Brooks Hartsock, Brian Doty, Jonathan Folkerts

Abstract: Low energy anti-neutrinos detected from reactors or other sources have typically used the conversion of an anti-neutrino on Hydrogen, producing a positron and a free neutron. This neutron is subsequently captured on a secondary element with a large neutron capture cross section such as gadolinium or cadmium. We have studied the anti-neutrino conversion and suggest other elements that have a compar… ▽ More Low energy anti-neutrinos detected from reactors or other sources have typically used the conversion of an anti-neutrino on Hydrogen, producing a positron and a free neutron. This neutron is subsequently captured on a secondary element with a large neutron capture cross section such as gadolinium or cadmium. We have studied the anti-neutrino conversion and suggest other elements that have a comparable cross section for anti-neutrino reactions. With most neutron captures on gadolinium, it is possible to get two or three delayed gamma signals of known energy to occur. Experiments like ATLAS can make measurements with timing on the order of 25 ns. With electronics like these, this leads to the possibility of having a triple delayed coincidence using the positron annihilation on atomic shell electrons as the starting signal. We have also found an isotope of tungsten, $^{183}$W that offers a large anti-neutrino interaction cross section of $1.19 \times 10^{-46}$ m$^2$ and an anti-neutrino threshold energy for the production of the ground state at 2.094 MeV and the production of the first excited state at 2.167 MeV. This reaction makes a nuclear m1 excited state of $^{183}$Ta$^*$ that emits a signature secondary gamma pulse of 73 keV with a 106 ns half-life. This offers a new delayed coincidence technique that can be used to identify anti-neutrinos with lower background noise. This allows for less shielding than required for modern inverse beta decay detectors. △ Less

Submitted 17 June, 2024; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: Provisional Patent submitted, submitted to JINST for publication, revised version answers all referee questions, final version as accepted

One- and two-argument equation of state parametrizations with continuous sound speed for neutron star simulations

Authors: Gaël Servignat, Philip J. Davis, Jérôme Novak, Micaela Oertel, José A. Pons

Abstract: We describe two fitting schemes that aim to represent the high-density part of realistic equations of state for numerical simulations such as neutron star oscillations. The low-density part of the equation of state is represented by an arbitrary polytropic crust, and we propose a generic procedure to stitch any desired crust to the high-density fit, which is performed on the internal energy, press… ▽ More We describe two fitting schemes that aim to represent the high-density part of realistic equations of state for numerical simulations such as neutron star oscillations. The low-density part of the equation of state is represented by an arbitrary polytropic crust, and we propose a generic procedure to stitch any desired crust to the high-density fit, which is performed on the internal energy, pressure and sound speed of barotropic equations of state that describe cold neutron stars in $β$-equilibrium. An extension of the fitting schemes to equations of state with an additional compositional argument is proposed. In particular we develop a formalism that ensures the existence of a $β$-equilibrium at low densities. An additional feature of this low-density model is that it can be, in principle, applied to any parametrization. The performance of the fits is checked on mass, radius and tidal deformability as well as on the dynamical radial oscillation frequencies. To that end, we use a pseudospectral isolated neutron star evolution code based on a non-conservative form of the hydrodynamical equations. A comparison to existing parametrizations is proposed, as far as possible, and to published radial frequency values in the literature. The static and dynamic quantities are well reproduced by the fitting schemes. Our results suggest that, even though the radius is very sensitive to the choice of the crust, this choice has little influence on the oscillation frequencies of a neutron star. △ Less

Submitted 17 May, 2024; v1 submitted 5 November, 2023; originally announced November 2023.

Comments: 21 pages, 12 figures. This version matches the accepted version, with typos corrected

Waveform Modelling for the Laser Interferometer Space Antenna

Authors: LISA Consortium Waveform Working Group, Niayesh Afshordi, Sarp Akçay, Pau Amaro Seoane, Andrea Antonelli, Josu C. Aurrekoetxea, Leor Barack, Enrico Barausse, Robert Benkel, Laura Bernard, Sebastiano Bernuzzi, Emanuele Berti, Matteo Bonetti, Béatrice Bonga, Gabriele Bozzola, Richard Brito, Alessandra Buonanno, Alejandro Cárdenas-Avendaño, Marc Casals, David F. Chernoff, Alvin J. K. Chua, Katy Clough, Marta Colleoni, Mekhi Dhesi, Adrien Druart , et al. (121 additional authors not shown)

Abstract: LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmologic… ▽ More LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome. △ Less

Submitted 20 December, 2023; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Detectability of a phase transition in neutron star matter with third-generation gravitational wave interferometers

Authors: C. Mondal, M. Antonelli, F. Gulminelli, M. Mancini, J. Novak, M. Oertel

Abstract: Possible strong first-order hadron-quark phase transitions in neutron star interiors leave an imprint on gravitational waves, which could be detected with planned third-generation interferometers. Given a signal from the late inspiral of a binary neutron star (BNS) coalescence, %the possibility of assessing the presence of such a phase transition depends on the precision that can be attained in th… ▽ More Possible strong first-order hadron-quark phase transitions in neutron star interiors leave an imprint on gravitational waves, which could be detected with planned third-generation interferometers. Given a signal from the late inspiral of a binary neutron star (BNS) coalescence, %the possibility of assessing the presence of such a phase transition depends on the precision that can be attained in the determination of the tidal deformability parameter, as well as on the model used to describe the hybrid star equation of state. For the latter, we employ here a phenomenological meta-modelling of the equation of state that largely spans the parameter space associated with both the low density phase and the quark high density compatible with current constraints. We show that with a network of third-generation detectors, a single loud BNS event might be sufficient to infer the presence of a phase transition at low baryon densities with an average Bayes factor $B\approx 100$, up to a luminosity distance ($\mathcal{D}_L \lesssim$ 300 Mpc). △ Less

Submitted 11 May, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Journal ref: Mon.Not.Roy.Astron.Soc. 524 (2023) 3, 3464-3473

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO's and Advanced Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1680 additional authors not shown)

Abstract: We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate t… ▽ More We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO's and Advanced Virgo's third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim 0.6$ (at 90% confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions we are unable to rule out $f_\mathrm{PBH} = 1$. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound $f_{\mathrm{DBH}} < 10^{-5}$ on the fraction of atomic dark matter collapsed into black holes. △ Less

Submitted 26 January, 2024; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: https://dcc.ligo.org/P2200139

A Black-Hole Excision Scheme for General Relativistic Core-Collapse Supernova Simulations

Authors: B. Sykes, B. Mueller, I. Cordero-Carrión, P. Cerdá-Durán, J. Novak

Abstract: Fallback supernovae and the collapsar scenario for long-gamma ray burst and hypernovae have received considerable interest as pathways to black-hole formation and extreme transient events. Consistent simulations of these scenarios require a general relativistic treatment and need to deal appropriately with the formation of a singularity. Free evolution schemes for the Einstein equations can handle… ▽ More Fallback supernovae and the collapsar scenario for long-gamma ray burst and hypernovae have received considerable interest as pathways to black-hole formation and extreme transient events. Consistent simulations of these scenarios require a general relativistic treatment and need to deal appropriately with the formation of a singularity. Free evolution schemes for the Einstein equations can handle the formation of black holes by means of excision or puncture schemes. However, in constrained schemes, which offer distinct advantages in long-term numerical stability in stellar collapse simulations over well above $10^{4}$ light-crossing time scales, the dynamical treatment of black-hole spacetimes is more challenging. Building on previous work on excision in conformally flat spacetimes, we here present the implementation of a black-hole excision scheme for supernova simulations with the CoCoNuT-FMT neutrino transport code. We describe in detail a choice of boundary conditions that ensures long-time numerical stability, and also address upgrades to the hydrodynamics solver that are required to stably evolve the relativistic accretion flow onto the black hole. The scheme is currently limited to a spherically symmetric metric, but the hydrodynamics can be treated multi-dimensionally. For demonstration, we present a spherically symmetric simulation of black-hole formation in an $85 M_\odot$ star, as well as a two-dimensional simulation of the fallback explosion of the same progenitor. These extend past 9s and 0.3s after black-hole formation, respectively. △ Less

Submitted 22 May, 2024; v1 submitted 23 October, 2022; originally announced October 2022.

Comments: 19 pages, 12 figures. Accepted by PRD. Updated to correct typographical error in funding acknowledgements

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to bala… ▽ More We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more. △ Less

Submitted 2 January, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, Open Access Journal PDF

Report number: LIGO-P2100110-v13

Journal ref: The Astrophysical Journal Letters, 941, L30 (2022)

Standard cooling of rapidly rotating isolated neutron stars in 2D

Authors: Mikhail V. Beznogov, Jérôme Novak, Dany Page, Adriana R. Raduta

Abstract: We study the thermal evolution of axisymmetric rotating neutron stars in full general relativity. To this aim we develop "NSCool 2D Rot", a major upgrade of the 1D neutron stars thermal evolution code "NSCool" by D. Page. As a first application of our new code we address the standard cooling of isolated neutron stars with rotation frequencies up to the mass shedding limit. We investigate the effec… ▽ More We study the thermal evolution of axisymmetric rotating neutron stars in full general relativity. To this aim we develop "NSCool 2D Rot", a major upgrade of the 1D neutron stars thermal evolution code "NSCool" by D. Page. As a first application of our new code we address the standard cooling of isolated neutron stars with rotation frequencies up to the mass shedding limit. We investigate the effects of the equation of state (EOS) by considering different combinations of core and crust EOSs. The results indicate complex time-dependent evolution of temperature distribution throughout the whole volume of the star and, in particular, in the crust. We show that most of that complexity can be attributed to the formation of a "heat blob" in the crust and to the latitude dependence of the heat diffusion timescale through the crust. △ Less

Submitted 9 December, 2022; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: 18 pages, 9 figures, 2 tables; ApJ, in press

Journal ref: ApJ, 942, 72 (2023)

Concept for a Space-based Near-Solar Neutrino Detector

Authors: N. Solomey, J. Folkerts, H. Meyer, C. Gimar, J. Novak, B. Doty, T. English, L. Buchele, A. Nelsen, R. McTaggart, M. Christl

Abstract: The concept of putting a neutrino detector in close orbit of the sun has been unexplored until very recently. The primary scientific return is to vastly enhance our understanding of the solar interior, which is a major NASA goal. Preliminary calculations show that such a spacecraft, if properly shielded, can operate in space environments while taking data from neutrino interactions. These interact… ▽ More The concept of putting a neutrino detector in close orbit of the sun has been unexplored until very recently. The primary scientific return is to vastly enhance our understanding of the solar interior, which is a major NASA goal. Preliminary calculations show that such a spacecraft, if properly shielded, can operate in space environments while taking data from neutrino interactions. These interactions can be distinguished from random background rates of solar electromagnetic emissions, galactic charged cosmic-rays, and gamma-rays by using a double pulsed signature. Early simulations of this project have shown this veto schema to be successful in eliminating background and identifying the neutrino interaction signal in upwards of 75% of gamma ray interactions and nearly 100% of other interactions. Hence, we propose a new instrument to explore and study our sun. Due to inverse square scaling, this instrument has the potential to outperform earth-based experiments in several domains such as making measurements not accessible from the earth's orbit. △ Less

Submitted 20 July, 2023; v1 submitted 1 June, 2022; originally announced June 2022.

CompOSE Reference Manual

Authors: S. Typel, M. Oertel, T. Klähn, D. Chatterjee, V. Dexheimer, C. Ishizuka, M. Mancini, J. Novak, H. Pais, C. Providencia, A. Raduta, M. Servillat, L. Tolos

Abstract: CompOSE (CompStar Online Supernovae Equations of State) is an online repository of equations of state (EoS) for use in nuclear physics and astrophysics, e.g., in the description of compact stars or the simulation of core-collapse supernovae and neutron-star mergers, see http://compose.obspm.fr. The main services, offered via the website, are: a collection of data tables in a flexible and easily ex… ▽ More CompOSE (CompStar Online Supernovae Equations of State) is an online repository of equations of state (EoS) for use in nuclear physics and astrophysics, e.g., in the description of compact stars or the simulation of core-collapse supernovae and neutron-star mergers, see http://compose.obspm.fr. The main services, offered via the website, are: a collection of data tables in a flexible and easily extendable data format for different EoS types and related physical quantities with extensive documentation and referencing; software for download to extract and to interpolate these data and to calculate additional quantities; webtools to generate EoS tables that are customized to the needs of the users and to illustrate dependencies of various EoS quantities in graphical form. This manual is an update of previous versions that are available on the CompOSE website, at arXiv:1307.5715 [astro-ph.SR], and that was originally published in the journal "Physics of Particles and Nuclei" with doi:10.1134/S1063779615040061. It contains a detailed description of the service, containing a general introduction as well as instructions for potential contributors and for users. Short versions of the manual for EoS users and providers will also be available as separate publications. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 91 pages, 20 tables. This manual is an update of arXiv:1307.571

Proto-neutron star evolution with improved charged-current neutrino-nucleon interactions

Authors: A Pascal, J Novak, M Oertel

Abstract: We perform simulations of the Kelvin-Helmholtz cooling phase of proto-neutron stars with a new numerical code in spherical symmetry and using the quasi-static approximation. We use for the first time the full set of charged-current neutrino-nucleon reactions, including neutron decay and modified Urca processes, together with the energy-dependent numerical representation for the inclusion of nuclea… ▽ More We perform simulations of the Kelvin-Helmholtz cooling phase of proto-neutron stars with a new numerical code in spherical symmetry and using the quasi-static approximation. We use for the first time the full set of charged-current neutrino-nucleon reactions, including neutron decay and modified Urca processes, together with the energy-dependent numerical representation for the inclusion of nuclear correlations with random-phase approximation. Moreover, convective motions are taken into account within the mixing-length theory. As we vary the assumptions for computing neutrino-nucleon reaction rates, we show that the dominant effect on the cooling timescale, neutrino signal and composition of the neutrino-driven wind comes from the inclusion of convective motion. Computation of nuclear correlations within the random phase approximation, as compared to mean field approach, has a relatively small impact. △ Less

Submitted 6 January, 2022; originally announced January 2022.

Journal ref: Monthly Notices of the Royal Astronomical Society, In press

Electron Captures and Stability of White Dwarfs

Authors: N. Chamel, L. Perot, A. F. Fantina, D. Chatterjee, S. Ghosh, J. Novak, M. Oertel

Abstract: Electron captures by atomic nuclei in dense matter are among the most important processes governing the late evolution of stars, limiting in particular the stability of white dwarfs. Despite considerable progress in the determination of the equation of state of dense Coulomb plasmas, the threshold electron Fermi energies are still generally estimated from the corresponding $Q$ values in vacuum. Mo… ▽ More Electron captures by atomic nuclei in dense matter are among the most important processes governing the late evolution of stars, limiting in particular the stability of white dwarfs. Despite considerable progress in the determination of the equation of state of dense Coulomb plasmas, the threshold electron Fermi energies are still generally estimated from the corresponding $Q$ values in vacuum. Moreover, most studies have focused on nonmagnetized matter. However, some white dwarfs are endowed with magnetic fields reaching $10^9$ G. Even more extreme magnetic fields might exist in super Chandrasekhar white dwarfs, the progenitors of overluminous type Ia supernovae like SN 2006gz and SN 2009dc. The roles of the dense stellar medium and magnetic fields on the onset of electron captures and on the structure of white dwarfs are briefly reviewed. New analytical formulas are derived to evaluate the threshold density for the onset of electron captures for arbitrary magnetic fields. Their influence on the structure of white dwarfs is illustrated by simple analytical formulas and numerical calculations. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 20 pages, 5 figures, contribution to the Sixteenth Marcel Grossmann Meeting - MG16

Journal ref: The Sixteenth Marcel Grossmann Meeting, pp. 4488-4507 (2023)

Structure of ultra-magnetised neutron stars

Authors: Debarati Chatterjee, Jerome Novak, Micaela Oertel

Abstract: In this review we discuss self-consistent methods to calculate the global structure of strongly magnetised neutron stars within the general-relativistic framework. We outline why solutions in spherical symmetry cannot be applied to strongly magnetised compact stars, and elaborate on a consistent formalism to compute rotating magnetised neutron star models. We also discuss an application of the abo… ▽ More In this review we discuss self-consistent methods to calculate the global structure of strongly magnetised neutron stars within the general-relativistic framework. We outline why solutions in spherical symmetry cannot be applied to strongly magnetised compact stars, and elaborate on a consistent formalism to compute rotating magnetised neutron star models. We also discuss an application of the above full numerical solution for studying the influence of strong magnetic fields on the radius and crust thickness of magnetars. The above technique is also applied to construct a "universal" magnetic field profile inside the neutron star, that may be useful for studies in nuclear physics. The methodology developed here is particularly useful to interpret multi-messenger astrophysical data of strongly magnetised neutron stars. △ Less

Submitted 31 August, 2021; originally announced August 2021.

Comments: 15 pages, 13 figures

Journal ref: Eur. Phys. J. A 57, 249 (2021)

What can be learned from a proto-neutron star's mass and radius?

Authors: Edwan Preau, Aurélien Pascal, Jérôme Novak, Micaela Oertel

Abstract: We make extensive numerical studies of masses and radii of proto-neutron stars during the first second after their birth in core-collapse supernova events. We use a quasi-static approach for the computation of proto-neutron star structure, built on parameterized entropy and electron fraction profiles, that are then evolved with neutrino cooling processes. We vary the equation of state of nuclear m… ▽ More We make extensive numerical studies of masses and radii of proto-neutron stars during the first second after their birth in core-collapse supernova events. We use a quasi-static approach for the computation of proto-neutron star structure, built on parameterized entropy and electron fraction profiles, that are then evolved with neutrino cooling processes. We vary the equation of state of nuclear matter, the proto-neutron star mass and the parameters of the initial profiles, to take into account our ignorance of the supernova progenitor properties. We show that if masses and radii of a proto-neutron star can be determined in the first second after the birth, e.g. from gravitational wave emission, no information could be obtained on the corresponding cold neutron star and therefore on the cold nuclear equation of state. Similarly, it seems unlikely that any property of the proto-neutron star equation of state (hot and not beta-equilibrated) could be determined either, mostly due to the lack of information on the entropy, or equivalently temperature, distribution in such objects. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Improved neutrino-nucleon interactions in dense and hot matter for numerical simulations

Authors: Micaela Oertel, Aurélien Pascal, Marco Mancini, Jerome Novak

Abstract: Neutrinos play an important role in compact star astrophysics: neutrino-heating is one of the main ingredients in core-collapse supernovae, neutrino-matter interactions determine the composition of matter in binary neutron star mergers and have among others a strong impact on conditions for heavy element nucleosynthesis and neutron star cooling is dominated by neutrino emission except for ve… ▽ More Neutrinos play an important role in compact star astrophysics: neutrino-heating is one of the main ingredients in core-collapse supernovae, neutrino-matter interactions determine the composition of matter in binary neutron star mergers and have among others a strong impact on conditions for heavy element nucleosynthesis and neutron star cooling is dominated by neutrino emission except for very old stars. Many works in the last decades have shown that in dense matter medium effects considerably change the neutrino-matter interaction rates, whereas many astrophysical simulations use analytic approximations which are often far from reproducing more complete calculations. In this work we present a scheme which allows to incorporate improved rates, for charged current interactions, into simulations and show as an example some results for core-collapse supernovae, where a noticeable difference is found in the location of the neutrinospheres of the low-energy neutrinos in the early post-bounce phase. △ Less

Submitted 3 September, 2020; v1 submitted 4 March, 2020; originally announced March 2020.

Journal ref: Phys. Rev. C 102, 035802 (2020)

Impact of electron capture rates on nuclei far from stability on core-collapse supernovae

Authors: Aurélien Pascal, Simon Giraud, Anthea Fantina, Francesca Gulminelli, Jerome Novak, Micaela Oertel, Adriana Raduta

Abstract: The impact of electron-capture (EC) cross sections on neutron-rich nuclei on the dynamics of core-collapse during infall and early post-bounce is studied performing spherically symmetric simulations in general relativity using a multigroup scheme for neutrino transport and full nuclear distributions in extended nuclear statistical equilibrium models. We thereby vary the prescription for EC rates o… ▽ More The impact of electron-capture (EC) cross sections on neutron-rich nuclei on the dynamics of core-collapse during infall and early post-bounce is studied performing spherically symmetric simulations in general relativity using a multigroup scheme for neutrino transport and full nuclear distributions in extended nuclear statistical equilibrium models. We thereby vary the prescription for EC rates on individual nuclei, the nuclear interaction for the EoS, the mass model for the nuclear statistical equilibrium distribution and the progenitor model. In agreement with previous works, we show that the individual EC rates are the most important source of uncertainty in the simulations, while the other inputs only marginally influence the results. A recently proposed analytic formula to extrapolate microscopic results on stable nuclei for EC rates to the neutron rich region, with a functional form motivated by nuclear-structure data and parameters fitted from large scale shell model calculations, is shown to lead to a sizable (16%) reduction of the electron fraction at bounce compared to more primitive prescriptions for the rates, leading to smaller inner core masses and slower shock propagation. We show that the EC process involves $\approx$ 130 different nuclear species around 86 Kr mainly in the N = 50 shell closure region, and establish a list of the most important nuclei to be studied in order to constrain the global rates. △ Less

Submitted 2 September, 2019; v1 submitted 12 June, 2019; originally announced June 2019.

Journal ref: Phys. Rev. C 101, 015803 (2020)

Magnetic field distribution in magnetars

Authors: Debarati Chatterjee, Jerome Novak, Micaela Oertel

Abstract: Using an axisymmetric numerical code, we perform an extensive study of the magnetic field configurations in non-rotating neutron stars, varying the mass, magnetic field strength and the equation of state. We find that the monopolar (spherically symmetric) part of the norm of the magnetic field can be described by a single profile, that we fit by a simple eighth-order polynomial, as a function of t… ▽ More Using an axisymmetric numerical code, we perform an extensive study of the magnetic field configurations in non-rotating neutron stars, varying the mass, magnetic field strength and the equation of state. We find that the monopolar (spherically symmetric) part of the norm of the magnetic field can be described by a single profile, that we fit by a simple eighth-order polynomial, as a function of the star's radius. This new generic profile applies remarkably well to all magnetized neutron star configurations built on hadronic equations of state. We then apply this profile to build magnetized neutron stars in spherical symmetry, using a modified Tolman-Oppenheimer-Volkov (TOV) system of equations. This new formalism produces slightly better results in terms of mass-radius diagrams than previous attempts to add magnetic terms to these equations. However, we show that such approaches are less accurate than usual, non-magnetized TOV models, and that consistent models must depart from spherical symmetry. Thus, our `universal' magnetic field profile is intended to serve as a tool for nuclear physicists to obtain estimates of magnetic field inside neutron stars, as a function of radial depth, in order to deduce its influence on composition and related properties. It possesses the advantage of being based on magnetic field distributions from realistic self-consistent computations, which are solutions of Maxwell's equations. △ Less

Submitted 11 June, 2019; v1 submitted 6 August, 2018; originally announced August 2018.

Journal ref: Phys. Rev. C 99, 055811 (2019)

French SKA White Book - The French Community towards the Square Kilometre Array

Authors: F. Acero, J. -T. Acquaviva, R. Adam, N. Aghanim, M. Allen, M. Alves, R. Ammanouil, R. Ansari, A. Araudo, E. Armengaud, B. Ascaso, E. Athanassoula, D. Aubert, S. Babak, A. Bacmann, A. Banday, K. Barriere, F. Bellossi, J. -P. Bernard, M. G. Bernardini, M. Béthermin, E. Blanc, L. Blanchet, J. Bobin, S. Boissier , et al. (153 additional authors not shown)

Abstract: The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive p… ▽ More The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration. △ Less

Submitted 28 March, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferrière, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Deschênes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

A new temperature dependent hyperonic equation of state: application to rotating neutron star models and I-Q-relations

Authors: Miguel Marques, Micaela Oertel, Matthias Hempel, Jérôme Novak

Abstract: In this work we present a newly constructed equation of state (EoS) --applicable to stellar core collapse and neutron star mergers--, including the entire baryon octet. Our EoS is compatible with the main constraints from nuclear physics and, in particular, with a maximum mass for cold beta-equilibrated neutron stars of 2 solar masses in agreement with recent observations. As an application of our… ▽ More In this work we present a newly constructed equation of state (EoS) --applicable to stellar core collapse and neutron star mergers--, including the entire baryon octet. Our EoS is compatible with the main constraints from nuclear physics and, in particular, with a maximum mass for cold beta-equilibrated neutron stars of 2 solar masses in agreement with recent observations. As an application of our new EoS, we compute numerical stationary models for rapidly (rigidly) rotating hot neutron stars. We consider maximum masses of hot stars, such as proto-neutron stars or hypermassive neutron stars in the post-merger phase of binary neutron star coalescence. The universality of I-Q-relations at nonzero temperature for fast rotating models, comparing a purely nuclear EoS with its counterparts containing Lambda-hyperons or the entire baryon octet, respectively, is discussed, too. We find that the I-Q universality is broken when thermal effects become important, whatever the value of entropy gradients in our models. Thus, the use of I-Q relations for the analysis of proto-neutron stars or merger remnant data, including gravitational wave signals from the last stages of binary neutron star mergers, should be regarded with care. △ Less

Submitted 4 October, 2017; v1 submitted 9 June, 2017; originally announced June 2017.

Comments: 20 pages, 13 figures, v2: additional discussion of EoS properties and I-Q fit, some references added, accepted for publication in PRC

Journal ref: Phys. Rev. C 96, 045806 (2017)

On the maximum mass of magnetised white dwarfs

Authors: D. Chatterjee, A. F. Fantina, N. Chamel, J. Novak, M. Oertel

Abstract: We develop a detailed and self-consistent numerical model for extremely-magnetised white dwarfs, which have been proposed as progenitors of overluminous Type Ia supernovae. This model can describe fully-consistent equilibria of magnetic stars in axial symmetry, with rotation, general-relativistic effects and realistic equations of state (including electron-ion interactions and taking into account… ▽ More We develop a detailed and self-consistent numerical model for extremely-magnetised white dwarfs, which have been proposed as progenitors of overluminous Type Ia supernovae. This model can describe fully-consistent equilibria of magnetic stars in axial symmetry, with rotation, general-relativistic effects and realistic equations of state (including electron-ion interactions and taking into account Landau quantisation of electrons due to the magnetic field). We study the influence of each of these ingredients onto the white dwarf structure and, in particular, on their maximum mass. We perform an extensive stability analysis of such objects, with their highest surface magnetic fields reaching $\sim 10^{13}~G$ (at which point the star adopts a torus-like shape). We confirm previous speculations that although very massive strongly magnetised white dwarfs could potentially exist, the onset of electron captures and pycnonuclear reactions may severely limit their stability. Finally, the emission of gravitational waves by these objects is addressed, showing no possibility of detection by the currently planned space-based detector eLISA. △ Less

Submitted 13 October, 2016; originally announced October 2016.

Comments: submitted to MNRAS on 24.2.2016

Global numerical simulations of the rise of vortex-mediated pulsar glitches in full general relativity

Authors: Aurélien Sourie, Nicolas Chamel, Jérôme Novak, Micaela Oertel

Abstract: In this paper, we study in detail the role of general relativity on the global dynamics of giant pulsar glitches as exemplified by Vela. For this purpose, we carry out numerical simulations of the spin up triggered by the sudden unpinning of superfluid vortices. In particular, we compute the exchange of angular momentum between the core neutron superfluid and the rest of the star within a two-flui… ▽ More In this paper, we study in detail the role of general relativity on the global dynamics of giant pulsar glitches as exemplified by Vela. For this purpose, we carry out numerical simulations of the spin up triggered by the sudden unpinning of superfluid vortices. In particular, we compute the exchange of angular momentum between the core neutron superfluid and the rest of the star within a two-fluid model including both (non-dissipative) entrainment effects and (dissipative) mutual friction forces. Our simulations are based on a quasi-stationary approach using realistic equations of state (EoSs). We show that the evolution of the angular velocities of both fluids can be accurately described by an exponential law. The associated characteristic rise time $τ_{\text{r}}$, which can be precisely computed from stationary configurations only, has a form similar to that obtained in the Newtonian limit. However, general relativity changes the structure of the star and leads to additional couplings between the fluids due to frame-dragging effects. As a consequence, general relativity can have a large impact on the actual value of $τ_{\text{r}}$: the errors incurred by using Newtonian gravity are thus found to be as large as $\sim 40 \%$ for the models considered. Values of the rise time are calculated for Vela and compared with current observational limits. Finally, we study the amount of gravitational waves emitted during a glitch. Simple expressions are obtained for the corresponding characteristic amplitudes and frequencies. The detectability of glitches through gravitational wave observatories is briefly discussed. △ Less

Submitted 3 November, 2016; v1 submitted 27 July, 2016; originally announced July 2016.

Comments: 19 pages, 12 figures, minor changes to match version to be published in MNRAS

Theory of Gravitational Waves

Authors: Alexandre Le Tiec, Jérôme Novak

Abstract: The existence of gravitational radiation is a natural prediction of any relativistic description of the gravitational interaction. In this chapter, we focus on gravitational waves, as predicted by Einstein's general theory of relativity. First, we introduce those mathematical concepts that are necessary to properly formulate the physical theory, such as the notions of manifold, vector, tensor, met… ▽ More The existence of gravitational radiation is a natural prediction of any relativistic description of the gravitational interaction. In this chapter, we focus on gravitational waves, as predicted by Einstein's general theory of relativity. First, we introduce those mathematical concepts that are necessary to properly formulate the physical theory, such as the notions of manifold, vector, tensor, metric, connection and curvature. Second, we motivate, formulate and then discuss Einstein's equation, which relates the geometry of spacetime to its matter content. Gravitational waves are later introduced as solutions of the linearized Einstein equation around flat spacetime. These waves are shown to propagate at the speed of light and to possess two polarization states. Gravitational waves can interact with matter, allowing for their direct detection by means of laser interferometers. Finally, Einstein's quadrupole formulas are derived and used to show that nonspherical compact objects moving at relativistic speeds are powerful gravitational wave sources. △ Less

Submitted 13 July, 2016; originally announced July 2016.

Comments: 40 pages, 10 figures, 2 tables; to appear in An Overview of Gravitational Waves: Theory and Detection, edited by G. Auger and E. Plagnol (World Scientific, 2016)

Numerical models for stationary superfluid neutron stars in general relativity with realistic equations of state

Authors: Aurélien Sourie, Micaela Oertel, Jérôme Novak

Abstract: We present a numerical model for uniformly rotating superfluid neutron stars, for the first time with realistic microphysics including entrainment, in a fully general relativistic framework. We compute stationary and axisymmetric configurations of neutron stars composed of two fluids, namely superfluid neutrons and charged particles (protons and electrons), rotating with different rates around a c… ▽ More We present a numerical model for uniformly rotating superfluid neutron stars, for the first time with realistic microphysics including entrainment, in a fully general relativistic framework. We compute stationary and axisymmetric configurations of neutron stars composed of two fluids, namely superfluid neutrons and charged particles (protons and electrons), rotating with different rates around a common axis. Both fluids are coupled by entrainment, a non-dissipative interaction which in case of a non-vanishing relative velocity between the fluids, causes the fluid momenta being not aligned with the respective fluid velocities. We extend the formalism by Comer and Joynt (2003) in order to calculate the equation of state (EoS) and entrainment parameters for an arbitrary relative velocity. The resulting entrainment matrix fulfills all necessary sum rules and in the limit of small relative velocity our results agree with Fermi liquid theory ones, derived to lowest order in the velocity. This formalism is applied to two new nuclear equations of state which are implemented in the numerical model. We are able to obtain precise equilibrium configurations. Resulting density profiles and moments of inertia are discussed employing both EoSs, showing the impact of entrainment and the dependence on the EoS. △ Less

Submitted 3 November, 2016; v1 submitted 18 February, 2016; originally announced February 2016.

Comments: 18 pages, 10 figures, minor changes to match published version in PRD, a typo present in Eq.A1 in the published version has been corrected

Journal ref: Phys. Rev. D 93, 083004(2016)

The explosion mechanism of core-collapse supernovae: progress in supernova theory and experiments

Authors: Thierry Foglizzo, Rémi Kazeroni, Jérôme Guilet, Frédéric Masset, Matthias González, Brendan K. Krueger, Jérôme Novak, Micaela Oertel, Jérôme Margueron, Julien Faure, Noël Martin, Patrick Blottiau, Bruno Peres, Gilles Durand

Abstract: The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus… ▽ More The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus on the mechanism which transforms the collapse into an explosion. Progress in understanding this mechanism is reviewed with particular attention to its asymmetric character. We highlight a series of successful studies connecting observations of supernova remnants and pulsars properties to the theory of core-collapse using numerical simulations. The encouraging results from first principles models in axisymmetric simulations is tempered by new puzzles in 3D. The diversity of explosion paths and the dependence on the pre-collapse stellar structure is stressed, as well as the need to gain a better understanding of hydrodynamical and MHD instabilities such as SASI and neutrino-driven convection. The shallow water analogy of shock dynamics is presented as a comparative system where buoyancy effects are absent. This dynamical system can be studied numerically and also experimentally with a water fountain. The potential of this complementary research tool for supernova theory is analyzed. We also review its potential for public outreach in science museums. △ Less

Submitted 6 February, 2015; v1 submitted 6 January, 2015; originally announced January 2015.

Comments: 19 pages, 6 figures, invited review accepted for publication in PASA

Consistent neutron star models with magnetic field dependent equations of state

Authors: Debarati Chatterjee, Thomas Elghozi, Jerome Novak, Micaela Oertel

Abstract: We present a self-consistent model for the study of the structure of a neutron star in strong magnetic fields. Starting from a microscopic Lagrangian, this model includes the effect of the magnetic field on the equation of state, the interaction of the electromagnetic field with matter (magnetisation), and anisotropies in the energy-momentum tensor, as well as general relativistic aspects. We buil… ▽ More We present a self-consistent model for the study of the structure of a neutron star in strong magnetic fields. Starting from a microscopic Lagrangian, this model includes the effect of the magnetic field on the equation of state, the interaction of the electromagnetic field with matter (magnetisation), and anisotropies in the energy-momentum tensor, as well as general relativistic aspects. We build numerical axisymmetric stationary models and show the applicability of the approach with one example quark matter equation of state (EoS) often employed in the recent literature for studies of strongly magnetised neutron stars. For this EoS, the effect of inclusion of magnetic field dependence or the magnetisation do not increase the maximum mass significantly in contrast to what has been claimed by previous studies. △ Less

Submitted 23 October, 2014; originally announced October 2014.

Comments: 13 pages, 8 figures

$f(R)$ gravity: scalar perturbations in the late Universe

Authors: Maxim Eingorn, Jan Novák, Alexander Zhuk

Abstract: In this paper we study scalar perturbations of the metric for nonlinear $f(R)$ models. We consider the Universe at the late stage of its evolution and deep inside the cell of uniformity. We investigate the astrophysical approach in the case of Minkowski spacetime background and two cases in the cosmological approach, the large scalaron mass approximation and the quasi-static approximation, getting… ▽ More In this paper we study scalar perturbations of the metric for nonlinear $f(R)$ models. We consider the Universe at the late stage of its evolution and deep inside the cell of uniformity. We investigate the astrophysical approach in the case of Minkowski spacetime background and two cases in the cosmological approach, the large scalaron mass approximation and the quasi-static approximation, getting explicit expressions for scalar perturbations for both these cases. In the most interesting quasi-static approximation, the scalar perturbation functions depend on both the nonlinearity function $f(R)$ and the scale factor $a$. Hence, we can study the dynamical behavior of the inhomogeneities (e.g., galaxies and dwarf galaxies) including into consideration their gravitational attraction and the cosmological expansion, and also taking into account the effects of nonlinearity. Our investigation is valid for functions $f(R)$ which have stable de Sitter points in future with respect to the present time, that is typical for the most popular $f(R)$ models. △ Less

Submitted 2 September, 2014; v1 submitted 21 January, 2014; originally announced January 2014.

Comments: 9 pages, no figures

Journal ref: The European Physical Journal C 74, 3005 (2014)

General relativistic neutrino transport using spectral methods

Authors: Bruno Peres, Andrew Jason Penner, Jerome Novak, Silvano Bonazzola

Abstract: We present a new code, Lorene's Ghost (for Lorene's gravitational handling of spectral transport) developed to treat the problem of neutrino transport in supernovae with the use of spectral methods. First, we derive the expression for the nonrelativistic Liouville operator in doubly spherical coordinates (r, theta, phi, epsilon, Theta, Phi)$, and further its general relativistic counterpart. We us… ▽ More We present a new code, Lorene's Ghost (for Lorene's gravitational handling of spectral transport) developed to treat the problem of neutrino transport in supernovae with the use of spectral methods. First, we derive the expression for the nonrelativistic Liouville operator in doubly spherical coordinates (r, theta, phi, epsilon, Theta, Phi)$, and further its general relativistic counterpart. We use the 3 + 1 formalism with the conformally flat approximation for the spatial metric, to express the Liouville operator in the Eulerian frame. Our formulation does not use any approximations when dealing with the angular arguments (theta, phi, Theta, Phi), and is fully energy-dependent. This approach is implemented in a spherical shell, using either Chebyshev polynomials or Fourier series as decomposition bases. It is here restricted to simplified collision terms (isoenergetic scattering) and to the case of a static fluid. We finish this paper by presenting test results using basic configurations, including general relativistic ones in the Schwarzschild metric, in order to demonstrate the convergence properties, the conservation of particle number and correct treatment of some general-relativistic effects of our code. The use of spectral methods enables to run our test cases in a six-dimensional setting on a single processor. △ Less

Submitted 29 January, 2014; v1 submitted 5 July, 2013; originally announced July 2013.

Comments: match published version

Journal ref: Class. Quantum Grav. 31 (2014) 045012

Influence of pions and hyperons on stellar black hole formation

Authors: Bruno Peres, Micaela Oertel, Jerome Novak

Abstract: We present numerical simulations of stellar core-collapse with spherically symmetric, general relativistic hydrodynamics up to black hole formation. Using the CoCoNuT code, with a newly developed grey leakage scheme for the neutrino treatment, we investigate the effects of including pions and Λ-hyperons into the equation of state at high densities and temperatures on the black hole formation proce… ▽ More We present numerical simulations of stellar core-collapse with spherically symmetric, general relativistic hydrodynamics up to black hole formation. Using the CoCoNuT code, with a newly developed grey leakage scheme for the neutrino treatment, we investigate the effects of including pions and Λ-hyperons into the equation of state at high densities and temperatures on the black hole formation process. Results show non-negligible differences between the models with reference equation of state without any additional particles and models with the extended ones. For the latter, the maximum masses supported by the proto-neutron star are smaller and the collapse to a black hole occurs earlier. A phase transition to hyperonic matter is observed when the progenitor allows for a high enough accretion rate onto the proto-neutron star. Rough estimates of neutrino luminosity from these collapses are given, too. △ Less

Submitted 10 January, 2013; v1 submitted 28 October, 2012; originally announced October 2012.

Comments: 22 pages, 10 figures. Minor changes

Journal ref: Phys. Rev. D 87, 043006 (2013)

Constraining the nature of the Galactic center black hole Sgr A* with present and future observations

Authors: F. H. Vincent, E. Gourgoulhon, O. Straub, M. Abramowicz, J. Novak, T. Paumard, G. Perrin

Abstract: The Galactic center is an ideal laboratory to study strong-field general relativistic phenomena, as the supermassive black hole Sgr A* has the biggest angular Schwarzschild radius among all black holes. This article presents three different ways of using the immediate surroundings of Sgr A* as a probe in order to either constrain its spin, or even test the very nature of this compact object. The Galactic center is an ideal laboratory to study strong-field general relativistic phenomena, as the supermassive black hole Sgr A* has the biggest angular Schwarzschild radius among all black holes. This article presents three different ways of using the immediate surroundings of Sgr A* as a probe in order to either constrain its spin, or even test the very nature of this compact object. △ Less

Submitted 3 September, 2012; originally announced September 2012.

Comments: 4 pages, 2 figures, Proceedings of the French Astronomy and Astrophysics Society Conference

3+1 geodesic equation and images in numerical spacetimes

Authors: Frederic H. Vincent, Eric Gourgoulhon, Jérôme Novak

Abstract: The equations governing null and timelike geodesics are derived within the 3+1 formalism of general relativity. In addition to the particle's position, they encompass an evolution equation for the particle's energy leading to a 3+1 expression of the redshift factor for photons. An important application is the computation of images and spectra in spacetimes arising from numerical relativity, via th… ▽ More The equations governing null and timelike geodesics are derived within the 3+1 formalism of general relativity. In addition to the particle's position, they encompass an evolution equation for the particle's energy leading to a 3+1 expression of the redshift factor for photons. An important application is the computation of images and spectra in spacetimes arising from numerical relativity, via the ray-tracing technique. This is illustrated here by images of numerically computed stationary neutron stars and dynamical neutron stars collapsing to a black hole. △ Less

Submitted 24 October, 2012; v1 submitted 20 August, 2012; originally announced August 2012.

Comments: 20 pages, 4 figures, accepted by CQG, with substantial improvements from referees remarks

Journal ref: Class. Quantum Grav., 29, 245005 (2012)

An extended equation of state for core-collapse simulations

Authors: M. Oertel, A. F. Fantina, J. Novak

Abstract: In stellar core-collapse events matter is heated and compressed to densities above nuclear matter saturation density. For progenitors stars with masses above about 25 solar masses, which eventually form a black hole, the temperatures and densities reached during the collapse are so high that a traditional description in terms of electrons, nuclei, and nucleons is no longer adequate. We present her… ▽ More In stellar core-collapse events matter is heated and compressed to densities above nuclear matter saturation density. For progenitors stars with masses above about 25 solar masses, which eventually form a black hole, the temperatures and densities reached during the collapse are so high that a traditional description in terms of electrons, nuclei, and nucleons is no longer adequate. We present here an improved equation of state which contains in addition pions and hyperons. They become abundant in the high temperature and density regime. We study the different constraints on such an equation of state, coming from both hyperonic data and observations of neutron star properties. In order to test the zero-temperature versions, we perform numerical simulations of the collapse of a neutron star with such additional particles to a black hole. We discuss the influence of the additional particles on the thermodynamic properties within the hot versions of the equation of state and we show that in regimes relevant to core-collapse and black hole formation, the effects of pions and hyperons on pressure, internal energy and sound speed are not negligible. △ Less

Submitted 13 February, 2012; originally announced February 2012.

Comments: 20 pages

Journal ref: Phys.Rev. C85 (2012) 055806

Testing MOND in the Solar System

Authors: Luc Blanchet, Jerome Novak

Abstract: The Modified Newtonian Dynamics (MOND) generically predicts a violation of the strong version of the equivalence principle. As a result the gravitational dynamics of a system depends on the external gravitational field in which the system is embedded. This so-called external field effect is shown to imply the existence of an anomalous quadrupolar correction, along the direction of the external gal… ▽ More The Modified Newtonian Dynamics (MOND) generically predicts a violation of the strong version of the equivalence principle. As a result the gravitational dynamics of a system depends on the external gravitational field in which the system is embedded. This so-called external field effect is shown to imply the existence of an anomalous quadrupolar correction, along the direction of the external galactic field, in the gravitational potential felt by planets in the Solar System. We compute this effect by a numerical integration of the MOND equation in the presence of an external field, and deduce the secular precession of the perihelion of planets induced by this effect. We find that the precession effect is rather large for outer gaseous planets, and in the case of Saturn is comparable to, and in some cases marginally excluded by published residuals of precession permitted by the best planetary ephemerides. △ Less

Submitted 29 May, 2011; originally announced May 2011.

Comments: 8 pages, two figures, to appear in the proceedings of the "46th Rencontres de Moriond and GPhyS Colloquium on Gravitational Waves and Experimental Gravity", la Thuile, March 20-27, 2011

External field effect of modified Newtonian dynamics in the Solar system

Authors: Luc Blanchet, Jerome Novak

Abstract: The Modified Newtonian Dynamics (MOND) has been formulated as a modification of the Poisson equation for the Newtonian gravitational field. This theory generically predicts a violation of the strong version of the equivalence principle, and as a result the gravitational dynamics of a system depends on the external gravitational field in which the system is embedded. This so-called external field e… ▽ More The Modified Newtonian Dynamics (MOND) has been formulated as a modification of the Poisson equation for the Newtonian gravitational field. This theory generically predicts a violation of the strong version of the equivalence principle, and as a result the gravitational dynamics of a system depends on the external gravitational field in which the system is embedded. This so-called external field effect has been recently shown to imply the existence of an anomalous quadrupolar correction, along the direction of the external Galactic field, in the gravitational potential felt by planets in the Solar System. In this paper we confirm the existence of this effect by a numerical integration of the MOND equation in the presence of an external field, and compute the secular precession of the perihelion of planets induced by this effect. We find that the precession effect is rather large for outer gaseous planets, and in the case of Saturn is comparable to published residuals of precession obtained by Saturn range tracking data. The effect is much smaller for inner planets, but in the case of the Earth it appears to be in conflict for most of the MOND functions $μ(y)$ with the very good constraint on the perihelion precession obtained from Jupiter VLBI data. The MOND functions that are compatible with this constraint appear to have a very rapid transition from the MONDian regime to the Newtonian one. △ Less

Submitted 3 January, 2011; v1 submitted 7 October, 2010; originally announced October 2010.

Comments: 15 pages, 8 figures, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 412, 4 (2011) 2530-2542

A new spectral apparent horizon finder for 3D numerical relativity

Authors: Lap-Ming Lin, Jerome Novak

Abstract: We present a new spectral-method-based algorithm for finding apparent horizons in three-dimensional space-like hypersurfaces without symmetries. While there are already a wide variety of algorithms for finding apparent horizons, our new algorithm does not suffer from the same weakness as previous spectral apparent horizon finders: namely the monopolar coefficient ($\ell=0$ in terms of the spheri… ▽ More We present a new spectral-method-based algorithm for finding apparent horizons in three-dimensional space-like hypersurfaces without symmetries. While there are already a wide variety of algorithms for finding apparent horizons, our new algorithm does not suffer from the same weakness as previous spectral apparent horizon finders: namely the monopolar coefficient ($\ell=0$ in terms of the spherical harmonics decomposition) needed to be determined by a root-finding procedure. Hence, this leads to a much faster and more robust spectral apparent horizon finder. The finder is tested with the Kerr-Schild and Brill-Lindquist data. Our finder is accurate and is as efficient as the currently fastest methods developed recently by Schnetter (2003 Class. Quantum Grav. {\bf 20}, 4719) and Thornburg (2004 Class. Quantum Grav. {\bf 21}, 743). At typical resolutions it takes only 0.5 second to find the apparent horizon of a Kerr-Schild black hole with $a=0.9M$ to the accuracy $\sim 10^{-5}$ for the fractional error in the horizon's location on a 2 GHz processor. △ Less

Submitted 2 May, 2007; v1 submitted 6 February, 2007; originally announced February 2007.

Comments: Minor changes and references added. Published in Class. Quantum Grav

Journal ref: Class.Quant.Grav.24:2665-2676,2007

Rotating star initial data for a constrained scheme in numerical relativity

Authors: Lap-Ming Lin, Jerome Novak

Abstract: A new numerical code for computing stationary axisymmetric rapidly rotating stars in general relativity is presented. The formulation is based on a fully constrained-evolution scheme for 3+1 numerical relativity using the Dirac gauge and maximal slicing. We use both the polytropic and MIT bag model equations of state to demonstrate that the code can construct rapidly rotating neutron star and st… ▽ More A new numerical code for computing stationary axisymmetric rapidly rotating stars in general relativity is presented. The formulation is based on a fully constrained-evolution scheme for 3+1 numerical relativity using the Dirac gauge and maximal slicing. We use both the polytropic and MIT bag model equations of state to demonstrate that the code can construct rapidly rotating neutron star and strange star models. We compare numerical models obtained by our code and a well-established code, which uses a different gauge condition, and show that the two codes agree to high accuracy. △ Less

Submitted 30 May, 2006; v1 submitted 13 March, 2006; originally announced March 2006.

Comments: Minor changes and one figure added. Version accepted for publication in Class. Quant. Grav

Journal ref: Class.Quant.Grav.23:4545,2006

Relativistic numerical models for stationary superfluid Neutron Stars

Authors: Reinhard Prix, Jerome Novak, G. L. Comer

Abstract: We have developed a theoretical model and a numerical code for stationary rotating superfluid neutron stars in full general relativity. The underlying two-fluid model is based on Carter's covariant multi-fluid hydrodynamic formalism. The two fluids, representing the superfluid neutrons on one hand, and the protons and electrons on the other, are restricted to uniform rotation around a common axi… ▽ More We have developed a theoretical model and a numerical code for stationary rotating superfluid neutron stars in full general relativity. The underlying two-fluid model is based on Carter's covariant multi-fluid hydrodynamic formalism. The two fluids, representing the superfluid neutrons on one hand, and the protons and electrons on the other, are restricted to uniform rotation around a common axis, but are allowed to have different rotation rates. We have performed extensive tests of the numerical code, including quantitative comparisons to previous approximative results for these models. The results presented here are the first ``exact'' calculations of such models in the sense that no approximations (other than that inherent in a discretized numerical treatment) are used. Using this code we reconfirm the existence of prolate-oblate shaped configurations. We studied the dependency of the Kepler rotation limit and of the mass-density relation on the relative rotation rate. We further demonstrate how one can simulate a (albeit fluid) neutron-star ``crust'' by letting one fluid extend further outwards than the other, which results in interesting cases where the Kepler limit is actually determined by the outermost but slower fluid. △ Less

Submitted 6 February, 2005; v1 submitted 5 October, 2004; originally announced October 2004.

Comments: 18 pages, 9 figures, RevTex4; to appear in Phys.Rev.D; v2: improved figures, more discussion of physical results

Journal ref: Phys.Rev. D71 (2005) 043005

"Mariage des Maillages": A new numerical approach for 3D relativistic core collapse simulations

Authors: Harald Dimmelmeier, Jerome Novak, Jose A. Font, Jose M. Ibanez, Ewald Mueller

Abstract: We present a new 3D general relativistic hydrodynamics code for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. It uses spectral methods for solving the metric equations, assuming the conformal flatness approximation for the three-metric. The matter equations are solved by high-resolution shock-capturing schemes. We… ▽ More We present a new 3D general relativistic hydrodynamics code for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. It uses spectral methods for solving the metric equations, assuming the conformal flatness approximation for the three-metric. The matter equations are solved by high-resolution shock-capturing schemes. We demonstrate that the combination of a finite difference grid and a spectral grid can be successfully accomplished. This "Mariage des Maillages" (French for grid wedding) approach results in high accuracy of the metric solver and allows for fully 3D applications using computationally affordable resources, and ensures long term numerical stability of the evolution. We compare our new approach to two other, finite difference based, methods to solve the metric equations. A variety of tests in 2D and 3D is presented, involving highly perturbed neutron star spacetimes and (axisymmetric) stellar core collapse, demonstrating the ability to handle spacetimes with and without symmetries in strong gravity. These tests are also employed to assess gravitational waveform extraction, which is based on the quadrupole formula. △ Less

Submitted 2 March, 2005; v1 submitted 8 July, 2004; originally announced July 2004.

Comments: 29 pages, 16 figures; added more information about convergence tests and grid setup

Journal ref: Phys.Rev. D71 (2005) 064023

Stationary structure of relativistic superfluid neutron stars

Authors: R. Prix, J. Novak, G. L. Comer

Abstract: We describe recent progress in the numerical study of the structure of rapidly rotating superfluid neutron star models in full general relativity. The superfluid neutron star is described by a model of two interpenetrating and interacting fluids, one representing the superfluid neutrons and the second consisting of the remaining charged particles (protons, electrons, muons). We consider general… ▽ More We describe recent progress in the numerical study of the structure of rapidly rotating superfluid neutron star models in full general relativity. The superfluid neutron star is described by a model of two interpenetrating and interacting fluids, one representing the superfluid neutrons and the second consisting of the remaining charged particles (protons, electrons, muons). We consider general stationary configurations where the two fluids can have different rotation rates around a common rotation axis. The previously discovered existence of configurations with one fluid in a prolate shape is confirmed. △ Less

Submitted 29 November, 2002; originally announced November 2002.

Comments: 5 pages, 2 figures. Conference proceedings for the 26th Spanish Relativity Meeting (ERE 2002), Menorca, Spain, 22-24 Sept. 2002

Velocity-induced collapses of stable neutron stars

Authors: Jerome Novak

Abstract: The collapse of spherical neutron stars is studied in General Relativity. The initial state is a stable neutron star to which an inward radial kinetic energy has been added through some velocity profile. For two different equations of state and two different shapes of velocity profiles, it is found that neutron stars can collapse to black holes for high enough inward velocities, provided that th… ▽ More The collapse of spherical neutron stars is studied in General Relativity. The initial state is a stable neutron star to which an inward radial kinetic energy has been added through some velocity profile. For two different equations of state and two different shapes of velocity profiles, it is found that neutron stars can collapse to black holes for high enough inward velocities, provided that their masses are higher than some minimal value, depending on the equation of state. For a polytropic equation of state of the form $p=Kρ^γ$, with $γ= 2$ it is found to be $1.16 (\frac{K}{0.1})^{0.5} \msol$, whereas for a more realistic one, it reads $0.36 \msol $. In some cases of collapse forming a black hole, part of the matter composing the initial neutron star can be ejected through a shock, leaving only a fraction of the initial mass to form a black hole. Therefore, black holes of very small masses can be formed and, in particular, the mass scaling relation, as a function of initial velocity, takes the form discovered by Choptuik for critical collapses. △ Less

Submitted 12 July, 2001; originally announced July 2001.

Comments: 8 pages, 7 figures, uses aa.cls, accepted for publication in Astronomy and Astrophysics

Journal ref: Astron.Astrophys.376:606-613,2001

Gravitational waves from the collapse and bounce of a stellar core in tensor-scalar gravity

Authors: J. Novak, J. M. Ibanez

Abstract: Tensor-scalar theory of gravity allows the generation of gravitational waves from astrophysical sources, like Supernovæ, even in the spherical case. That motivated us to study the collapse of a degenerate stellar core, within tensor-scalar gravity, leading to the formation of a neutron star through a bounce and the formation of a shock. We discuss in this paper the effects of the scalar field on… ▽ More Tensor-scalar theory of gravity allows the generation of gravitational waves from astrophysical sources, like Supernovæ, even in the spherical case. That motivated us to study the collapse of a degenerate stellar core, within tensor-scalar gravity, leading to the formation of a neutron star through a bounce and the formation of a shock. We discuss in this paper the effects of the scalar field on the evolution of the system, as well as the appearance of strong non-perturbative effects of this scalar field (the so-called ``spontaneous scalarization''). As a main result, we describe the resulting gravitational monopolar radiation (form and amplitude) and discuss the possibility of its detection by the gravitational detectors currently under construction, taking into account the existing constraints on the scalar field. From the numerical point of view it is worthy to point out that we have developed a combined code which uses pseudo-spectral methods, for the evolution of the scalar field, and High Resolution Shock-Capturing schemes, for the evolution of the hydrodynamical system. Although this code has been used to integrate the field equations of that theory of gravity, in the spherically symmetric case, a by-product of the present work is to gain experience for an ulterior extension to multidimensional problems in Numerical Relativity of such numerical strategy. △ Less

Submitted 16 November, 1999; originally announced November 1999.

Comments: 30 pages, 17 figures, uses AASLATEX v4.0 Accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J. 533 (2000) 392-405

Neutron star transition to strong-scalar-field state in tensor scalar gravity

Authors: J. Novak

Abstract: Spherical neutron star models are studied within tensor-scalar theories of gravity. Particularly, it is shown that, under some conditions on the second derivative of the coupling function and on star's mass, for a given star there exist two strong-scalar-field solutions as well as the usual weak-field one. This last solution happens to be unstable and a star, becoming massive enough to allow for… ▽ More Spherical neutron star models are studied within tensor-scalar theories of gravity. Particularly, it is shown that, under some conditions on the second derivative of the coupling function and on star's mass, for a given star there exist two strong-scalar-field solutions as well as the usual weak-field one. This last solution happens to be unstable and a star, becoming massive enough to allow for all three solutions, evolves to reach one of the strong field configurations. This transition is dynamically computed and it appears that the star radiates away the difference in energy between both states (a few $10^{-3} M_\odot c^2$) as gravitational radiation. Since part of the energy ($\sim 10^{-5} M_\odot c^2$) is injected into the star as kinetic energy, the velocity of star's surface can reach up to $10^{-2} c$. The waveform of this monopolar radiation is shown as well as the oscillations undergone by the star. These oscillations are also studied within the slowly-rotating approximation, in order to estimate an order of magnitude of the resulting quadrupolar radiation. △ Less

Submitted 4 June, 1998; originally announced June 1998.

Comments: 14 pages, 10 figures, accepted for publication in Physical Review D

Journal ref: Phys.Rev. D58 (1998) 064019

Spherical neutron star collapse toward a black hole in tensor-scalar theory of gravity

Authors: J. Novak

Abstract: Complete tensor-scalar and hydrodynamic equations are presented and integrated, for a self-gravitating perfect fluid. The initial conditions describe unstable-equilibrium neutron star configuration, with a polytropic equation of state. They are necessary in order to follow the gravitational collapse (including full hydrodynamics) of this star toward a black hole and to study the resulting scalar… ▽ More Complete tensor-scalar and hydrodynamic equations are presented and integrated, for a self-gravitating perfect fluid. The initial conditions describe unstable-equilibrium neutron star configuration, with a polytropic equation of state. They are necessary in order to follow the gravitational collapse (including full hydrodynamics) of this star toward a black hole and to study the resulting scalar gravitational wave. The amplitude of this wave, as well as the radiated energy dramatically increase above some critical value of the parameter of the coupling function, due to the spontaneous scalarization, an effect not present in Brans-Dicke theory. In most cases, the pressure of the collapsing fluid does not have a significant impact on the resulting signal. These kind of sources are not likely to be observed by future laser interferometric detectors (such as VIRGO or LIGO) of gravitational waves, if they are located at more than a few 100 kpc. However, spontaneous scalarization could be constrained if such a gravitational collapse is detected by its quadrupolar gravitational signal, since this latter is quite lower than the monopolar one. △ Less

Submitted 23 January, 1998; v1 submitted 18 July, 1997; originally announced July 1997.

Comments: 23 pages, LaTeX, 21 PostScript figures, uses RevTeX, accepted for publication in Physical Review D Numerical appendix added, sec.III rewritten, figures changed

Journal ref: Phys.Rev. D57 (1998) 4789-4801

ROTATING NEUTRON STAR MODELS WITH MAGNETIC FIELD

Authors: M. Bocquet, S. Bonazzola, E. Gourgoulhon, J. Novak

Abstract: We present the first numerical solutions of the coupled Einstein-Maxwell equations describing rapidly rotating neutron stars endowed with a magnetic field. These solutions are fully relativistic and self-consistent, all the effects of the electromagnetic field on the star's equilibrium (Lorentz force, spacetime curvature generated by the electromagnetic stress-energy) being taken into account. T… ▽ More We present the first numerical solutions of the coupled Einstein-Maxwell equations describing rapidly rotating neutron stars endowed with a magnetic field. These solutions are fully relativistic and self-consistent, all the effects of the electromagnetic field on the star's equilibrium (Lorentz force, spacetime curvature generated by the electromagnetic stress-energy) being taken into account. The magnetic field is axisymmetric and poloidal. Five dense matter equations of state are employed. The partial differential equation system is integrated by means of a pseudo-spectral method. Various tests passed by the numerical code are presented. The effects of the magnetic field on neutron stars structure are then investigated, especially by comparing magnetized and non-magnetized configurations with the same baryon number. The deformation of the star induced by the magnetic field is important only for huge values of B (B>10^{10} T). The maximum mass as well as the maximum rotational velocity are found to increase with the magnetic field. The maximum allowable poloidal magnetic field is of the order of 10^{14} T (10^{18} G) and is reached when the magnetic pressure is comparable to the fluid pressure at the centre of the star. For such values, the maximum mass of neutron stars is found to increase by 13 to 29 % (depending upon the EOS) with respect to the maximum mass of non-magnetized stars. △ Less

Submitted 24 March, 1995; originally announced March 1995.

Comments: 19 pages, 26 PostScript figures (available as uuencoded compressed tar file), uses psfig.tex, l-aa.sty (LaTeX style file for Astronomy & Astrophysics, vers. 3.0). Accepted on March 1995 for publication in Astronomy & Astrophysics

Report number: Preprint Observatoire de Paris

Journal ref: Astron. Astrophys. 301, 757 (1995)