-
Electron Heating in the Trans-Relativistic Perpendicular Shocks of Tilted Accretion Flows
Authors:
Lorenzo Sironi,
Aaron Tran
Abstract:
General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks -- where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin -- commonly display standing shocks, within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable t…
▽ More
General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks -- where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin -- commonly display standing shocks, within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have trans-relativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure is $β_\mathrm{pi1} \sim 1-8$), and low sonic Mach number (the ratio of shock speed to sound speed is $M_s \sim 1-5$). We study such shocks with two-dimensional particle-in-cell simulations and we quantify the efficiency and mechanisms of electron heating, for the special case of pre-shock magnetic fields perpendicular to the shock direction of propagation. We find that the post-shock electron temperature $T_\mathrm{e2}$ exceeds the adiabatic expectation $T_\mathrm{e2,ad}$ by an amount $T_\mathrm{e2}/T_\mathrm{e2,ad} - 1 \simeq 0.0016 M_s^{3.6}$, nearly independent of the plasma beta and of the pre-shock electron-to-ion temperature ratio $T_\mathrm{e1}/T_\mathrm{i1}$, which we vary from $0.1$ to unity. We investigate the heating physics for $M_s \sim 5-6$ and find that electron super-adiabatic heating is governed by magnetic pumping at $T_\mathrm{e1}/T_\mathrm{i1}=1$, whereas heating by $B$-parallel electric fields (i.e., parallel to the local magnetic field) dominates at $T_\mathrm{e1}/T_\mathrm{i1}=0.1$. Our results provide physically-motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows.
△ Less
Submitted 20 February, 2024;
originally announced February 2024.
-
Massive Star Cluster Formation I. High Star Formation Efficiency While Resolving Feedback of Individual Stars
Authors:
Brooke Polak,
Mordecai-Mark Mac Low,
Ralf S. Klessen,
Jia Wei Teh,
Claude Cournoyer-Cloutier,
Eric P. Andersson,
Sabrina M. Appel,
Aaron Tran,
Sean C. Lewis,
Maite J. C. Wilhelm,
Simon Portegies Zwart,
Simon C. O. Glover,
Long Wang,
Stephen L. W. McMillan
Abstract:
The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upg…
▽ More
The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgrade Torch by implementing the N-body code PeTar, thereby enabling Torch to handle massive clusters forming from $10^6\rm\, M_\odot$ clouds with $\ge10^5$ individual stars. We present results from Torch simulations of star clusters forming from $10^4, 10^5$, and $10^6\rm M_\odot$ turbulent, spherical gas clouds (named M4, M5, M6) of radius $R=11.7$ pc. We find that star formation is highly efficient and becomes more so at higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities $2.325\times 10^{1,2,3}\rm\, M_\odot\, pc^{-2}$, after a free-fall time of $t_{ff}=6.7,2.1,0.67$ Myr, we find that $\sim30\%$, 40%, and 60% of the cloud mass has formed into stars, respectively. The final integrated star formation efficiency is $32\%,\, 65\%$, and 85\% for M4, M5, and M6. Observations of nearby clusters similar to M4 have similar integrated star formation efficiencies of $\leq30\%$. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense ($\ge 10^2 \rm\,M_\odot\,pc^{-2}$) and massive ($\ge10^5\rm\, M_\odot$). In such environments, the dynamical time of the cloud becomes short enough that stellar feedback cannot act quickly enough to slow star formation.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Electron Heating in 2D Particle-in-Cell Simulations of Quasi-Perpendicular Low-Beta Shocks
Authors:
Aaron Tran,
Lorenzo Sironi
Abstract:
We measure the thermal electron energization in 1D and 2D particle-in-cell (PIC) simulations of quasi-perpendicular, low-beta ($β_p=0.25$) collisionless ion-electron shocks with mass ratio $m_i/m_e=200$, fast Mach number $\mathcal{M}_{ms}=1$-$4$, and upstream magnetic field angle $θ_{Bn} = 55$-$85^\circ$ from shock normal $\hat{\boldsymbol{n}}$. It is known that shock electron heating is described…
▽ More
We measure the thermal electron energization in 1D and 2D particle-in-cell (PIC) simulations of quasi-perpendicular, low-beta ($β_p=0.25$) collisionless ion-electron shocks with mass ratio $m_i/m_e=200$, fast Mach number $\mathcal{M}_{ms}=1$-$4$, and upstream magnetic field angle $θ_{Bn} = 55$-$85^\circ$ from shock normal $\hat{\boldsymbol{n}}$. It is known that shock electron heating is described by an ambipolar, $\boldsymbol{B}$-parallel electric potential jump, $Δφ_\parallel$, that scales roughly linearly with the electron temperature jump. Our simulations have $Δφ_\parallel/(0.5 m_i {u_\mathrm{sh}}^2) \sim 0.1$-$0.2$ in units of the pre-shock ions' bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measure $φ_\parallel$, including the use of de Hoffmann-Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate of $φ_\parallel$ in our low-$β_p$ shocks. We further focus on two $θ_{Bn}=65^\circ$ shocks: a $\mathcal{M}_s=4$ ($\mathcal{M}_A=1.8$) case with a long, $30 d_i$ precursor of whistler waves along $\hat{\boldsymbol{n}}$, and a $\mathcal{M}_s=7$ ($\mathcal{M}_A=3.2$) case with a shorter, $5d_i$ precursor of whistlers oblique to both $\hat{\boldsymbol{n}}$ and $\boldsymbol{B}$; $d_i$ is the ion skin depth. Within the precursors, $φ_\parallel$ has a secular rise towards the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of the $\mathcal{M}_s=4$, $θ_{Bn}=65^\circ$ case, $φ_\parallel$ shows a weak dependence on the electron plasma-to-cyclotron frequency ratio $ω_{pe}/Ω_{ce}$, and $φ_\parallel$ decreases by a factor of 2 as $m_i/m_e$ is raised to the true proton-electron value of 1836.
△ Less
Submitted 26 March, 2024; v1 submitted 31 August, 2023;
originally announced August 2023.
-
Electron-Ion Temperature Ratio in Astrophysical Shocks
Authors:
John C. Raymond,
Parviz Ghavamian,
Artem Bohdan,
Dongsu Ryu,
Jacek Niemiec,
Lorenzo Sironi,
Aaron Tran,
Elena Amato,
Masahiro Hoshino,
Martin Pohl,
Takanobu Amano,
Federico Fiuza
Abstract:
Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T$_e$/T$_p$ inferred from the H alpha profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Her…
▽ More
Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T$_e$/T$_p$ inferred from the H alpha profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Here we determine T$_e$/T$_p$ for supernova remnant shocks using H alpha profiles, shock speeds from proper motions, and electron temperatures from X-ray spectra. We also improve the estimates of sound speed and Alfvén speed used to determine Mach numbers. We find that the H alpha determinations are robust and that the discrepancies among supernova remnant shocks, solar wind shocks and computer-simulated shocks remain. We discuss some possible contributing factors, including shock precursors, turbulence and varying preshock conditions.
△ Less
Submitted 15 March, 2023;
originally announced March 2023.
-
Early Evolution and 3D Structure of Embedded Star Clusters
Authors:
Claude Cournoyer-Cloutier,
Alison Sills,
William E. Harris,
Sabrina M. Appel,
Sean C. Lewis,
Brooke Polak,
Aaron Tran,
Martijn J. C. Wilhelm,
Mordecai-Mark Mac Low,
Stephen L. W. McMillan,
Simon Portegies Zwart
Abstract:
We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at…
▽ More
We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at 0.0683 pc resolution for $\sim$ 2 Myr after the onset of star formation, with virial parameters $α$ = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial $α$ having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over $\sim$ 0.01 Myr timescales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses $\lesssim$ 1000 M$_{\odot}$ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation.
△ Less
Submitted 16 February, 2023;
originally announced February 2023.
-
Radiation shielding of protoplanetary discs in young star-forming regions
Authors:
Martijn J. C. Wilhelm,
Simon Portegies Zwart,
Claude Cournoyer-Cloutier,
Sean C. Lewis,
Brooke Polak,
Aaron Tran,
Mordecai-Mark Mac Low
Abstract:
Protoplanetary discs spend their lives in the dense environment of a star forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the AMUSE framework to couple the Torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two en…
▽ More
Protoplanetary discs spend their lives in the dense environment of a star forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the AMUSE framework to couple the Torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two environmental processes. We compare simulations with and without extinction of photoevaporation-driving radiation. We find that the majority of discs in our simulations are considerably shielded from photoevaporation-driving radiation for at least 0.5 Myr after the formation of the first massive stars. Radiation shielding increases disc lifetimes by an order of magnitude and can let a disc retain more solid material for planet formation. The reduction in external photoevaporation leaves discs larger and more easily dynamically truncated, although external photoevaporation remains the dominant mass loss process. Finally, we find that the correlation between disc mass and projected distance to the most massive nearby star (often interpreted as a sign of external photoevaporation) can be erased by the presence of less massive stars that dominate their local radiation field. Overall, we find that the presence and dynamics of gas in embedded clusters with massive stars is important for the evolution of protoplanetary discs.
△ Less
Submitted 7 February, 2023;
originally announced February 2023.
-
Early-Forming Massive Stars Suppress Star Formation and Hierarchical Cluster Assembly
Authors:
Sean C. Lewis,
Stephen L. W. McMillan,
Mordecai-Mark Mac Low,
Claude Cournoyer-Cloutier,
Brooke Polak,
Martijn J. C. Wilhelm,
Aaron Tran,
Alison Sills,
Simon Portegies Zwart,
Ralf S. Klessen,
Joshua E. Wall
Abstract:
Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation…
▽ More
Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suite \texttt{Torch}, combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer, $N$-body dynamics, and stellar feedback to model four initially identical $10^4$ M$_\odot$ giant molecular clouds with a Gaussian density profile peaking at $521.5 \mbox{ cm}^{-3}$. Using the \texttt{Torch} software suite through the \texttt{AMUSE} framework we modify three of the models to ensure that the first star that forms is very massive (50, 70, 100 M$_\odot$). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star forming region. The star formation rate is suppressed, reducing the total mass of stars formed. Our fiducial control model without an early massive star has a larger star formation rate and total efficiency by up to a factor of three and a higher average star formation efficiency per free-fall time by up to a factor of seven. Early-forming massive stars promote the buildup of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster.
△ Less
Submitted 28 February, 2023; v1 submitted 2 December, 2022;
originally announced December 2022.
-
Electron Re-acceleration via Ion Cyclotron Waves in the Intracluster Medium
Authors:
Aaron Tran,
Lorenzo Sironi,
Francisco Ley,
Ellen G. Zweibel,
Mario A. Riquelme
Abstract:
In galaxy clusters, the intracluster medium (ICM) is expected to host a diffuse, long-lived, and invisible population of "fossil" cosmic-ray electrons (CRe) with 1-100 MeV energies. These CRe, if re-accelerated by 100x in energy, can contribute synchrotron luminosity to cluster radio halos, relics, and phoenices. Re-acceleration may be aided by CRe scattering upon the ion-Larmor-scale waves that s…
▽ More
In galaxy clusters, the intracluster medium (ICM) is expected to host a diffuse, long-lived, and invisible population of "fossil" cosmic-ray electrons (CRe) with 1-100 MeV energies. These CRe, if re-accelerated by 100x in energy, can contribute synchrotron luminosity to cluster radio halos, relics, and phoenices. Re-acceleration may be aided by CRe scattering upon the ion-Larmor-scale waves that spawn when ICM is compressed, dilated, or sheared. We study CRe scattering and energy gain due to ion cyclotron (IC) waves generated by continuously-driven compression in 1D fully kinetic particle-in-cell simulations. We find that pitch-angle scattering of CRe by IC waves induces energy gain via magnetic pumping. In an optimal range of IC-resonant momenta, CRe may gain up to ~10-30% of their initial energy in one compress/dilate cycle with magnetic field amplification ~3-6x, assuming adiabatic decompression without further scattering and averaging over initial pitch angle.
△ Less
Submitted 24 February, 2023; v1 submitted 26 September, 2022;
originally announced September 2022.
-
A Heating Mechanism via Magnetic Pumping in the Intracluster Medium
Authors:
Francisco Ley,
Ellen G. Zweibel,
Mario Riquelme,
Lorenzo Sironi,
Drake Miller,
Aaron Tran
Abstract:
Turbulence driven by AGN activity, cluster mergers and galaxy motion constitutes an attractive energy source for heating the intracluster medium (ICM). How this energy dissipates into the ICM plasma remains unclear, given its low collisionality and high magnetization (precluding viscous heating by Coulomb processes). Kunz et al. 2011 proposed a viable heating mechanism based on the anisotropy of t…
▽ More
Turbulence driven by AGN activity, cluster mergers and galaxy motion constitutes an attractive energy source for heating the intracluster medium (ICM). How this energy dissipates into the ICM plasma remains unclear, given its low collisionality and high magnetization (precluding viscous heating by Coulomb processes). Kunz et al. 2011 proposed a viable heating mechanism based on the anisotropy of the plasma pressure (gyroviscous heating) under ICM conditions. The present paper builds upon that work and shows that particles can be gyroviscously heated by large-scale turbulent fluctuations via magnetic pumping. We study how the anisotropy evolves under a range of forcing frequencies, what waves and instabilities are generated and demonstrate that the particle distribution function acquires a high energy tail. For this, we perform particle-in-cell simulations where we periodically vary the mean magnetic field $\textbf{B}(t)$. When $\textbf{B}(t)$ grows (dwindles), a pressure anisotropy $P_{\perp}>P_{\parallel}$ ($P_{\perp}< P_{\parallel}$) builds up ($P_{\perp}$ and $P_{\parallel}$ are, respectively, the pressures perpendicular and parallel to $\textbf{B}(t)$). These pressure anisotropies excite mirror ($P_{\perp}>P_{\parallel}$) and oblique firehose ($P_{\parallel}>P_{\perp}$) instabilities, which trap and scatter the particles, limiting the anisotropy and providing a channel to heat the plasma. The efficiency of this mechanism depends on the frequency of the large-scale turbulent fluctuations and the efficiency of the scattering the instabilities provide in their nonlinear stage. We provide a simplified analytical heating model that captures the phenomenology involved. Our results show that this process can be relevant in dissipating and distributing turbulent energy at kinetic scales in the ICM.
△ Less
Submitted 31 August, 2022;
originally announced September 2022.
-
LIGO Detector Characterization in the Second and Third Observing Runs
Authors:
D. Davis,
J. S. Areeda,
B. K. Berger,
R. Bruntz,
A. Effler,
R. C. Essick,
R. P. Fisher,
P. Godwin,
E. Goetz,
A. F. Helmling-Cornell,
B. Hughey,
E. Katsavounidis,
A. P. Lundgren,
D. M. Macleod,
Z. Márka,
T. J. Massinger,
A. Matas,
J. McIver,
G. Mo,
K. Mogushi,
P. Nguyen,
L. K. Nuttall,
R. M. S. Schofield,
D. H. Shoemaker,
S. Soni
, et al. (262 additional authors not shown)
Abstract:
The characterization of the Advanced LIGO detectors in the second and third observing runs has increased the sensitivity of the instruments, allowing for a higher number of detectable gravitational-wave signals, and provided confirmation of all observed gravitational-wave events. In this work, we present the methods used to characterize the LIGO detectors and curate the publicly available datasets…
▽ More
The characterization of the Advanced LIGO detectors in the second and third observing runs has increased the sensitivity of the instruments, allowing for a higher number of detectable gravitational-wave signals, and provided confirmation of all observed gravitational-wave events. In this work, we present the methods used to characterize the LIGO detectors and curate the publicly available datasets, including the LIGO strain data and data quality products. We describe the essential role of these datasets in LIGO-Virgo Collaboration analyses of gravitational-waves from both transient and persistent sources and include details on the provenance of these datasets in order to support analyses of LIGO data by the broader community. Finally, we explain anticipated changes in the role of detector characterization and current efforts to prepare for the high rate of gravitational-wave alerts and events in future observing runs.
△ Less
Submitted 9 July, 2021; v1 submitted 27 January, 2021;
originally announced January 2021.
-
Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method
Authors:
Claude Cournoyer-Cloutier,
Aaron Tran,
Sean Lewis,
Joshua E. Wall,
William E. Harris,
Mordecai-Mark Mac Low,
Stephen L. W. McMillan,
Simon Portegies Zwart,
Alison Sills
Abstract:
The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled M…
▽ More
The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled MHD and direct N-body star cluster formation code Torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the $q \geq 0.1$ detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semi-major axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems.
△ Less
Submitted 11 November, 2020;
originally announced November 2020.
-
Population Properties of Compact Objects from the Second LIGO-Virgo Gravitational-Wave Transient Catalog
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva
, et al. (1316 additional authors not shown)
Abstract:
We report on the population of the 47 compact binary mergers detected with a false-alarm rate 1/yr in the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. We observe several characteristics of the merging binary black hole (BBH) population not discernible until now. First, we find that the primary mass spectrum contains structure beyond a power-law with a sharp high-mass cut-off; i…
▽ More
We report on the population of the 47 compact binary mergers detected with a false-alarm rate 1/yr in the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. We observe several characteristics of the merging binary black hole (BBH) population not discernible until now. First, we find that the primary mass spectrum contains structure beyond a power-law with a sharp high-mass cut-off; it is more consistent with a broken power law with a break at $39.7^{+20.3}_{-9.1}\,M_\odot$, or a power law with a Gaussian feature peaking at $33.1^{+4.0}_{-5.6}\,M_\odot$ (90\% credible interval). While the primary mass distribution must extend to $\sim65\,M_\odot$ or beyond, only $2.9^{+3.5}_{1.7}\%$ of systems have primary masses greater than $45\,M_\odot$. Second, we find that a fraction of BBH systems have component spins misaligned with the orbital angular momentum, giving rise to precession of the orbital plane. Moreover, 12% to 44% of BBH systems have spins tilted by more than $90^\circ$, giving rise to a negative effective inspiral spin parameter $χ_\mathrm{eff}$. Under the assumption that such systems can only be formed by dynamical interactions, we infer that between 25% and 93% of BBH with non-vanishing $|χ_\mathrm{eff}| > 0.01$ are dynamically assembled. Third, we estimate merger rates, finding $\mathcal{R}_\text{BBH} = 23.9^{+14.3}_{8.6}$ Gpc$^{-3}$ yr$^{-1}$ for BBH and $\mathcal{R}_\text{BNS}= 320^{+490}_{-240}$ Gpc$^{-3}$ yr$^{-1}$ for binary neutron stars. We find that the BBH rate likely increases with redshift ($85\%$ credibility), but not faster than the star-formation rate ($86\%$ credibility). Additionally, we examine recent exceptional events in the context of our population models, finding that the asymmetric masses of GW190412 and the high component masses of GW190521 are consistent with our models, but the low secondary mass of GW190814 makes it an outlier.
△ Less
Submitted 25 February, 2021; v1 submitted 27 October, 2020;
originally announced October 2020.
-
Tests of General Relativity with Binary Black Holes from the second LIGO-Virgo Gravitational-Wave Transient Catalog
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva
, et al. (1322 additional authors not shown)
Abstract:
Gravitational waves enable tests of general relativity in the highly dynamical and strong-field regime. Using events detected by LIGO-Virgo up to 1 October 2019, we evaluate the consistency of the data with predictions from the theory. We first establish that residuals from the best-fit waveform are consistent with detector noise, and that the low- and high-frequency parts of the signals are in ag…
▽ More
Gravitational waves enable tests of general relativity in the highly dynamical and strong-field regime. Using events detected by LIGO-Virgo up to 1 October 2019, we evaluate the consistency of the data with predictions from the theory. We first establish that residuals from the best-fit waveform are consistent with detector noise, and that the low- and high-frequency parts of the signals are in agreement. We then consider parametrized modifications to the waveform by varying post-Newtonian and phenomenological coefficients, improving past constraints by factors of ${\sim}2$; we also find consistency with Kerr black holes when we specifically target signatures of the spin-induced quadrupole moment. Looking for gravitational-wave dispersion, we tighten constraints on Lorentz-violating coefficients by a factor of ${\sim}2.6$ and bound the mass of the graviton to $m_g \leq 1.76 \times 10^{-23} \mathrm{eV}/c^2$ with 90% credibility. We also analyze the properties of the merger remnants by measuring ringdown frequencies and damping times, constraining fractional deviations away from the Kerr frequency to $δ\hat{f}_{220} = 0.03^{+0.38}_{-0.35}$ for the fundamental quadrupolar mode, and $δ\hat{f}_{221} = 0.04^{+0.27}_{-0.32}$ for the first overtone; additionally, we find no evidence for postmerger echoes. Finally, we determine that our data are consistent with tensorial polarizations through a template-independent method. When possible, we assess the validity of general relativity based on collections of events analyzed jointly. We find no evidence for new physics beyond general relativity, for black hole mimickers, or for any unaccounted systematics.
△ Less
Submitted 16 June, 2021; v1 submitted 27 October, 2020;
originally announced October 2020.
-
GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run
Authors:
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
S. Akcay,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva,
S. B. Anderson
, et al. (1327 additional authors not shown)
Abstract:
We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave event…
▽ More
We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave events. At this threshold, we expect a contamination fraction of less than 10%. Of these, 26 candidate events were reported previously in near real-time through GCN Notices and Circulars; 13 are reported here for the first time. The catalog contains events whose sources are black hole binary mergers up to a redshift of ~0.8, as well as events whose components could not be unambiguously identified as black holes or neutron stars. For the latter group, we are unable to determine the nature based on estimates of the component masses and spins from gravitational wave data alone. The range of candidate events which are unambiguously identified as binary black holes (both objects $\geq 3~M_\odot$) is increased compared to GWTC-1, with total masses from $\sim 14~M_\odot$ for GW190924_021846 to $\sim 150~M_\odot$ for GW190521. For the first time, this catalog includes binary systems with significantly asymmetric mass ratios, which had not been observed in data taken before April 2019. We also find that 11 of the 39 events detected since April 2019 have positive effective inspiral spins under our default prior (at 90% credibility), while none exhibit negative effective inspiral spin. Given the increased sensitivity of Advanced LIGO and Advanced Virgo, the detection of 39 candidate events in ~26 weeks of data (~1.5 per week) is consistent with GWTC-1.
△ Less
Submitted 8 March, 2021; v1 submitted 27 October, 2020;
originally announced October 2020.
-
Gravitational-wave constraints on the equatorial ellipticity of millisecond pulsars
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
G. Allen,
A. Allocca,
P. A. Altin,
A. Amato,
S. Anand,
A. Ananyeva
, et al. (1311 additional authors not shown)
Abstract:
We present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first t…
▽ More
We present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first time we are able to match (for PSR J0437-4715) or surpass (for PSR J0711-6830) the indirect limits on gravitational-wave emission from recycled pulsars inferred from their observed spin-downs, and constrain their equatorial ellipticities to be less than $10^{-8}$. For each of the five pulsars, we perform targeted searches that assume a tight coupling between the gravitational-wave and electromagnetic signal phase evolution. We also present constraints on PSR J0711-6830, the Crab pulsar and the Vela pulsar from a search that relaxes this assumption, allowing the gravitational-wave signal to vary from the electromagnetic expectation within a narrow band of frequencies and frequency derivatives.
△ Less
Submitted 13 October, 2020; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Electron Heating in Perpendicular Low-Beta Shocks
Authors:
Aaron Tran,
Lorenzo Sironi
Abstract:
Collisionless shocks heat electrons in the solar wind, interstellar blast waves, and hot gas permeating galaxy clusters. How much shock heating goes to electrons instead of ions, and what plasma physics controls electron heating? We simulate 2-D perpendicular shocks with a fully kinetic particle-in-cell code. For magnetosonic Mach number $\mathcal{M}_\mathrm{ms} \sim 1$-$10$ and plasma beta…
▽ More
Collisionless shocks heat electrons in the solar wind, interstellar blast waves, and hot gas permeating galaxy clusters. How much shock heating goes to electrons instead of ions, and what plasma physics controls electron heating? We simulate 2-D perpendicular shocks with a fully kinetic particle-in-cell code. For magnetosonic Mach number $\mathcal{M}_\mathrm{ms} \sim 1$-$10$ and plasma beta $β_\mathrm{p} \lesssim 4$, the post-shock electron/ion temperature ratio $T_\mathrm{e}/T_\mathrm{i}$ decreases from $1$ to $0.1$ with increasing $\mathcal{M}_\mathrm{ms}$. In a representative $\mathcal{M}_\mathrm{ms}=3.1$, $β_\mathrm{p}=0.25$ shock, electrons heat above adiabatic compression in two steps: ion-scale $E_\parallel = \vec{E} \cdot \hat{b}$ accelerates electrons into streams along $\vec{B}$, which then relax via two-stream-like instability. The $\vec{B}$-parallel heating is mostly induced by waves; $\vec{B}$-perpendicular heating is mostly adiabatic compression by quasi-static fields.
△ Less
Submitted 21 August, 2020; v1 submitted 25 February, 2020;
originally announced February 2020.
-
Energy Dependence of Synchrotron X-Ray Rims in Tycho's Supernova Remnant
Authors:
Aaron Tran,
Brian J. Williams,
Robert Petre,
Sean M. Ressler,
Stephen P. Reynolds
Abstract:
Several young supernova remnants exhibit thin X-ray bright rims of synchrotron radiation at their forward shocks. Thin rims require strong magnetic field amplification beyond simple shock compression if rim widths are only limited by electron energy losses. But, magnetic field damping behind the shock could produce similarly thin rims with less extreme field amplification. Variation of rim width w…
▽ More
Several young supernova remnants exhibit thin X-ray bright rims of synchrotron radiation at their forward shocks. Thin rims require strong magnetic field amplification beyond simple shock compression if rim widths are only limited by electron energy losses. But, magnetic field damping behind the shock could produce similarly thin rims with less extreme field amplification. Variation of rim width with energy may thus discriminate between competing influences on rim widths. We measured rim widths around Tycho's supernova remnant in 5 energy bands using an archival 750 ks Chandra observation. Rims narrow with increasing energy and are well described by either loss-limited or damped scenarios, so X-ray rim width-energy dependence does not uniquely specify a model. But, radio counterparts to thin rims are not loss-limited and better reflect magnetic field structure. Joint radio and X-ray modeling favors magnetic damping in Tycho's SNR with damping lengths ~1--5% of remnant radius and magnetic field strengths ~50--400 $μ$G assuming Bohm diffusion. X-ray rim widths are ~1% of remnant radius, somewhat smaller than inferred damping lengths. Electron energy losses are important in all models of X-ray rims, suggesting that the distinction between loss-limited and damped models is blurred in soft X-rays. All loss-limited and damping models require magnetic fields $\gtrsim$ 20 $μ$G, affirming the necessity of magnetic field amplification beyond simple compression.
△ Less
Submitted 2 September, 2015;
originally announced September 2015.
-
Fueling AGN II: Spatially Resolved Molecular Inflows and Outflows
Authors:
R. I. Davies,
W. Maciejewski,
E. K. S. Hicks,
E. Emsellem,
P. Erwin,
L. Burtscher,
G. Dumas,
M. Lin,
M. A. Malkan,
F. Mueller-Sanchez,
G. Orban de Xivry,
D. J. Rosario,
A. Schnorr-Mueller,
A. Tran
Abstract:
We analyse the 2-dimensional distribution and kinematics of the stars as well as molecular and ionised gas in the central few hundred parsecs of 5 active and 5 matched inactive galaxies. The equivalent widths of the Br-gamma line indicate there is no on-going star formation in their nuclei, although recent (terminated) starbursts are possible in the active galaxies. The stellar velocity fields sho…
▽ More
We analyse the 2-dimensional distribution and kinematics of the stars as well as molecular and ionised gas in the central few hundred parsecs of 5 active and 5 matched inactive galaxies. The equivalent widths of the Br-gamma line indicate there is no on-going star formation in their nuclei, although recent (terminated) starbursts are possible in the active galaxies. The stellar velocity fields show no signs of non-circular motions, while the 1-0S(1) H_2 kinematics exhibit significant deviations from simple circular rotation. In the active galaxies the H_2 kinematics reveal inflow and outflow superimposed on disk rotation. Steady-state circumnuclear inflow is seen in three AGN, and hydrodynamical models indicate it can be driven by a large scale bar. In three of the five AGN, molecular outflows are spatially resolved. The outflows are oriented such that they intersect, or have an edge close to, the disk - which may be the source of molecular gas in the outflow. The relatively low speeds imply the gas will fall back onto the disk; and with moderate outflow rates, they will have only a local impact on the host galaxy. H_2 was detected in two inactive galaxies. These exhibit chaotic circumnuclear dust morphologies and have molecular structures that are counter-rotating with respect to the main gas component, which could lead to gas inflow in the near future. In our sample, all four galaxies with chaotic dust morphology in the circumnuclear region exist in moderately dense groups with 10-15 members where accretion of stripped gas can easily occur.
△ Less
Submitted 9 July, 2014;
originally announced July 2014.
-
Electromgnetic-gravitational cross-sections in external elctromagnetic fields
Authors:
Hoang Ngoc Long,
Dang Van Soa,
Tuan A. Tran
Abstract:
The classical processes: the conversion of photons into gravitons in the static electromagnetic fields are considered by using Feynman perturbation techniques. The differential cross sections are presented for the conversion in the electric field of the flat condesor and the magnetic field of the selenoid. A numerical evaluation shows that the cross sections may have the observable value in the…
▽ More
The classical processes: the conversion of photons into gravitons in the static electromagnetic fields are considered by using Feynman perturbation techniques. The differential cross sections are presented for the conversion in the electric field of the flat condesor and the magnetic field of the selenoid. A numerical evaluation shows that the cross sections may have the observable value in the present technical scenario.
△ Less
Submitted 3 October, 1994;
originally announced October 1994.
