We report on the population of 47 compact binary mergers detected with a falsealarm rate of lt 1mathrmyr1 in the second LIGOVirgo GravitationalWave Transient Catalog. We observe several characteristics of the merging binary black hole BBH population not discernible until now. First the primary mass spectrum contains structure beyond a power law with a sharp highmass cutoff it is more consistent with a broken power law with a break at 39.79.120.3Modot  or a power law with a Gaussian feature peaking at 33.15.64.0Modot  90 credible interval. While the primary mass distribution must extend to sim 65Modot  or beyond only 2.91.73.5   of systems have primary masses greater than 45Modot  . 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  44  of BBH systems have spins tilted by more than 90 giving rise to a negative effective inspiral spin parameter chi mathrmeff . Under the assumption that such systems can only be formed by dynamical interactions we infer that between 25 and 93 of BBHs with nonvanishing  chi mathrmeff gt 0.01 are dynamically assembled. Third we estimate merger rates finding  mathcal R mathrmBBH23.98.614.3mathrmGpc3mathrmyr1 for BBHs and  mathcal R mathrmBNS320240490mathrmGpc3mathrmyr1 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.
['Gravitational waves', 'Astrophysical black holes', 'Stellar mass black holes', 'Massive stars', 'Compact objects', 'Bayesian statistics', '678', '98', '1611', '732', '288', '1900', 'Astrophysics - High Energy Astrophysical Phenomena', 'General Relativity and Quantum Cosmology']
Population Properties of Compact Objects from the Second LIGOVirgo GravitationalWave Transient Catalog
{'No Header': "Population properties of compact objects from the second LIGO-Virgo Gravitational-Wave Transient Catalog  \nR. Abbott, 1 T. D. Abbott, 2 S. Abraham, 3 F. Acernese, 4, 5 K. Ackley, 6 A. Adams, 7 C. Adams, 8 R. X. Adhikari, 1 V. B. Adya, 9 C. Affeldt, 10, 11 M. Agathos, 12, 13 K. Agatsuma, 14 N. Aggarwal, 15 O. D. Aguiar, 16 L. Aiello, 17, 18 A. Ain, 19, 20 P. Ajith, 21 G. Allen, 22 A. Allocca, 19 P. A. Altin, 9 A. Amato, 23 S. Anand, 1 A. Ananyeva, 1 S. B. Anderson, 1 W. G. Anderson, 24 S. V. Angelova, 25 S. Ansoldi, 26, 27 J. M. Antelis, 28 S. Antier, 29 S. Appert, 1 K. Arai, 1 M. C. Araya, 1 J. S. Areeda, 30 M. Arène, 29 N. Arnaud, 31, 32 S. M. Aronson, 33 K. G. Arun, 34 Y. Asali, 35 S. Ascenzi, 17, 36 G. Ashton, 6 S. M. Aston, 8 P. Astone, 37 F. Aubin, 38 P. Aufmuth, 10, 11 K. AultONeal, 28 C. Austin, 2 V. Avendano, 39 S. Babak, 29 F. Badaracco, 17, 18 M. K. M. Bader, 40 S. Bae, 41 A. M. Baer, 7 S. Bagnasco, 42 J. Baird, 29 M. Ball, 43 G. Ballardin, 32 S. W. Ballmer, 44 A. Bals, 28 A. Balsamo, 7 G. Baltus, 45 S. Banagiri, 46 D. Bankar, 3 R. S. Bankar, 3 J. C. Barayoga, 1 C. Barbieri, 47, 48, 49 B. C. Barish, 1 D. Barker, 50 P. Barneo, 51 S. Barnum, 52 F. Barone, 53, 5 B. Barr, 54 L. Barsotti, 52 M. Barsuglia, 29 D. Barta, 55 J. Bartlett, 50 I. Bartos, 33 R. Bassiri, 56 A. Basti, 20, 19 M. Bawaj, 57, 58 J. C. Bayley, 54 M. Bazzan, 59, 60 B. R. Becher, 61 B. Bécsy, 62 V. M. Bedakihale, 63 M. Bejger, 64 I. Belahcene, 31 D. Beniwal, 65 M. G. Benjamin, 28 T. F. Bennett, 66 J. D. Bentley, 14 F. Bergamin, 10, 11 B. K. Berger, 56 G. Bergmann, 10, 11 S. Bernuzzi, 13 C. P. L. Berry, 15 D. Bersanetti, 67 A. Bertolini, 40 J. Betzwieser, 8 R. Bhandare, 68 A. V. Bhandari, 3 D. Bhattacharjee, 69 J. Bidler, 30 I. A. Bilenko, 70 G. Billingsley, 1 R. Birney, 71 O. Birnholtz, 72 S. Biscans, 1, 52 M. Bischi, 73, 74 S. Biscoveanu, 52 A. Bisht, 10, 11 M. Bitossi, 32, 19 M.-A. Bizouard, 75 J. K. Blackburn, 1 J. Blackman, 76 C. D. Blair, 77 D. G. Blair, 77 R. M. Blair, 50 O. Blanch, 78 F. Bobba, 79, 80 N. Bode, 10, 11 M. Boer, 75 Y. Boetzel, 81 G. Bogaert, 75 M. Boldrini, 82, 37 F. Bondu, 83 E. Bonilla, 56 R. Bonnand, 38 P. Booker, 10, 11 B. A. Boom, 40 R. Bork, 1 V. Boschi, 19 S. Bose, 3 V. Bossilkov, 77 V. Boudart, 45 Y. Bouffanais, 59, 60 A. Bozzi, 32 C. Bradaschia, 19 P. R. Brady, 24 A. Bramley, 8 M. Branchesi, 17, 18 J. E. Brau, 43 M. Breschi, 13 T. Briant, 84 J. H. Briggs, 54 F. Brighenti, 73, 74 A. Brillet, 75 M. Brinkmann, 10, 11 P. Brockill, 24 A. F. Brooks, 1 J. Brooks, 32 D. D. Brown, 65 S. Brunett, 1 G. Bruno, 85 R. Bruntz, 7 A. Buikema, 52 T. Bulik, 86 H. J. Bulten, 40, 87 A. Buonanno, 88, 89 R. Buscicchio, 14 D. Buskulic, 38 R. L. Byer, 56 M. Cabero, 10, 11 L. Cadonati, 90 M. Caesar, 91 G. Cagnoli, 23 C. Cahillane, 1 J. Calderón Bustillo, 6 J. D. Callaghan, 54 T. A. Callister, 92 E. Calloni, 93, 5 J. B. Camp, 94 M. Canepa, 95, 67 K. C. Cannon, 96 H. Cao, 65 J. Cao, 97 G. Carapella, 79, 80 F. Carbognani, 32 M. F. Carney, 15 M. Carpinelli, 98, 99 G. Carullo, 20, 19 T. L. Carver, 100 J. Casanueva Diaz, 32 C. Casentini, 101, 36 S. Caudill, 40 M. Cavaglià, 69 F. Cavalier, 31 R. Cavalieri, 32 G. Cella, 19 P. Cerdá-Durán, 102 E. Cesarini, 36 W. Chaibi, 75 K. Chakravarti, 3 C.-L. Chan, 103 C. Chan, 96 K. Chandra, 104 P. Chanial, 32 S. Chao, 105 P. Charlton, 106 E. A. Chase, 15 E. Chassande-Mottin, 29 D. Chatterjee, 24 D. Chattopadhyay, 107 M. Chaturvedi, 68 K. Chatziioannou, 92 A. Chen, 103 H. Y. Chen, 108 X. Chen, 77 Y. Chen, 76 H.-P. Cheng, 33 C. K. Cheong, 103 H. Y. Chia, 33 F. Chiadini, 109, 80 R. Chierici, 110 A. Chincarini, 67 A. Chiummo, 32 G. Cho, 111 H. S. Cho, 112 M. Cho, 89 S. Choate, 91 N. Christensen, 75 Q. Chu, 77 S. Chua, 84 K. W. Chung, 113 S. Chung, 77 G. Ciani, 59, 60 P. Ciecielag, 64 M. Cieślar, 64 M. Cifaldi, 101, 36 A. A. Ciobanu, 65 R. Ciolfi, 114, 60 F. Cipriano, 75 A. Cirone, 95, 67 F. Clara, 50 E. N. Clark, 115 J. A. Clark, 90 L. Clarke, 116 P. Clearwater, 117 S. Clesse, 85 F. Cleva, 75 E. Coccia, 17, 18 P.-F. Cohadon, 84 D. E. Cohen, 31 M. Colleoni, 118 C. G. Collette, 119 C. Collins, 14 M. Colpi, 47, 48 M. Constancio Jr., 16 L. Conti, 60 S. J. Cooper, 14 P. Corban, 8 T. R. Corbitt, 2 I. Cordero-Carrión, 120 S. Corezzi, 58, 57 K. R. Corley, 35 N. Cornish, 62 D. Corre, 31 A. Corsi, 121 S. Cortese, 32 C. A. Costa, 16 R. Cotesta, 88 M. W. Coughlin, 46, 1 S. B. Coughlin, 15, 100 J.-P. Coulon, 75 S. T. Countryman, 35 P. Couvares, 1 P. B. Covas, 118 D. M. Coward, 77 M. J. Cowart, 8 D. C. Coyne, 1 R. Coyne, 122 J. D. E. Creighton, 24 T. D. Creighton, 123 M. Croquette, 84 S. G. Crowder, 124 J.R. Cudell, 45 T. J. Cullen, 2 A. Cumming, 54 R. Cummings, 54 L. Cunningham, 54 E. Cuoco, 32, 125 M. Curylo, 86 T. Dal Canton, 31, 88 G. Dálya, 126 A. Dana, 56 L. M. DaneshgaranBajastani, 66 B. D'Angelo, 95, 67 S. L. Danilishin, 127 S. D'Antonio, 36 K. Danzmann, 10, 11 C. Darsow-Fromm, 128 A. Dasgupta, 63 L. E. H. Datrier, 54 V. Dattilo, 32 I. Dave, 68 M. Davier, 31 G. S. Davies, 129 D. Davis, 1 E. J. Daw, 130 R. Dean, 91 D. DeBra, 56 M. Deenadayalan, 3 J. Degallaix, 131 M. De Laurentis, 93, 5 S. Deléglise, 84 V. Del Favero, 132 F. De Lillo, 85 N. De Lillo, 54 W. Del Pozzo, 20, 19 L. M. DeMarchi, 15 F. De Matteis, 101, 36 V. D'Emilio, 100 N. Demos, 52 T. Denker, 10, 11 T. Dent, 129 A. Depasse, 85 R. De Pietri, 133, 134 R. De Rosa, 93, 5 C. De Rossi, 32 R. DeSalvo, 135, 80 O. de Varona, 10, 11 S. Dhurandhar, 3 M. C. Díaz, 123 M. Diaz-Ortiz Jr., 33 N. A. Didio, 44 T. Dietrich, 40 L. Di Fiore, 5 C. DiFronzo, 14 C. Di Giorgio, 79, 80 F. Di Giovanni, 102 M. Di Giovanni, 136, 137 T. Di Girolamo, 93, 5 A. Di Lieto, 20, 19 B. Ding, 119 S. Di Pace, 82, 37 I. Di Palma, 82, 37 F. Di Renzo, 20, 19 A. K. Divakarla, 33 A. Dmitriev, 14 Z. Doctor, 43 L. D'Onofrio, 93, 5 F. Donovan, 52 K. L. Dooley, 100 S. Doravari, 3 I. Dorrington, 100 T. P. Downes, 24 M. Drago, 17, 18 J. C. Driggers, 50 Z. Du, 97 J.-G. Ducoin, 31 P. Dupej, 54 O. Durante, 79, 80 D. D'Urso, 98, 99 P.-A. Duverne, 31 S. E. Dwyer, 50 P. J. Easter, 6 G. Eddolls, 54 B. Edelman, 43 T. B. Edo, 130 O. Edy, 138 A. Effler, 8 J. Eichholz, 9 S. S. Eikenberry, 33 M. Eisenmann, 38 R. A. Eisenstein, 52 A. Ejlli, 100 L. Errico, 93, 5 R. C. Essick, 108 H. Estellés, 118 D. Estevez, 38 Z. B. Etienne, 139 T. Etzel, 1 M. Evans, 52  \nT. M. Evans, 8 B. E. Ewing, 140 V. Fafone, 101, 36, 17 H. Fair, 44 S. Fairhurst, 100 X. Fan, 97 A. M. Farah, 108 S. Farinon, 67 B. Farr, 43 W. M. Farr, 141, 92 E. J. Fauchon-Jones, 100 M. Favata, 39 M. Fays, 45, 130 M. Fazio, 142 J. Feicht, 1 M. M. Fejer, 56 F. Feng, 29 E. Fenyvesi, 55, 143 D. L. Ferguson, 90 A. Fernandez-Galiana, 52 I. Ferrante, 20, 19 T. A. Ferreira, 16 F. Fidecaro, 20, 19 P. Figura, 86 I. Fiori, 32 D. Fiorucci, 17, 18 M. Fishbach, 108 R. P. Fisher, 7 J. M. Fishner, 52 R. Fittipaldi, 144, 80 M. Fitz-Axen, 46 V. Fiumara, 145, 80 R. Flaminio, 38, 146 E. Floden, 46 E. Flynn, 30 H. Fong, 96 J. A. Font, 102, 147 P. W. F. Forsyth, 9 J.-D. Fournier, 75 S. Frasca, 82, 37 F. Frasconi, 19 Z. Frei, 126 A. Freise, 14 R. Frey, 43 V. Frey, 31 P. Fritschel, 52 V. V. Frolov, 8 G. G. Fronzé, 42 P. Fulda, 33 M. Fyffe, 8 H. A. Gabbard, 54 B. U. Gadre, 88 S. M. Gaebel, 14 J. R. Gair, 88 J. Gais, 103 S. Galaudage, 6 R. Gamba, 13 D. Ganapathy, 52 A. Ganguly, 21 S. G. Gaonkar, 3 B. Garaventa, 67, 95 C. García-Quirós, 118 F. Garufi, 93, 5 B. Gateley, 50 S. Gaudio, 28 V. Gayathri, 33 G. Gemme, 67 A. Gennai, 19 D. George, 22 J. George, 68 L. Gergely, 148 S. Ghonge, 90 Abhirup Ghosh, 88 Archisman Ghosh, 40, 149, 150, 151 S. Ghosh, 24, 39 B. Giacomazzo, 47, 48, 49 L. Giacoppo, 82, 37 J. A. Giaime, 2, 8 K. D. Giardina, 8 D. R. Gibson, 71 C. Gier, 25 K. Gill, 35 P. Giri, 19, 20 J. Glanzer, 2 A. E. Gleckl, 30 P. Godwin, 140 E. Goetz, 152 R. Goetz, 33 N. Gohlke, 10, 11 B. Goncharov, 6 G. González, 2 A. Gopakumar, 153 S. E. Gossan, 1 M. Gosselin, 20, 19 R. Gouaty, 38 B. Grace, 9 A. Grado, 154, 5 M. Granata, 131 V. Granata, 79 A. Grant, 54 S. Gras, 52 P. Grassia, 1 C. Gray, 50 R. Gray, 54 G. Greco, 73, 74 A. C. Green, 33 R. Green, 100 E. M. Gretarsson, 28 H. L. Griggs, 90 G. Grignani, 58, 57 A. Grimaldi, 136, 137 E. Grimes, 28 S. J. Grimm, 17, 18 H. Grote, 100 S. Grunewald, 88 P. Gruning, 31 J. G. Guerrero, 30 G. M. Guidi, 73, 74 A. R. Guimaraes, 2 G. Guixé, 51 H. K. Gulati, 63 Y. Guo, 40 Anchal Gupta, 1 Anuradha Gupta, 140 P. Gupta, 40, 155 E. K. Gustafson, 1 R. Gustafson, 156 F. Guzman, 115 L. Haegel, 29 O. Halim, 18, 17 E. D. Hall, 52 E. Z. Hamilton, 100 G. Hammond, 54 M. Haney, 81 M. M. Hanke, 10, 11 J. Hanks, 50 C. Hanna, 140 O. A. Hannuksela, 103 O. Hannuksela, 155, 40 H. Hansen, 50 T. J. Hansen, 28 J. Hanson, 8 T. Harder, 75 T. Hardwick, 2 K. Haris, 40, 155, 21 J. Harms, 17, 18 G. M. Harry, 157 I. W. Harry, 138 D. Hartwig, 128 R. K. Hasskew, 8 C.-J. Haster, 52 K. Haughian, 54 F. J. Hayes, 54 J. Healy, 132 A. Heidmann, 84 M. C. Heintze, 8 J. Heinze, 10, 11 J. Heinzel, 158 H. Heitmann, 75 F. Hellman, 159 P. Hello, 31 A. F. Helmling-Cornell, 43 G. Hemming, 32 M. Hendry, 54 I. S. Heng, 54 E. Hennes, 40 J. Hennig, 10, 11 M. H. Hennig, 10, 11 F. Hernandez Vivanco, 6 M. Heurs, 10, 11 S. Hild, 127 P. Hill, 25 A. S. Hines, 115 S. Hochheim, 10, 11 E. Hofgard, 56 D. Hofman, 131 J. N. Hohmann, 128 A. M. Holgado, 22 N. A. Holland, 9 I. J. Hollows, 130 Z. J. Holmes, 65 K. Holt, 8 D. E. Holz, 108 P. Hopkins, 100 C. Horst, 24 J. Hough, 54 E. J. Howell, 77 C. G. Hoy, 100 D. Hoyland, 14 Y. Huang, 52 M. T. Hübner, 6 A. D. Huddart, 116 E. A. Huerta, 22 B. Hughey, 28 V. Hui, 38 S. Husa, 118 S. H. Huttner, 54 B. M. Hutzler, 2 R. Huxford, 140 T. Huynh-Dinh, 8 B. Idzkowski, 86 A. Iess, 101, 36 S. Imperato, 15 H. Inchauspe, 33 C. Ingram, 65 G. Intini, 82, 37 M. Isi, 52 B. R. Iyer, 21 V. JaberianHamedan, 77 T. Jacqmin, 84 S. J. Jadhav, 160 S. P. Jadhav, 3 A. L. James, 100 K. Jani, 90 K. Janssens, 161 N. N. Janthalur, 160 P. Jaranowski, 162 D. Jariwala, 33 R. Jaume, 118 A. C. Jenkins, 113 M. Jeunon, 46 J. Jiang, 33 G. R. Johns, 7 A. W. Jones, 14 D. I. Jones, 163 J. D. Jones, 50 P. Jones, 14 R. Jones, 54 R. J. G. Jonker, 40 L. Ju, 77 J. Junker, 10, 11 C. V. Kalaghatgi, 100 V. Kalogera, 15 B. Kamai, 1 S. Kandhasamy, 3 G. Kang, 41 J. B. Kanner, 1 S. J. Kapadia, 21 D. P. Kapasi, 9 C. Karathanasis, 78 S. Karki, 69 R. Kashyap, 140 M. Kasprzack, 1 W. Kastaun, 10, 11 S. Katsanevas, 32 E. Katsavounidis, 52 W. Katzman, 8 K. Kawabe, 50 F. Kéfélian, 75 D. Keitel, 118 J. S. Key, 164 S. Khadka, 56 F. Y. Khalili, 70 I. Khan, 17, 36 S. Khan, 100 E. A. Khazanov, 165 N. Khetan, 17, 18 M. Khursheed, 68 N. Kijbunchoo, 9 C. Kim, 166 G. J. Kim, 90 J. C. Kim, 167 K. Kim, 168 W. S. Kim, 169 Y.-M. Kim, 170 C. Kimball, 15 P. J. King, 50 M. Kinley-Hanlon, 54 R. Kirchhoff, 10, 11 J. S. Kissel, 50 L. Kleybolte, 128 S. Klimenko, 33 T. D. Knowles, 139 E. Knyazev, 52 P. Koch, 10, 11 S. M. Koehlenbeck, 10, 11 G. Koekoek, 40, 171 S. Koley, 40 M. Kolstein, 78 K. Komori, 52 V. Kondrashov, 1 A. Kontos, 61 N. Koper, 10, 11 M. Korobko, 128 W. Z. Korth, 1 M. Kovalam, 77 D. B. Kozak, 1 C. Krämer, 10, 11 V. Kringel, 10, 11 N. V. Krishnendu, 10, 11 A. Królak, 172, 173 G. Kuehn, 10, 11 A. Kumar, 160 P. Kumar, 174 Rahul Kumar, 50 Rakesh Kumar, 63 K. Kuns, 52 S. Kwang, 24 B. D. Lackey, 88 D. Laghi, 20, 19 E. Lalande, 175 T. L. Lam, 103 A. Lamberts, 75, 176 M. Landry, 50 B. B. Lane, 52 R. N. Lang, 52 J. Lange, 132 B. Lantz, 56 R. K. Lanza, 52 I. La Rosa, 38 A. Lartaux-Vollard, 31 P. D. Lasky, 6 M. Laxen, 8 A. Lazzarini, 1 C. Lazzaro, 60, 59 P. Leaci, 82, 37 S. Leavey, 10, 11 Y. K. Lecoeuche, 50 H. M. Lee, 168 H. W. Lee, 167 J. Lee, 111 K. Lee, 56 J. Lehmann, 10, 11 E. Leon, 30 N. Leroy, 31 N. Letendre, 38 Y. Levin, 6 A. Li, 1 J. Li, 97 K. J. L. Li, 103 T. G. F. Li, 103 X. Li, 76 F. Linde, 177, 40 S. D. Linker, 66 J. N. Linley, 54 T. B. Littenberg, 178 J. Liu, 10, 11 X. Liu, 24 M. Llorens-Monteagudo, 102 R. K. L. Lo, 1 A. Lockwood, 179 L. T. London, 52 A. Longo, 180, 181 M. Lorenzini, 101, 36 V. Loriette, 182 M. Lormand, 8 G. Losurdo, 19 J. D. Lough, 10, 11 C. O. Lousto, 132 G. Lovelace, 30 H. Lück, 10, 11 D. Lumaca, 101, 36 A. P. Lundgren, 138 Y. Ma, 76 R. Macas, 100 M. MacInnis, 52 D. M. Macleod, 100 I. A. O. MacMillan, 1 A. Macquet, 75 I. Magaña Hernandez, 24 F. Magaña-Sandoval, 33 C. Magazzù, 19 R. M. Magee, 140 E. Majorana, 37 I. Maksimovic, 182 S. Maliakal, 1 A. Malik, 68 N. Man, 75 V. Mandic, 46 V. Mangano, 82, 37 G. L. Mansell, 50, 52 M. Manske, 24 M. Mantovani, 32 M. Mapelli, 59, 60 F. Marchesoni, 183, 57 F. Marion, 38 S. Márka, 35 Z. Márka, 35 C. Markakis, 12 A. S. Markosyan, 56 A. Markowitz, 1 E. Maros, 1 A. Marquina, 120 S. Marsat, 29 F. Martelli, 73, 74 I. W. Martin, 54 R. M. Martin, 39 M. Martinez, 78 V. Martinez, 23 D. V. Martynov, 14 H. Masalehdan, 128 K. Mason, 52 E. Massera, 130 A. Masserot, 38 T. J. Massinger, 52 M. Masso-Reid, 54 S. Mastrogiovanni, 29 A. Matas, 88 M. Mateu-Lucena, 118 F. Matichard, 1, 52 M. Matiushechkina, 10, 11 N. Mavalvala, 52 E. Maynard, 2 J. J. McCann, 77 R. McCarthy, 50 D. E. McClelland, 9 S. McCormick, 8 L. McCuller, 52 S. C. McGuire, 184 C. McIsaac, 138 J. McIver, 152 D. J. McManus, 9 T. McRae, 9 S. T. McWilliams, 139 D. Meacher, 24  \nG. D. Meadors, 6 M. Mehmet, 10, 11 A. K. Mehta, 88 A. Melatos, 117 D. A. Melchor, 30 G. Mendell, 50 A. Menendez-Vazquez, 78 R. A. Mercer, 24 L. Mereni, 131 K. Merfeld, 43 E. L. Merilh, 50 J. D. Merritt, 43 M. Merzougui, 75 S. Meshkov, 1 C. Messenger, 54 C. Messick, 185 R. Metzdorff, 84 P. M. Meyers, 117 F. Meylahn, 10, 11 A. Mhaske, 3 A. Miani, 136, 137 H. Miao, 14 I. Michaloliakos, 33 C. Michel, 131 H. Middleton, 117 L. Milano, 93, 5 A. L. Miller, 33, 85 S. Miller, 186, 92 M. Millhouse, 117 J. C. Mills, 100 E. Milotti, 187, 27 M. C. Milovich-Goff, 66 O. Minazzoli, 75, 188 Y. Minenkov, 36 Ll. M. Mir, 78 A. Mishkin, 33 C. Mishra, 189 T. Mistry, 130 S. Mitra, 3 V. P. Mitrofanov, 70 G. Mitselmakher, 33 R. Mittleman, 52 G. Mo, 52 K. Mogushi, 69 S. R. P. Mohapatra, 52 S. R. Mohite, 24 I. Molina, 30 M. Molina-Ruiz, 159 M. Mondin, 66 M. Montani, 73, 74 C. J. Moore, 14 D. Moraru, 50 F. Morawski, 64 G. Moreno, 50 S. Morisaki, 96 B. Mours, 190 C. M. Mow-Lowry, 14 S. Mozzon, 138 F. Muciaccia, 82, 37 Arunava Mukherjee, 54 D. Mukherjee, 140 Soma Mukherjee, 123 Subroto Mukherjee, 63 N. Mukund, 10, 11 A. Mullavey, 8 J. Munch, 65 E. A. Muñiz, 44 P. G. Murray, 54 S. L. Nadji, 10, 11 A. Nagar, 191, 42, 192 I. Nardecchia, 101, 36 L. Naticchioni, 37 R. K. Nayak, 193 B. F. Neil, 77 J. Neilson, 135, 80 G. Nelemans, 194 T. J. N. Nelson, 8 M. Nery, 10, 11 A. Neunzert, 164 K. Y. Ng, 52 S. Ng, 65 C. Nguyen, 29 P. Nguyen, 43 T. Nguyen, 52 S. A. Nichols, 2 S. Nissanke, 149, 40 F. Nocera, 32 M. Noh, 152 C. North, 100 D. Nothard, 195 L. K. Nuttall, 138 J. Oberling, 50 B. D. O'Brien, 33 J. O'Dell, 116 G. Oganesyan, 17, 18 G. H. Ogin, 196 J. J. Oh, 169 S. H. Oh, 169 F. Ohme, 10, 11 H. Ohta, 96 M. A. Okada, 16 C. Olivetto, 32 P. Oppermann, 10, 11 R. J. Oram, 8 B. O'Reilly, 8 R. G. Ormiston, 46 N. Ormsby, 7 L. F. Ortega, 33 R. O'Shaughnessy, 132 S. Ossokine, 88 C. Osthelder, 1 D. J. Ottaway, 65 H. Overmier, 8 B. J. Owen, 121 A. E. Pace, 140 G. Pagano, 20, 19 M. A. Page, 77 G. Pagliaroli, 17, 18 A. Pai, 104 S. A. Pai, 68 J. R. Palamos, 43 O. Palashov, 165 C. Palomba, 37 H. Pan, 105 P. K. Panda, 160 T. H. Pang, 40, 155 C. Pankow, 15 F. Pannarale, 82, 37 B. C. Pant, 68 F. Paoletti, 19 A. Paoli, 32 A. Paolone, 37, 197 W. Parker, 8, 184 D. Pascucci, 40 A. Pasqualetti, 32 R. Passaquieti, 20, 19 D. Passuello, 19 M. Patel, 7 B. Patricelli, 20, 19 E. Payne, 6 T. C. Pechsiri, 33 M. Pedraza, 1 M. Pegoraro, 60 A. Pele, 8 S. Penn, 198 A. Perego, 136, 137 C. J. Perez, 50 C. Périgois, 38 A. Perreca, 136, 137 S. Perriès, 110 J. Petermann, 128 D. Petterson, 1 H. P. Pfeiffer, 88 K. A. Pham, 46 K. S. Phukon, 40, 177, 3 O. J. Piccinni, 82, 37 M. Pichot, 75 M. Piendibene, 20, 19 F. Piergiovanni, 73, 74 L. Pierini, 82, 37 V. Pierro, 135, 80 G. Pillant, 32 F. Pilo, 19 L. Pinard, 131 I. M. Pinto, 135, 80, 191 K. Piotrzkowski, 85 M. Pirello, 50 M. Pitkin, 199 E. Placidi, 82 W. Plastino, 180, 181 C. Pluchar, 115 R. Poggiani, 20, 19 E. Polini, 38 D. Y. T. Pong, 103 S. Ponrathnam, 3 P. Popolizio, 32 E. K. Porter, 29 A. Poverman, 61 J. Powell, 107 M. Pracchia, 38 A. K. Prajapati, 63 K. Prasai, 56 R. Prasanna, 160 G. Pratten, 14 T. Prestegard, 24 M. Principe, 135, 191, 80 G. A. Prodi, 200, 137 L. Prokhorov, 14 P. Prosposito, 101, 36 A. Puecher, 40, 155 M. Punturo, 57 F. Puosi, 19, 20 P. Puppo, 37 M. Pürrer, 88 H. Qi, 100 V. Quetschke, 123 P. J. Quinonez, 28 R. Quitzow-James, 69 F. J. Raab, 50 G. Raaijmakers, 149, 40 H. Radkins, 50 N. Radulesco, 75 P. Raffai, 126 H. Rafferty, 201 S. X. Rail, 175 S. Raja, 68 C. Rajan, 68 B. Rajbhandari, 121 M. Rakhmanov, 123 K. E. Ramirez, 123 T. D. Ramirez, 30 A. Ramos-Buades, 118 J. Rana, 140 K. Rao, 15 P. Rapagnani, 82, 37 U. D. Rapol, 202 B. Ratto, 28 V. Raymond, 100 M. Razzano, 20, 19 J. Read, 30 T. Regimbau, 38 L. Rei, 67 S. Reid, 25 D. H. Reitze, 1, 33 P. Rettegno, 203, 42 F. Ricci, 82, 37 C. J. Richardson, 28 J. W. Richardson, 1 L. Richardson, 115 P. M. Ricker, 22 G. Riemenschneider, 203, 42 K. Riles, 156 M. Rizzo, 15 N. A. Robertson, 1, 54 F. Robinet, 31 A. Rocchi, 36 J. A. Rocha, 30 S. Rodriguez, 30 R. D. Rodriguez-Soto, 28 L. Rolland, 38 J. G. Rollins, 1 V. J. Roma, 43 M. Romanelli, 83 R. Romano, 4, 5 C. L. Romel, 50 A. Romero, 78 I. M. Romero-Shaw, 6 J. H. Romie, 8 S. Ronchini, 17, 18 C. A. Rose, 24 D. Rose, 30 K. Rose, 195 M. J. B. Rosell, 185 D. Rosińska, 86 S. G. Rosofsky, 22 M. P. Ross, 179 S. Rowan, 54 S. J. Rowlinson, 14 Santosh Roy, 3 Soumen Roy, 204 P. Ruggi, 32 K. Ryan, 50 S. Sachdev, 140 T. Sadecki, 50 M. Sakellariadou, 113 O. S. Salafia, 49, 48, 47 L. Salconi, 32 M. Saleem, 34 A. Samajdar, 40, 155 E. J. Sanchez, 1 J. H. Sanchez, 30 L. E. Sanchez, 1 N. Sanchis-Gual, 205 J. R. Sanders, 206 K. A. Santiago, 39 E. Santos, 75 T. R. Saravanan, 3 N. Sarin, 6 B. Sassolas, 131 B. S. Sathyaprakash, 140, 100 O. Sauter, 38 R. L. Savage, 50 V. Savant, 3 D. Sawant, 104 S. Sayah, 131 D. Schaetzl, 1 P. Schale, 43 M. Scheel, 76 J. Scheuer, 15 A. Schindler-Tyka, 33 P. Schmidt, 14 R. Schnabel, 128 R. M. S. Schofield, 43 A. Schönbeck, 128 E. Schreiber, 10, 11 B. W. Schulte, 10, 11 B. F. Schutz, 100, 10 O. Schwarm, 196 E. Schwartz, 100 J. Scott, 54 S. M. Scott, 9 M. Seglar-Arroyo, 38 E. Seidel, 22 D. Sellers, 8 A. S. Sengupta, 204 N. Sennett, 88 D. Sentenac, 32 V. Sequino, 93, 5 A. Sergeev, 165 Y. Setyawati, 10, 11 T. Shaffer, 50 M. S. Shahriar, 15 S. Sharifi, 2 A. Sharma, 17, 18 P. Sharma, 68 P. Shawhan, 89 H. Shen, 22 M. Shikauchi, 96 R. Shink, 175 D. H. Shoemaker, 52 D. M. Shoemaker, 90 K. Shukla, 159 S. ShyamSundar, 68 M. Sieniawska, 64 D. Sigg, 50 L. P. Singer, 94 D. Singh, 140 N. Singh, 86 A. Singha, 127 A. Singhal, 17, 37 A. M. Sintes, 118 V. Sipala, 98, 99 V. Skliris, 100 B. J. J. Slagmolen, 9 T. J. Slaven-Blair, 77 J. Smetana, 14 J. R. Smith, 30 R. J. E. Smith, 6 S. N. Somala, 207 E. J. Son, 169 S. Soni, 2 B. Sorazu, 54 V. Sordini, 110 F. Sorrentino, 67 N. Sorrentino, 20, 19 R. Soulard, 75 T. Souradeep, 202, 3 E. Sowell, 121 A. P. Spencer, 54 M. Spera, 59, 60, 15 A. K. Srivastava, 63 V. Srivastava, 44 K. Staats, 15 C. Stachie, 75 D. A. Steer, 29 M. Steinke, 10, 11 J. Steinlechner, 127, 54 S. Steinlechner, 127 D. Steinmeyer, 10, 11 S. P. Stevenson, 107 G. Stolle-McAllister, 195 D. J. Stops, 14 M. Stover, 195 K. A. Strain, 54 G. Stratta, 208, 74 A. Strunk, 50 R. Sturani, 209 A. L. Stuver, 91 J. Südbeck, 128 S. Sudhagar, 3 V. Sudhir, 52 H. G. Suh, 24 T. Z. Summerscales, 210 H. Sun, 77 L. Sun, 1 S. Sunil, 63 A. Sur, 64 J. Suresh, 96 P. J. Sutton, 100 B. L. Swinkels, 40 M. J. Szczepańczyk, 33 M. Tacca, 40 S. C. Tait, 54 C. Talbot, 6 A. J. Tanasijczuk, 85 D. B. Tanner, 33 D. Tao, 1 A. Tapia, 30 E. N. Tapia San Martin, 40 J. D. Tasson, 158 R. Taylor, 1 R. Tenorio, 118 L. Terkowski, 128 M. P. Thirugnanasambandam, 3 L. Thomas, 14 M. Thomas, 8 P. Thomas, 50 J. E. Thompson, 100 S. R. Thondapu, 68  \nK. A. Thorne, 8 E. Thrane, 6 Shubhanshu Tiwari, 81 Srishti Tiwari, 153 V. Tiwari, 100 K. Toland, 54 A. E. Tolley, 138 M. Tonelli, 20, 19 Z. Tornasi, 54 A. Torres-Forné, 88 C. I. Torrie, 1 I. Tosta e Melo, 98, 99 D. Töyrä, 9 A. T. Tran, 124 A. Trapananti, 183, 57 F. Travasso, 57, 183 G. Traylor, 8 M. C. Tringali, 86 A. Tripathee, 156 A. Trovato, 29 R. J. Trudeau, 1 D. S. Tsai, 105 K. W. Tsang, 40, 211, 155 M. Tse, 52 R. Tso, 76 L. Tsukada, 96 D. Tsuna, 96 T. Tsutsui, 96 M. Turconi, 75 A. S. Ubhi, 14 R. P. Udall, 90 K. Ueno, 96 D. Ugolini, 201 C. S. Unnikrishnan, 153 A. L. Urban, 2 S. A. Usman, 108 A. C. Utina, 127 H. Vahlbruch, 10, 11 G. Vajente, 1 A. Vajpeyi, 6 G. Valdes, 2 M. Valentini, 136, 137 V. Valsan, 24 N. van Bakel, 40 M. van Beuzekom, 40 J. F. J. van den Brand, 171, 87, 40 C. Van Den Broeck, 155, 40 D. C. Vander-Hyde, 44 L. van der Schaaf, 40 J. V. van Heijningen, 77 M. Vardaro, 177, 40 A. F. Vargas, 117 V. Varma, 76 S. Vass, 1 M. Vasúth, 55 A. Vecchio, 14 G. Vedovato, 60 J. Veitch, 54 P. J. Veitch, 65 K. Venkateswara, 179 J. Venneberg, 10, 11 G. Venugopalan, 1 D. Verkindt, 38 Y. Verma, 68 D. Veske, 35 F. Vetrano, 73 A. Viceré, 73, 74 A. D. Viets, 212 V. Villa-Ortega, 129 J.-Y. Vinet, 75 S. Vitale, 52 T. Vo, 44 H. Vocca, 58, 57 C. Vorvick, 50 S. P. Vyatchanin, 70 A. R. Wade, 9 L. E. Wade, 195 M. Wade, 195 R. C. Walet, 40 M. Walker, 7 G. S. Wallace, 25 L. Wallace, 1 S. Walsh, 24 J. Z. Wang, 156 S. Wang, 22 W. H. Wang, 123 Y. F. Wang, 103 R. L. Ward, 9 J. Warner, 50 M. Was, 38 N. Y. Washington, 1 J. Watchi, 119 B. Weaver, 50 L. Wei, 10, 11 M. Weinert, 10, 11 A. J. Weinstein, 1 R. Weiss, 52 F. Wellmann, 10, 11 L. Wen, 77 P. Weßels, 10, 11 J. W. Westhouse, 28 K. Wette, 9 J. T. Whelan, 132 D. D. White, 30 L. V. White, 44 B. F. Whiting, 33 C. Whittle, 52 D. M. Wilken, 10, 11 D. Williams, 54 M. J. Williams, 54 A. R. Williamson, 138 J. L. Willis, 1 B. Willke, 10, 11 D. J. Wilson, 115 M. H. Wimmer, 10, 11 W. Winkler, 10, 11 C. C. Wipf, 1 G. Woan, 54 J. Woehler, 10, 11 J. K. Wofford, 132 I. C. F. Wong, 103 J. Wrangel, 10, 11 J. L. Wright, 54 D. S. Wu, 10, 11 D. M. Wysocki, 132 L. Xiao, 1 H. Yamamoto, 1 L. Yang, 142 Y. Yang, 33 Z. Yang, 46 M. J. Yap, 9 D. W. Yeeles, 100 A. Yoon, 7 Hang Yu, 76 Haocun Yu, 52 S. H. R. Yuen, 103 A. Zadrożny, 173 M. Zanolin, 28 T. Zelenova, 32 J.-P. Zendri, 60 M. Zevin, 15 J. Zhang, 77 L. Zhang, 1 R. Zhang, 33 T. Zhang, 14 C. Zhao, 77 G. Zhao, 119 M. Zhou, 15 Z. Zhou, 15 X. J. Zhu, 6 A. B. Zimmerman, 185 M. E. Zucker, 1, 52 1  \nand J. Zweizig  \nThe LIGO Scientific Collaboration and the Virgo Collaboration  \n1 LIGO, California Institute of Technology, Pasadena, CA 91125, USA  \n2 Louisiana State University, Baton Rouge, LA 70803, USA  \n3 Inter-University Centre for Astronomy and Astrophysics, Pune 411007, India  \n4 Dipartimento di Farmacia, Università di Salerno, I-84084 Fisciano, Salerno, Italy  \n5 INFN, Sezione di Napoli, Complesso Universitario di Monte S.Angelo, I-80126 Napoli, Italy  \n6 OzGrav, School of Physics & Astronomy, Monash University, Clayton 3800, Victoria, Australia  \n7 Christopher Newport University, Newport News, VA 23606, USA  \n8 LIGO Livingston Observatory, Livingston, LA 70754, USA  \n9 OzGrav, Australian National University, Canberra, Australian Capital Territory 0200, Australia  \n10 Max Planck Institute for Gravitational Physics (Albert Einstein Institute), D-30167 Hannover, Germany  \n11 Leibniz Universität Hannover, D-30167 Hannover, Germany  \n12 University of Cambridge, Cambridge CB2 1TN, United Kingdom  \n13 Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universität Jena, D-07743 Jena, Germany  \n14 University of Birmingham, Birmingham B15 2TT, United Kingdom  \n15 Center for Interdisciplinary Exploration & Research in Astrophysics (CIERA), Northwestern University, Evanston, IL 60208, USA  \n16 Instituto Nacional de Pesquisas Espaciais, 12227-010 São José dos Campos, São Paulo, Brazil  \n17 Gran Sasso Science Institute (GSSI), I-67100 L'Aquila, Italy  \n18 INFN, Laboratori Nazionali del Gran Sasso, I-67100 Assergi, Italy  \n19 INFN, Sezione di Pisa, I-56127 Pisa, Italy  \n20 Università di Pisa, I-56127 Pisa, Italy  \n21 International Centre for Theoretical Sciences, Tata Institute of Fundamental Research, Bengaluru 560089, India  \n22  \nNCSA, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA  \n23 Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France  \n24 University of Wisconsin-Milwaukee, Milwaukee, WI 53201, USA  \n25 SUPA, University of Strathclyde, Glasgow G1 1XQ, United Kingdom  \n26 Dipartimento di Matematica e Informatica, Università di Udine, I-33100 Udine, Italy  \n27 INFN, Sezione di Trieste, I-34127 Trieste, Italy  \n28 Embry-Riddle Aeronautical University, Prescott, AZ 86301, USA  \n29 Université de Paris, CNRS, Astroparticule et Cosmologie, F-75013 Paris, France  \n30 California State University Fullerton, Fullerton, CA 92831, USA  \n- 31 Université Paris-Saclay, CNRS/IN2P3, IJCLab, 91405 Orsay, France  \n32 European Gravitational Observatory (EGO), I-56021 Cascina, Pisa, Italy  \n53  \n84  \n34  \n33 University of Florida, Gainesville, FL 32611, USA  \nChennai Mathematical Institute, Chennai 603103, India  \n35 Columbia University, New York, NY 10027, USA  \n36 INFN, Sezione di Roma Tor Vergata, I-00133 Roma, Italy  \n37 INFN, Sezione di Roma, I-00185 Roma, Italy  \n38 Laboratoire d'Annecy de Physique des Particules (LAPP), Univ. Grenoble Alpes, Université Savoie Mont Blanc, CNRS/IN2P3, F-74941  \nAnnecy, France  \n39 Montclair State University, Montclair, NJ 07043, USA  \n40 Nikhef, Science Park 105, 1098 XG Amsterdam, Netherlands  \n41 Korea Institute of Science and Technology Information, Daejeon 34141, South Korea  \n42 INFN Sezione di Torino, I-10125 Torino, Italy  \n43 University of Oregon, Eugene, OR 97403, USA  \n44  \nSyracuse University, Syracuse, NY 13244, USA  \n45 Université de Liège, B-4000 Liège, Belgium  \n46 University of Minnesota, Minneapolis, MN 55455, USA  \n47 Università degli Studi di Milano-Bicocca, I-20126 Milano, Italy  \n48 INFN, Sezione di Milano-Bicocca, I-20126 Milano, Italy  \n49 INAF, Osservatorio Astronomico di Brera sede di Merate, I-23807 Merate, Lecco, Italy  \n50 LIGO Hanford Observatory, Richland, WA 99352, USA  \n51 Institut de Ciències del Cosmos, Universitat de Barcelona, C/ Martí i Franquès 1, Barcelona, 08028, Spain  \n52 LIGO, Massachusetts Institute of Technology, Cambridge, MA 02139, USA  \nDipartimento di Medicina, Chirurgia e Odontoiatria 'Scuola Medica Salernitana,' Università di Salerno, I-84081 Baronissi, Salerno, Italy  \n54 SUPA, University of Glasgow, Glasgow G12 8QQ, United Kingdom  \n55 Wigner RCP, RMKI, H-1121 Budapest, Konkoly Thege Miklós út 29-33, Hungary  \n56 Stanford University, Stanford, CA 94305, USA  \n57 INFN, Sezione di Perugia, I-06123 Perugia, Italy  \n58 Università di Perugia, I-06123 Perugia, Italy  \n59 Università di Padova, Dipartimento di Fisica e Astronomia, I-35131 Padova, Italy  \n60 INFN, Sezione di Padova, I-35131 Padova, Italy  \n61 Bard College, 30 Campus Rd, Annandale-On-Hudson, NY 12504, USA  \n62 Montana State University, Bozeman, MT 59717, USA  \n63 Institute for Plasma Research, Bhat, Gandhinagar 382428, India  \n64 Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, 00-716, Warsaw, Poland  \n65 OzGrav, University of Adelaide, Adelaide, South Australia 5005, Australia  \n66 California State University, Los Angeles, 5151 State University Dr, Los Angeles, CA 90032, USA  \n67 INFN, Sezione di Genova, I-16146 Genova, Italy  \n68 RRCAT, Indore, Madhya Pradesh 452013, India  \n69 Missouri University of Science and Technology, Rolla, MO 65409, USA  \n70 Faculty of Physics, Lomonosov Moscow State University, Moscow 119991, Russia  \n71 SUPA, University of the West of Scotland, Paisley PA1 2BE, United Kingdom  \n72 Bar-Ilan University, Ramat Gan, 5290002, Israel  \n73 Università degli Studi di Urbino 'Carlo Bo', I-61029 Urbino, Italy  \n74 INFN, Sezione di Firenze, I-50019 Sesto Fiorentino, Firenze, Italy  \n75 Artemis, Université Côte d'Azur, Observatoire Côte d'Azur, CNRS, F-06304 Nice, France  \n76 Caltech CaRT, Pasadena, CA 91125, USA  \n77 OzGrav, University of Western Australia, Crawley, Western Australia 6009, Australia  \n78 Institut de Física d'Altes Energies (IFAE), Barcelona Institute of Science and Technology, and ICREA, E-08193 Barcelona, Spain  \n79 Dipartimento di Fisica 'E.R. Caianiello,' Università di Salerno, I-84084 Fisciano, Salerno, Italy  \n80 INFN, Sezione di Napoli, Gruppo Collegato di Salerno, Complesso Universitario di Monte S. Angelo, I-80126 Napoli, Italy  \n81 Physik-Institut, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland  \n82 Università di Roma 'La Sapienza', I-00185 Roma, Italy  \n83 Univ Rennes, CNRS, Institut FOTON - UMR6082, F-3500 Rennes, France  \nLaboratoire Kastler Brossel, Sorbonne Université, CNRS, ENS-Université PSL, Collège de France, F-75005 Paris, France  \n85 Université catholique de Louvain, B-1348 Louvain-la-Neuve, Belgium  \n86 Astronomical Observatory Warsaw University, 00-478 Warsaw, Poland  \n87 VU University Amsterdam, 1081 HV Amsterdam, Netherlands  \n88 Max Planck Institute for Gravitational Physics (Albert Einstein Institute), D-14476 Potsdam-Golm, Germany  \n89 University of Maryland, College Park, MD 20742, USA  \n90 School of Physics, Georgia Institute of Technology, Atlanta, GA 30332, USA  \n91 Villanova University, 800 Lancaster Ave, Villanova, PA 19085, USA  \n92 Center for Computational Astrophysics, Flatiron Institute, New York, NY 10010, USA  \n93 Università di Napoli 'Federico II', Complesso Universitario di Monte S.Angelo, I-80126 Napoli, Italy  \n94 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA  \n95 Dipartimento di Fisica, Università degli Studi di Genova, I-16146 Genova, Italy  \n96 RESCEU, University of Tokyo, Tokyo, 113-0033, Japan.  \n97 Tsinghua University, Beijing 100084, China  \n98 Università degli Studi di Sassari, I-07100 Sassari, Italy  \n99 INFN, Laboratori Nazionali del Sud, I-95125 Catania, Italy  \n100 Gravity Exploration Institute, Cardiff University, Cardiff CF24 3AA, United Kingdom  \n101 Università di Roma Tor Vergata, I-00133 Roma, Italy  \n102 Departamento de Astronomía y Astrofísica, Universitat de València, E-46100 Burjassot, València, Spain  \n103 The Chinese University of Hong Kong, Shatin, NT, Hong Kong  \n104 Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India  \n105 National Tsing Hua University, Hsinchu City, 30013 Taiwan, Republic of China  \n106 Charles Sturt University, Wagga Wagga, New South Wales 2678, Australia  \n107 OzGrav, Swinburne University of Technology, Hawthorn VIC 3122, Australia  \n108 University of Chicago, Chicago, IL 60637, USA  \n109 Dipartimento di Ingegneria Industriale (DIIN), Università di Salerno, I-84084 Fisciano, Salerno, Italy  \n110 Institut de Physique des 2 Infinis de Lyon, CNRS/IN2P3, Université de Lyon, Université Claude Bernard Lyon 1, F-69622 Villeurbanne,", 'France': "111 Seoul National University, Seoul 08826, South Korea  \n112 Pusan National University, Busan 46241, South Korea  \n113 King's College London, University of London, London WC2R 2LS, United Kingdom  \n114 INAF, Osservatorio Astronomico di Padova, I-35122 Padova, Italy  \n115 University of Arizona, Tucson, AZ 85721, USA  \n116 Rutherford Appleton Laboratory, Didcot OX11 0DE, United Kingdom  \n117 OzGrav, University of Melbourne, Parkville, Victoria 3010, Australia  \n118 Universitat de les Illes Balears, IAC3-IEEC, E-07122 Palma de Mallorca, Spain  \n119 Université Libre de Bruxelles, Brussels 1050, Belgium  \n120 Departamento de Matemáticas, Universitat de València, E-46100 Burjassot, València, Spain  \n121 Texas Tech University, Lubbock, TX 79409, USA  \n122 University of Rhode Island, Kingston, RI 02881, USA  \n123 The University of Texas Rio Grande Valley, Brownsville, TX 78520, USA  \n124 Bellevue College, Bellevue, WA 98007, USA  \n125 Scuola Normale Superiore, Piazza dei Cavalieri, 7 - 56126 Pisa, Italy  \n126 MTA-ELTE Astrophysics Research Group, Institute of Physics, Eötvös University, Budapest 1117, Hungary  \n127 Maastricht University, 6200 MD, Maastricht, Netherlands  \n128 Universität Hamburg, D-22761 Hamburg, Germany  \n129 IGFAE, Campus Sur, Universidade de Santiago de Compostela, 15782 Spain  \n130 The University of Sheffield, Sheffield S10 2TN, United Kingdom  \n131 Laboratoire des Matériaux Avancés (LMA), Institut de Physique des 2 Infinis de Lyon, CNRS/IN2P3, Université de Lyon, F-69622 Villeurbanne, France  \n132 Rochester Institute of Technology, Rochester, NY 14623, USA  \n133 Dipartimento di Scienze Matematiche, Fisiche e Informatiche, Università di Parma, I-43124 Parma, Italy  \n134 INFN, Sezione di Milano Bicocca, Gruppo Collegato di Parma, I-43124 Parma, Italy  \n135 Dipartimento di Ingegneria, Università del Sannio, I-82100 Benevento, Italy  \n136 Università di Trento, Dipartimento di Fisica, I-38123 Povo, Trento, Italy  \n137 INFN, Trento Institute for Fundamental Physics and Applications, I-38123 Povo, Trento, Italy  \n138 University of Portsmouth, Portsmouth, PO1 3FX, United Kingdom  \n139 West Virginia University, Morgantown, WV 26506, USA  \n140 The Pennsylvania State University, University Park, PA 16802, USA  \n141 Stony Brook University, Stony Brook, NY 11794, USA  \n142 Colorado State University, Fort Collins, CO 80523, USA  \n143 Institute for Nuclear Research, Hungarian Academy of Sciences, Bem t'er 18/c, H-4026 Debrecen, Hungary  \n194  \n144 CNR-SPIN, c/o Università di Salerno, I-84084 Fisciano, Salerno, Italy  \n145 Scuola di Ingegneria, Università della Basilicata, I-85100 Potenza, Italy  \n146 National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan  \n147 Observatori Astronòmic, Universitat de València, E-46980 Paterna, València, Spain  \n148  \nUniversity of Szeged, Dóm tér 9, Szeged 6720, Hungary  \n149 GRAPPA, Anton Pannekoek Institute for Astronomy and Institute for High-Energy Physics, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, Netherlands  \n150 Delta Institute for Theoretical Physics, Science Park 904, 1090 GL Amsterdam, Netherlands  \n151 Lorentz Institute, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, Netherlands  \n152 University of British Columbia, Vancouver, BC V6T 1Z4, Canada  \n153 Tata Institute of Fundamental Research, Mumbai 400005, India  \n154 INAF, Osservatorio Astronomico di Capodimonte, I-80131 Napoli, Italy  \n155 Department of Physics, Utrecht University, Princetonplein 1, 3584 CC Utrecht, Netherlands  \n156 University of Michigan, Ann Arbor, MI 48109, USA  \n157 American University, Washington, D.C. 20016, USA  \n158 Carleton College, Northfield, MN 55057, USA  \n159 University of California, Berkeley, CA 94720, USA  \n160 Directorate of Construction, Services & Estate Management, Mumbai 400094 India  \n161 Universiteit Antwerpen, Prinsstraat 13, 2000 Antwerpen, Belgium  \n162 University of Bialystok, 15-424 Bialystok, Poland  \n163 University of Southampton, Southampton SO17 1BJ, United Kingdom  \n164 University of Washington Bothell, Bothell, WA 98011, USA  \n165 Institute of Applied Physics, Nizhny Novgorod, 603950, Russia  \n166 Ewha Womans University, Seoul 03760, South Korea  \n167 Inje University Gimhae, South Gyeongsang 50834, South Korea  \n168 Korea Astronomy and Space Science Institute, Daejeon 34055, South Korea  \n169 National Institute for Mathematical Sciences, Daejeon 34047, South Korea  \n170 Ulsan National Institute of Science and Technology, Ulsan 44919, South Korea  \n171 Maastricht University, P.O. Box 616, 6200 MD Maastricht, Netherlands  \n172 Institute of Mathematics, Polish Academy of Sciences, 00656 Warsaw, Poland  \n173 National Center for Nuclear Research, 05-400 Świerk-Otwock, Poland  \n174 Cornell University, Ithaca, NY 14850, USA  \n175 Université de Montréal/Polytechnique, Montreal, Quebec H3T 1J4, Canada  \n176 Laboratoire Lagrange, Université Côte d'Azur, Observatoire Côte d'Azur, CNRS, F-06304 Nice, France  \n177 Institute for High-Energy Physics, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, Netherlands  \n178  \nNASA Marshall Space Flight Center, Huntsville, AL 35811, USA  \n179 University of Washington, Seattle, WA 98195, USA  \n180 Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, I-00146 Roma, Italy  \n181 INFN, Sezione di Roma Tre, I-00146 Roma, Italy  \n182 ESPCI, CNRS, F-75005 Paris, France  \n183 Università di Camerino, Dipartimento di Fisica, I-62032 Camerino, Italy  \n184 Southern University and A&M College, Baton Rouge, LA 70813, USA  \n185 Department of Physics, University of Texas, Austin, TX 78712, USA  \n186 Department of Physics, Smith College, Northampton, MA 01063, USA  \n187 Dipartimento di Fisica, Università di Trieste, I-34127 Trieste, Italy  \n188 Centre Scientifique de Monaco, 8 quai Antoine Ier, MC-98000, Monaco  \n189 Indian Institute of Technology Madras, Chennai 600036, India  \n190 Institut Pluridisciplinaire Hubert CURIEN, 23 rue du loess - BP28 67037 Strasbourg cedex 2, France  \n191 Museo Storico della Fisica e Centro Studi e Ricerche 'Enrico Fermi', I-00184 Roma, Italy  \n192  \nInstitut des Hautes Etudes Scientifiques, F-91440 Bures-sur-Yvette, France  \n193 Indian Institute of Science Education and Research, Kolkata, Mohanpur, West Bengal 741252, India  \nDepartment of Astrophysics/IMAPP, Radboud University Nijmegen, P.O. Box 9010, 6500 GL Nijmegen, Netherlands  \n195 Kenyon College, Gambier, OH 43022, USA  \n196 Whitman College, 345 Boyer Avenue, Walla Walla, WA 99362 USA  \n197 Consiglio Nazionale delle Ricerche - Istituto dei Sistemi Complessi, Piazzale Aldo Moro 5, I-00185 Roma, Italy  \n198 Hobart and William Smith Colleges, Geneva, NY 14456, USA  \n199 Lancaster University, Lancaster LA1 4YW, United Kingdom  \n200 Università di Trento, Dipartimento di Matematica, I-38123 Povo, Trento, Italy  \n201 Trinity University, San Antonio, TX 78212, USA  \n202 Indian Institute of Science Education and Research, Pune, Maharashtra 411008, India  \n203 Dipartimento di Fisica, Università degli Studi di Torino, I-10125 Torino, Italy  \n204 Indian Institute of Technology, Palaj, Gandhinagar, Gujarat 382355, India  \n205 Centro de Astrofísica e Gravitação (CENTRA), Departamento de Física, Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal  \n206 Marquette University, 11420 W. Clybourn St., Milwaukee, WI 53233, USA  \n207 Indian Institute of Technology Hyderabad, Sangareddy, Khandi, Telangana 502285, India  \n208 INAF, Osservatorio di Astrofisica e Scienza dello Spazio, I-40129 Bologna, Italy  \n209 International Institute of Physics, Universidade Federal do Rio Grande do Norte, Natal RN 59078-970, Brazil  \n210 Andrews University, Berrien Springs, MI 49104, USA  \n211 Van Swinderen Institute for Particle Physics and Gravity, University of Groningen, Nijenborgh 4, 9747 AG Groningen, Netherlands 212 Concordia University Wisconsin, Mequon, WI 53097, USA", 'ABSTRACT': 'We report on the population of the 47 compact binary mergers detected with a false-alarm rate < 1yr -1 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, 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 glyph[circledot] , or a power law with a Gaussian feature peaking at 33 . 1 +4 . 0 -5 . 6 M glyph[circledot] (90% credible interval). While the primary mass distribution must extend to ∼ 65 M glyph[circledot] or beyond, only 2 . 9 +3 . 5 -1 . 7 % of systems have primary masses greater than 45 M glyph[circledot] . 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 · , giving rise to a negative effective inspiral spin parameter χ 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 | χ eff | > 0 . 01 are dynamically assembled. Third, we estimate merger rates, finding R BBH = 23 . 9 +14 . 3 -8 . 6 Gpc -3 yr -1 for BBH and R 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.  \nKeywords: gravitational waves', '1. INTRODUCTION': "We analyze the population properties of black holes (BHs) and neutron stars (NSs) in compact binary systems using data from the LIGO-Virgo Gravitational-wave Transient Catalog 2 (GWTC-2; Abbott et al. 2020c). The GWTC-2 catalog combines observations from the first two observing runs (O1 and O2; Abbott et al. 2019b) and the first half of the third observing run (O3a; Abbott et al. 2020c) of the Advanced LIGO (Aasi et al. 2015) and Advanced Virgo (Acernese et al. 2015) gravitational-wave observatories. With the 39 additional candidates from O3a, we have more than quadrupled the number of events from O1 and O2, published in the first LIGO-Virgo Transient Catalog (GWTC-1; Abbott et al. 2019b). Counting only events with false alarm rate (FAR) < 1 yr -1 (as opposed to the less conservative FAR threshold of < 2 yr -1  \nin GWTC-2), the new combined catalog includes: two binary NS (BNS) events, 44 confident binary black hole (BBH) events, and one NS-BH (NSBH) candidate, which may be a BBH-a topic we discuss below. We define BBH events as systems where both masses are above 3 M glyph[circledot] at 90% credibility. These 47 events are listed in Table 1. Our chosen FAR threshold ensures a relatively pure sample with only ∼ 1 noise event (see Section 2) and excludes three marginal triggers presented in GWTC2. Two of these excluded events are BBH candidates (GW190719\\_215514 and GW190909\\_114149) and one is an NSBH candidate (GW190426\\_152155).  \nThe latest observations include a number of individually remarkable events, which invite theoretical speculation while providing challenges to existing models. The observation of high-mass binaries such as  \nGW190521 (Abbott et al. 2020e), which has a primary mass m 1 ∼ 85 M glyph[circledot] , is in tension with the sharp mass cutoff m max = 40 . 8 +11 . 8 -4 . 4 M glyph[circledot] inferred from the GWTC1 detections, forcing us to reconsider models for the distribution of black hole (BH) masses in binary systems (Abbott et al. 2020e,f). Here and in the following, the primary mass m 1 refers to the bigger of the two component masses in the binary, while the secondary mass m 2 refers to the smaller of the two. Along with GW190521, GW190602\\_175927 and GW190519\\_153544 also have primary masses above 45 M glyph[circledot] at > 99% credibility. These high-mass binaries are interesting from a theoretical perspective since they occur in the predicted pair-instability gap (Barkat et al. 1967; Fowler & Hoyle 1964; Heger & Woosley 2002; Heger et al. 2003; Woosley & Heger 2015; Belczynski et al. 2016). Additionally, GWTC-2 includes the first systems with confidently asymmetric component masses, including GW190412 with mass ratio m 2 /m 1 ≡ q = 0 . 28 +0 . 12 -0 . 06 (Abbott et al. 2020a) and GW190814 (Abbott et al. 2020b), with q = 0 . 112 +0 . 008 -0 . 009 . Furthermore, the secondary mass of GW190814, m 2 = 2 . 59 +0 . 08 -0 . 09 M glyph[circledot] , is near the purported NS-BH gap (Bailyn et al. 1998; Özel et al. 2011; Farr et al. 2011), posing a challenge to our understanding of binary formation (Abbott et al. 2020b; Zevin et al. 2020; Mandel et al. 2021). We can gain insight into these exceptional events by studying them in the context of the larger population of compact binaries (Fishbach et al. 2020b).  \nIn particular, studying the enlarged population of BBH events enables us to investigate how compact binaries form. 1 Several evolutionary channels have been proposed to explain the origin of the compact binary mergers observed with Advanced LIGO and Advanced Virgo. The isolated binary channel may occur via common envelope evolution (e.g., Bethe & Brown 1998; Portegies Zwart & Yungelson 1998; Belczynski et al. 2002; Dominik et al. 2015), the remnants of Population III stars (e.g., Madau & Rees 2001; Inayoshi et al. 2017), or chemically homogeneous stellar evolution (e.g., Marchant et al. 2016; de Mink & Mandel 2016; Mandel & de Mink 2016). Evolution via common envelope predicts BBH systems with component masses up to ∼ 50 M glyph[circledot] , mass ratios in the range 0 . 3 glyph[lessorsimilar] q < 1 , and nearly aligned spins (Kalogera 2000; Mandel & O'Shaughnessy 2010; Dominik et al. 2013; Giacobbo et al. 2017; Eldridge et al. 2017; Olejak et al. 2020).  \nIn the chemically homogeneous scenario, BBH systems are expected to form with nearly equal mass components, in addition to aligned spins and component masses in the range ∼ 20 -50 M glyph[circledot] . In isolated formation scenarios, component BHs form via stellar collapse, and are thus not expected to occur within the pair-instability mass gap, ∼ 50 -120 M glyph[circledot] , but may populate either side of the gap.  \nAlternatively, BBH mergers could form via dynamical interactions in young stellar clusters, globular clusters, or nuclear star clusters (e.g., Kulkarni et al. 1993; Sigurdsson & Hernquist 1993; Portegies Zwart & McMillan 2000), triple systems (e.g. Antonini et al. 2014; Kimpson et al. 2016; Antonini et al. 2017; Vigna-Gómez et al. 2020) or the disks of active galactic nuclei (e.g. McKernan et al. 2012; Bartos et al. 2017; Stone et al. 2017; Fragione et al. 2019). Dynamical formation in dense stellar clusters typically produces an isotropic distribution of spin directions (e.g. Vitale et al. 2017b; Rodriguez et al. 2016), and may enable hierarchical mergers characterized by relatively high-mass binaries (e.g. Antonini & Rasio 2016; Mapelli 2016; McKernan et al. 2018; Rodriguez et al. 2018; Arca Sedda et al. 2020) and large BH spins (Fishbach et al. 2017; Gerosa & Berti 2017). Finally, BBH systems might originate as part of a primordial BH population in the early Universe (Carr & Hawking 1974; Carr et al. 2016). Most primordial BH models predict low spins and isotropic spin orientation (Fernandez & Profumo 2019).  \nBefore we continue, we summarize key questions addressed in the previous analysis of GWTC-1 data (Abbott et al. 2019a) and how the inclusion of O3a events affects our findings:  \n- 1. Are there BBH systems with component masses glyph[greaterorsimilar] 45 M glyph[circledot] ? Following O1 and O2, we inferred that 99% of BBH systems have primary mass below m 99% ≈ 45 M glyph[circledot] . Moreover, this limit was consistent with a sharp cutoff at m max = 40 . 8 +11 . 8 -4 . 4 M glyph[circledot] , hypothesized to correspond to the lower edge of the pair-instability mass gap expected from supernova theory (Woosley & Heger 2015; Heger & Woosley 2002; Heger et al. 2003; Fishbach & Holz 2017; Talbot & Thrane 2018). While the GWTC-2 events remain consistent with 97 . 1 +1 . 7 -3 . 5 % of BBH systems having primary masses below 45 M glyph[circledot] (inferred using the the Power Law + Peak mass model described in Section 3, or 'Model C' from Abbott et al. 2019a), high-mass detections such as GW190521, GW190602\\_175927 and GW190519\\_153544 imply a non-zero rate of BBH mergers beyond the ∼ 45 M glyph[circledot] mass limit. We infer that the merger rate for systems with primary masses in the range  \n- 45 M glyph[circledot] < m 1 < 100 M glyph[circledot] is 0 . 70 +0 . 65 -0 . 35 Gpc -3 yr -1 , consistent with estimates inferred from GWTC1 (Fishbach et al. 2020b).\n- 2. Is there a preference for nearly equal component masses? All of the GWTC-1 observations are consistent with equal-mass binaries, with mass ratios q ≡ m 2 /m 1 = 1 . In O3a, however, we detected two events with mass ratios bounded confidently away from unity: GW190814 and GW190412, though, most binaries are consistent with equal-mass. The NSBH candidate GW190426\\_152155, if real, also has unequal component masses q = 0 . 25 +0 . 41 -0 . 15 , but is above the FAR threshold for this work.\n- 3. Does the merger rate evolve with redshift? From GWTC-1 we inferred that the BBH merger rate is 53 . 2 +55 . 8 -28 . 2 Gpc -3 yr -1 , assuming a rate density that is uniform in comoving volume. Allowing for a rate that evolves with redshift (Fishbach et al. 2018), we found that the local merger rate is R BBH ( z = 0) = 19 . 7 +57 . 3 -15 . 9 Gpc -3 yr -1 , and that, while still consistent with a uniformin-comoving volume model, the rate density is likely increasing with redshift with 93% credibility (Abbott et al. 2019a). With GWTC-2, we are able to more tightly bound the BBH merger rate at R BBH = 23 . 9 +14 . 3 -8 . 6 Gpc -3 yr -1 (assuming the Power Law + Peak mass model and a constantin-comoving-volume merger rate), as well as its evolution with redshift. The data remain consistent with a merger rate that does not evolve with redshift, but prefer a rate that increases with redshift ( 85 % credibility). Using the Power Law redshift evolution model of Section 3, we find that the merger rate evolves slower than the naive expectation of (1 + z ) 2 . 7 from the local ( z glyph[lessorsimilar] 1 ) star formation rate (SFR; Madau & Dickinson 2014) at 86 % credibility.\n- 4. How fast do black holes spin? From GWTC-1, we inferred that the BH spin component aligned with the orbital angular momentum is typically small (Abbott et al. 2016a; Farr et al. 2017; Farr et al. 2018; Tiwari et al. 2018; Abbott et al. 2019a; Wysocki et al. 2019a; Fernandez & Profumo 2019; Roulet & Zaldarriaga 2019; Miller et al. 2020; Roulet et al. 2020). Among the GWTC-1 events, GW151226 is the only event to exhibit unambiguous signs of spin (Abbott et al. 2016b; Vitale et al. 2017a; Kimball et al. 2020a), while a few other events, including GW170729, show a mild preference for spin (Chatziioannou et al. 2019). We were unable to determine if this typical smallness was  \nbecause the spin vectors are misaligned, because the spin magnitudes are small, or both, although the GWTC-1 data weakly disfavors the scenario in which all spins are perfectly aligned (Farr et al. 2017; Tiwari et al. 2018; Abbott et al. 2019a; Biscoveanu et al. 2020). With additional data, we can now say confidently that some BBH systems have spins misaligned with the orbital angular momentum. A nonzero fraction of systems have measurable in-plane spin components, leading them to undergo precession of the orbital plane. Additionally, 12 % to 44 % of BBH systems merge with a negative effective inspiral spin parameter χ eff , see Eq. (5) below, implying that some component spins are tilted by more than 90 · relative to the orbital angular momentum axis. We refer to spins tilted more than 90 · as anti-aligned spins . 2 While some events identified in O3a have individually measurable nonzero spin, they occur infrequently, consistent with our previous estimates. We identify nine out of 44 BBH events in GWTC-2 with a positive effective inspiral spin parameter that excludes zero at greater than 95% credibility. 3 These nine events include both massive BBH systems like GW190519\\_153544, with a source primary mass m 1 = 66 . 0 +10 . 7 -12 . 0 M glyph[circledot] (90% credibility, uniform in redshifted mass prior) and less massive BBH systems like GW190728\\_064510, with m 1 = 12 . 3 +7 . 2 -2 . 2 M glyph[circledot] . No individual BBH events are observed with confidently negative effective inspiral spin parameter.  \n- 5. What is the minimum black hole mass? Using GWTC-1, we previously inferred that, if there is a low-mass cut-off in the BBH mass spectrum, it is somewhere below 9 M glyph[circledot] (Abbott et al. 2019a). With GWTC-2, we tighten the constraints on the minimum BH mass in a BBH system, finding m min < 6 . 6 M glyph[circledot] (90% credibility). Furthermore, we find that if the BH mass spectrum extends down to 3 M glyph[circledot] , it likely turns over at ∼ 7 . 8 +1 . 8 -2 . 0 M glyph[circledot] . This suggests that there may be a dearth of systems between NS and BH masses (Fishbach et al. 2020a). However, the O3a observation of GW190814 (Abbott et al. 2020b), with a sec-\n- 2 Our definition of anti-aligned does not imply perfect antialignment (tilt angles of exactly 180 · ) as the phrase is sometimes used to mean in the context of waveform modelling or numerical relativity.  \nondary mass m 2 = 2 . 59 +0 . 08 -0 . 09 M glyph[circledot] , complicates this picture. If the secondary mass is a BH, it would indicate that the BH mass spectrum extends below 3 M glyph[circledot] , to much lower masses than exhibited by the Galactic X-ray binary population (Bailyn et al. 1998; Özel et al. 2011; Farr et al. 2011, but see also Kreidberg et al. 2012; Thompson et al. 2019; Mandel et al. 2021). Alternatively, the secondary object in GW190814 could be an NS , but it would likely have to be significantly spinning to satisfy constraints on the maximum NS mass (Abbott et al. 2020b; Most et al. 2020; Tan et al. 2020; Essick & Landry 2020; Tews et al. 2020; Zhang & Li 2020). Either way, we find that GW190814 is an outlier with respect to the BBH population and the models we consider in this work. Unless stated otherwise, when presenting results on the BBH population, we exclude GW190814.  \nThe remainder of this paper is organized as follows. In Section 2 we describe the data used in this study and detail how we select events for analysis. In Section 3, we provide a high-level overview of our models for the distributions of binary mass, spins, and redshift. In Section 4, we describe the hierarchical method used to fit population models to data. In Section 5, we present key results and discuss the astrophysical implications. Concluding remarks are provided in Section 6. The Appendix provides additional details regarding the statistical method (Appendix A) and descriptions of models (Appendix B, D, E), outlier analyses and model checking (Appendix C), and other supplementary material, including a discussion of gravitational lensing (Appendix F). The data release for this paper is available online in Abbott et al. (2020d).", '2. DATA AND EVENT SELECTION': 'Searches for gravitational wave transients in the O3a data identified 39 candidate events with FAR below 2yr -1 (Abbott et al. 2020c). This FAR cut excludes two BBH candidates and one low-significance NSBH candidate. It is unlikely the results of our BBH analyses would differ qualitatively with the inclusion of these two additional BBH events because they are typical of other, more confident GWTC-2 detections, and including them would lead to a modest 5% increase in the catalog size. The event list was collated by combining results from several pipelines searching for compact binary mergers using archival data. The search pipelines include GstLAL (Sachdev et al. 2019; Hanna et al. 2020; Messick et al. 2017) and PyCBC (Nitz et al. 2019a; Allen et al. 2012; Allen 2005; Canton et al. 2014; Usman et al. 2016; Nitz et al. 2017), which use template-based matched  \nfiltering techniques, and cWB (Klimenko & Mitselmakher 2004; Klimenko et al. 2016), which uses a wavelet-based search that does not assume a physically parameterized signal model. These pipelines, along with two additional pipelines, MBTA (Adams et al. 2016), and SPIIR (Chu 2017), recovered most of the event candidates in lowlatency.  \nFor the population studies presented here, the event list is further restricted to the 36 events with FAR < 1 yr -1 as a means to increase the purity of the sample. This FAR threshold excludes the lower-significance triggers GW190426\\_152155 (a potential NSBH or BBH candidate), GW190719\\_215514 and GW190909\\_114149 that appear in Abbott et al. (2020c). At this FAR threshold, we expect ∼ 1 noise event, and therefore a contamination fraction of glyph[lessorsimilar] 3% . 4 In this work we assume that all of the event candidates that meet this FAR threshold are of astrophysical origin. For the population analysis of the GWTC-1 BBH events, the selection criteria used for inclusion in the study is the FAR and the probability p astro that the source is of astrophysical origin. This value was only computed for BBHs with FAR < 2yr -1 in Abbott et al. (2020c), so we simplify our criterion to only select on FAR. The set of GWTC-1 events that pass this FAR threshold is identical to the previous set chosen by FAR and p astro . Therefore, while the selection criteria here are different from our GWTC-1 analysis, the analyzed events are consistent.  \nIn addition to the 36 events from O3a which passed the FAR threshold, the 11 detections presented in GWTC1 (Abbott et al. 2019b) are also included in the event list used here to infer properties of the underlying population. All 47 events used in this analysis are tabulated in Table 1. Of these systems, 44 have both masses confidently in the BH mass range, with m 1 ≥ m 2 > 3 M glyph[circledot] . Unless stated otherwise, we restrict BBH population results to these 44 (see also Appendix C.3). Since our statistical framework relies on accurately quantifying the selection effects of our search, we only include events identified in GWTC-2, for which we have measured the search sensitivity; see Appendix A. In particular, the event list does not include candidates identified by independent analyses (Zackay et al. 2020; Venumadhav et al. 2020, 2019; Zackay et al. 2019; Nitz et al. 2019b, 2020; Magee et al. 2019) of the publicly released LIGO and Virgo data (Abbott et al. 2019c; Trovato et al. 2020). Galaudage et al. (2019) and Roulet et al. (2020) suggest that our results are unlikely to change significantly with the inclusion of these events, and in the future it may be possible to include events  \nfrom independent groups using a unified framework for the calculation of significance (Ashton & Thrane 2020; Pratten & Vecchio 2020; Roulet et al. 2020).  \nParameter estimation results for each candidate event (Abbott et al. 2020c) were obtained using the LALInference (Veitch et al. 2015; LIGO Scientific Collaboration 2018), RIFT (Lange et al. 2018; Wysocki et al. 2019b), and Bilby (Ashton et al. 2019; Romero-Shaw et al. 2020b) pipelines, the latter employing the dynesty nested sampling tool (Speagle 2020). The parameter estimation pipelines use Bayesian sampling methods to produce fair draws from the posterior distribution function of the source parameters, conditioned on the data and given a models for the signal and noise (Abbott et al. 2016c).  \nFor the BBH events previously published in GWTC-1, we use the published posterior samples inferred using the IMRPhenomPv2 (Hannam et al. 2014; Khan et al. 2016; Husa et al. 2016; Bohé et al. 2016) and SEOBNRv3 (Pan et al. 2014; Taracchini et al. 2014; Babak et al. 2017) waveform models. For the new BBH events of GWTC-2, we use waveform approximants that include higher-order multipole content, including the IMRPhenomPv3HM (Khan et al. 2020), NRSur7dq4 (Varma et al. 2019) and SEOBNRv4PHM (Bohé et al. 2017; Ossokine et al. 2020). For all events, we average over these waveform families, in contrast to our previously published parameter estimation results for GW190521, which highlighted results from one waveform, NRSur7dq4 (Abbott et al. 2020e). A more complete description of the parameter estimation methods and waveform models used can be found in Section V of Abbott et al. (2020c).', '3. POPULATION MODELS': 'In this section, we provide a high-level overview of the different population models used in this paper. Each model is given a nickname indicated in boldface. There are three categories of models: models for the shape of the mass distribution (Section 3.1), models for the spin distribution (Section 3.2), and models for the redshift distribution (Section 3.3). Readers interested in the astrophysical results may wish to skip ahead to Section 5. In Fig. 1, we provide graphical representations of each mass model described below. Additional details about each model are provided in Appendix B, including their mathematical definitions, lists of hyper-parameters, and their associated prior distributions.', '3.1. Black hole mass distribution': "- · Truncated (4 parameters; Appendix B.1). Our simplest mass model, the distribution of primary masses is a power-law with hard cut-offs at both  \nlow ( m min ) and high ( m max ) masses. The highmass cut-off corresponds to the lower edge of the theorized pair-instability gap in the BH mass spectrum (Woosley & Heger 2015; Heger & Woosley 2002; Heger et al. 2003). The mass ratio distribution is also assumed to be a power law (Kovetz et al. 2017; Fishbach & Holz 2017). In Abbott et al. (2019a), it is referred to as 'Model B.' This model struggles to accommodate high-mass events like GW190521, necessitating more complicated models.  \n- · Power Law + Peak (8 parameters; Appendix B.2). Similar to the Truncated model, but with two modifications. At low masses we implement a smoothing function to avoid a hard cut-off. At high masses, we include a Gaussian peak, initially introduced to model a pile-up from pulsational pair-instability supernovae (Talbot & Thrane 2018). This model is better able to accommodate the high-mass events that pose a challenge for the Truncated model. In Abbott et al. (2019a), it is referred to as 'Model C.'\n- · Broken Power Law (7 parameters; Appendix B.3). The same as Truncated except, instead of a single power-law between m min and m max , the model allows for a break in the powerlaw at some mass m break . This model includes the low-mass smoothing function used in the Power Law + Peak model. The high mass feature m break may correspond to the onset of pair-instability, and the second power-law can be thought of as either a gradual tapering off, or a subpopulation of BHs within the pair-instability mass gap. This model accommodates the high-mass events that pose a challenge for the Truncated while providing an alternative to the Power Law + Peak model.\n- · Multi Peak (11 parameters; Appendix B.4). This phenomenological model is similar to Power Law + Peak except we allow for two peaks. The resulting mass distribution can accommodate hierarchical BBH mergers (Miller & Hamilton 2002b; Gültekin et al. 2004; Antonini & Rasio 2016; Rodriguez et al. 2019), in which second-generation mergers create a second high-mass peak in the BH mass spectrum. We use this model to test if GWTC-2 exhibits evidence for hierarchical mergers.  \n3.2. Spin Distribution  \nFigure 1. Graphical representations of the various mass distributions described in Section 3.1. Multi Spin , a model of both mass and spin, is similar to the mass distribution of Power Law + Peak , with a sharp lower mass cutoff rather than the smooth low mass turn-on.  \n<!-- image -->  \n- · Default (4 parameters; Appendix D.1). This parameterization for the component BH spin magnitudes and tilts was previously used to explore the spin distribution of compact binaries in GWTC1 (Abbott et al. 2019a). The spin of each component BH in a binary is assumed to be independently drawn from the same underlying distribution. The dimensionless spin magnitude is described using a beta distribution as in Wysocki et al. (2019a). The spin tilt distribution from Talbot & Thrane (2017) is a mixture model comprising two components: an isotropic component designed to model dynamically assembled binaries, and a second component in which the spins are preferentially aligned with the orbital angular momentum, as expected for isolated field binaries; the fraction of binaries in the purely-aligned sub-population is denoted ζ . 5 For this latter component, the spin tilt angles are distributed as a truncated Gaussian, with width σ t , that peaks when the BH spin is aligned to the orbital angular momentum. We use this model in concert with the mass models described above.\n- · Gaussian (5 parameters; Appendix D.2). While the Default spin model is physically inspired, this model, based on that of Miller et al. (2020), allows us to fit the distribution of phenomenological spin parameters χ eff (the effective inspiral spin parameter, Eq. 5) and χ p (the effective precession spin parameter, Eq. 6), assuming that their distribution is jointly described as a bivariate Gaussian. The ensemble properties of χ eff and χ p allow us to con-  \nthat the BBHs in GWTC-2 exhibit general relativistic spin-induced precession of the orbital plane ( χ p > 0 ), and that some systems have component spins misaligned by more than 90 · ( χ eff < 0 ) relative to the orbital angular momentum.  \n- · Multi Spin (12 spin parameters, 10 mass parameters; Appendix D.3). This model allows for multiple subpopulations of BBH systems with distinct mass and spin distributions. Specifically, this model assumes a Truncated power-law mass distribution with the additional presence of a 2-D Gaussian subpopulation in m 1 and m 2 , truncated such that m 1 ≥ m 2 . While similar to the Power Law + Peak mass model, there is no smooth turn on and the mass ratio distribution is allowed to differ between each subpopulation. Most importantly, the two subpopulations have independently parameterized Default spin distributions. We use this model to test whether the BBH spin distribution varies as a function of mass as expected if higher-mass systems are the products of hierarchical mergers.", '3.3. Redshift evolution': '- · Non-Evolving (0 parameters). Our default model posits that the merger rate is uniform in comoving volume.\n- · Power-law Evolution (1 parameter; Appendix E). Following Fishbach et al. (2018), the merger rate density is described by a power-law in (1 + z ) where z is redshift. Given the finite range of Advanced LIGO and Advanced Virgo to BBH mergers, we only expect to constrain the redshift evolution at redshifts z glyph[lessorsimilar] 1 (Abbott et al. 2013). The farthest event in our analysis is likely GW190706\\_222641, at redshift z = 0 . 71 +0 . 32 -0 . 27 .  \nTable 1. A summary of the events included in this analysis. Asterisks (*) denote binaries in which both components lie in the NS mass range while a dagger (†) indicates a system in which the nature of the secondary component is unknown. Both of these classes are excluded from our analyses unless explicitly stated. From left to right, the columns show the event name, the 90% credible interval for primary source mass (units of M glyph[circledot] ), the 90% credible interval for secondary mass (units of M glyph[circledot] ), and the minimum available FAR. For events detected in more than one CBC detection pipeline, we report the lowest of the available FAR estimate. These credible intervals are obtained using a prior that is uniform in redshifted component mass and comoving volume, as in Table 6 of Abbott et al. (2020c).  \n×  \nWe adopt a hierarchical Bayesian approach, marginalizing over the properties of individual events to measure parameters of the population models described above; see, e.g., (Thrane & Talbot 2019; Mandel et al. 2019; Vitale 2020). Given data { d i } from N det gravitational-wave detections, the likelihood of the data given population parameters Λ is (Loredo 2004; Mandel et al. 2019; Thrane & Talbot 2019)  \nL ( { d }| Λ , N ) ∝ N N det e -Nξ (Λ) N det ∏ i =1 ∫ L ( d i | θ ) π ( θ | Λ) dθ. (1)  \nHere, N is the total number of events expected during the observation period (including both resolvable and unresolvable signals). Each event is described by a set of parameters θ , the likelihood of which is L ( d | θ ) . The conditional prior π ( θ | Λ) , meanwhile, is defined by the mass, spin, and redshift models described above in Sec. 3. It serves to quantify the predicted distribution of event parameters θ given the hyper-parameters Λ of the population model. An example of a hyper-parameter is the power-law index α governing primary mass distribution π ( m 1 | α ) ∝ m -α 1 for the Truncated mass model. Finally, ξ (Λ) is the fraction of binaries that we expect to successfully detect for a population described by hyperparameters Λ . The procedure for calculating ξ (Λ) is described in Appendix A. One of our primary goals in this paper is to constrain the population hyper-parameters describing the distribution of gravitational-wave signals. Given a log-uniform prior on N , one can marginalize Eq. (1) over N to obtain (Fishbach et al. 2018; Mandel et al. 2019; Vitale 2020)  \nL ( { d }| Λ) ∝ N det ∏ i =1 ∫ L ( d i | θ ) π ( θ | Λ) dθ ξ (Λ) . (2)  \nTo evaluate the single-event likelihood L ( d i | θ ) , we use posterior samples that are obtained using some default prior π ∅ ( θ ) . In this case, integrals over the likelihood can be replaced with weighted averages \' 〈 . . . 〉 \' over discrete samples. Equation (2), for example, becomes  \nL ( { d }| Λ) ∝ N det ∏ i =1 1 ξ (Λ) 〈 π ( θ | Λ) π ∅ ( θ ) 〉 samples , (3)  \nwhere the factor of π ∅ ( θ ) serves to divide out the prior used for initial parameter estimation. We evaluate the likelihoods for population models using the emcee , dynesty , and stan packages (Foreman-Mackey et al. 2013; Speagle 2020; Carpenter et al. 2017; Riddell et al. 2018). The likelihoods are implemented in a variety of  \nsoftware including GWPopulation (Talbot et al. 2019) and PopModels (Wysocki & O\'Shaughnessy 2017). The priors adopted for each of our hyper-parameters are described in Appendix B.  \nThroughout this paper, we find it useful to distinguish between the astrophysical distribution of a parameterthe distribution as it is in nature-and the observed distribution of a parameter-the distribution as it appears among detected events due to selection effects. The posterior population distribution for a given model represents our best guess for the astrophysical distribution of some source parameter θ , averaged over the posterior for population parameters Λ :  \np Λ ( θ ) = ∫ d Λ π ( θ | Λ) p (Λ |{ d } ) . (4)  \nThe subscript Λ indicates that we have marginalised over population parameters. Meanwhile, the posterior predictive distribution refers to the population-averaged distribution of source parameters θ conditioned on detection .  \nIn Appendix C, D and E, we provide posterior predictive checks for our population models. These checks consist of comparing simulated sets of N obs \'predicted" and \'observed" events. For every sample in the model hyper-posterior, we generate a set of predicted events by reweighting our found injections to the population model, and drawing N obs synthetic events. For each hyper-posterior sample, we then generate a set of observed events by reweighting the single-event posteriors of the N obs events to this population prior, and drawing one sample per event. Therefore, the inferred distribution of observed events depends on the population model considered. The first example of such a posterior predictive check is shown in Fig. 19. We calculate the cumulative distribution function for each set of predicted and observed events in the model hyper-posterior, and summarize these curves with 90% credibility bands. The uncertainty in the predicted bands comes from uncertainty in the population hyper-posterior, as well as Poisson fluctuations, because each cumulative distribution curve consists of N obs simulated events.  \nPhenomenological models such as we employ here can fail if their assumed form does not adequately represent reality or if priors are inappropriately chosen. Inferences from such models are inevitably prone to systematic error, given that the model is unlikely to be a perfect description of reality; for instance, the inferred local rate of BBH mergers is subject to a systematic error associated with our choice of model(s) for the BBH mass distribution. While such errors cannot be eliminated, we take steps to control and minimize their magnitude. First, we carry out posterior predictive checks to make sure  \nour models are consistent with the data. Next, where possible, we check for consistency between models with different assumptions, for example, looking for common features in the Power Law + Peak and Broken Power Law models. Finally, we carry out tests to make sure that our inferences are not artifacts of our model design, for example, by applying the models to simulated uninformative data. Additional information about these tests is provided in the Appendix.', '5.1. Mass Distribution': "In this subsection we report results obtained using the mass models described in Section 3.1 (see Fig. 1). 6 We employ the Default spin model and the Non-Evolving redshift model. The results shown here are marginalized over the hyperparameters of the spin distribution.  \nAtruncated power law fails to fit the high-mass BBH events. The Truncated model applied in Abbott et al. (2019a) measured the sharp high-mass cutoff to be m max = 40 . 8 +11 . 8 -4 . 4 M glyph[circledot] . When we fit this model to GWTC-2 data, we obtain m max = 78 . 5 +14 . 1 -9 . 4 M glyph[circledot] , in significant tension ( > 99% credibility) with the GWTC-1 result; see Fig. 2. The Truncated model struggles to accommodate the high-mass systems of GWTC-2. At 99% credibility, three events of GWTC-2 have m 1 > 45 M glyph[circledot] (regardless of whether we use a population-informed prior; Fishbach et al. 2020b). This tension remains (at the > 93% credibility level) even when we exclude the highest-mass event GW190521 (Abbott et al. 2020e,f). The poor fit of the Truncated model is further seen in the posterior predictive check of Fig. 19 in Appendix C, which shows that the Truncated model fails to capture the relative excess of observations with m 1 ∼ 30 M glyph[circledot] compared to the number of events with m 1 glyph[greaterorsimilar] 45 M glyph[circledot] . Our fit to the Truncated model overpredicts the number of observations with m 1 > 45 M glyph[circledot] relative to the number of observations with 30 M glyph[circledot] < m 1 < 45 M glyph[circledot] (98% credibility).  \nThe Power Law + Peak , Broken Power Law , and Multi Peak models provide better fits to the shape of the mass distribution, particularly at high masses. Although our updated fit to the mass distribution extends to higher masses than the GWTC-1 fit, we find that 97 . 1 +1 . 7 -3 . 5 % of BBH systems have primary masses  \nTable 2. Bayes factors for each mass model relative to the favored Power Law + Peak model, which gives the highest Bayesian evidence for GWTC-2. For models that have a smooth turn on at low masses parameterized by δ m , we also compare the corresponding sub-model with a sharp minimum mass cutoff ( δ m = 0 ). For the Power Law + Peak model which includes a fraction λ peak of systems in the Gaussian component, we compare the sub-model with λ peak = 0 . GW190814 is excluded from this analysis.  \n-  \nbelow 45 M glyph[circledot] ( Power Law + Peak model), consistent with the GWTC-1 estimate that 99% of primary masses lie below ∼ 45 M glyph[circledot] (Abbott et al. 2019a; Kimball et al. 2020a). In Table 2, we provide log Bayes factors ( log 10 B ) comparing the mass models; each Bayes factor is measured relative to the model with the highest Bayesian evidence: Power Law + Peak . In each case we use the Default spin model. For context, log 10 B > 1 . 5 is often interpreted as a strong preference for one model over another, and log 10 B > 2 as decisive evidence (Jeffreys 1961).  \nWhile Bayes factors depend on the choice of hyperparameter priors, it is nonetheless possible to see that the Truncated model is disfavored compared to the more complicated models. This inference is driven in part by the fact that-in our posterior predictive checks94% of the time, the Truncated model overpredicts the number of detections with m 1 > 50 M glyph[circledot] . See Appendix C for information about the posterior predictive checks. Meanwhile, there is not a strong preference for Power Law + Peak over Broken Power Law or Multi Peak . We currently lack the resolving power to determine whether the deviations from Truncated are best modeled as a break, a Gaussian peak, or two Gaussian peaks. As a further check, we carried out a follow-up analysis using a hybrid Broken Power Law + peak model, which indicated only modest support for a peak on top of the Broken Power Law distribution ( log 10 B = 0 . 79 ).  \nThere are features in the black hole mass spectrum beyond a power-law. Figure 3 shows the astrophysical merger rate density as a function of primary BH mass for the Truncated , Power Law + Peak ,  \n/circledot  \nTable 3. The m 1% and m 99% credible intervals (90%) for various mass models and combinations of events. These variables are defined such that, among the astrophysical BBH population, 1% of systems have primary masses m 1 ≤ m 1% , while 99% have primary masses m 1 ≤ m 99% . The Power Law + Peak , Multi Peak , and Broken Power Law models are preferred over the Truncated model.  \nFigure 2. Posterior for the maximum mass using GWTC-1 and fit to the Truncated model (blue), compared to the posterior obtained by adding events from O3a data (red). The two distributions are inconsistent, suggesting the Truncated model is inadequate. The tension between GWTC-1 and GWTC-2 is somewhat alleviated by the exclusion of the high-mass event GW190521 (green). However, there remain several other high-mass events in O3a. The black dashed lines show primary mass posteriors for the three events in which m 1 > 45 M glyph[circledot] at 99% credibility (we employ a prior that is uniform in redshifted masses). These events cause a significant shift in the m max posterior if we assume a simple power-law fits the data.  \n<!-- image -->  \nMulti Peak and Broken Power Law models. Figure 4, meanwhile, shows the posterior predictive distribution for primary masses, including selection effects. Corner plots showing the constraints for the parameters in each model are available in Appendix B; see Figs. 16, 17 and 18. In Appendix C, we show posterior predictive  \nchecks for each model (Fig. 19 and Fig. 23), comparing mock observations predicted by the model to the empirical distribution inferred from GWTC-2.  \nWe turn first to the Broken Power Law model (second panel in Fig. 3), which is characterized by two spectral indices, α 1 and α 2 , with p ( m 1 ) ∝ m -α 1 1 for m 1 < m break , and p ( m 1 ) ∝ m -α 2 1 above the break. We find the data prefer a break at m break = 39 . 7 +20 . 3 -9 . 1 M glyph[circledot] ; 90% credible bounds on the location of this break are denoted by the gray vertical band in the second panel of Fig. 3. For masses above the break, the Broken Power Law model prefers a significantly steeper powerlaw slope, from α 1 = 1 . 58 +0 . 82 -0 . 86 before to α 2 = 5 . 6 +4 . 1 -2 . 5 after. Figure 5 shows the joint posterior on α 1 and α 2 . We infer that α 2 > α 1 at credibility 98% . The break aligns with the cutoff m max inferred with GWTC1 data (Abbott et al. 2019a), and we speculate that the steep drop-off in the merger rate that occurs after m break may be an imprint of PPSN, which are expected to become important for BH masses around 35 M glyph[circledot] (Woosley & Heger 2015; Heger & Woosley 2002; Heger et al. 2003).  \nThe Power Law + Peak model (third panel of Fig. 3) produces a qualitatively similar fit to the one obtained from the Broken Power Law model. However, a key feature of the Power Law + Peak model is the Gaussian peak at 33 . 1 +4 . 0 -5 . 6 M glyph[circledot] , denoted by the gray vertical band in the third panel of Fig. 3. Evidence for a peak can be seen in Fig. 6, which shows the posterior for λ peak , the fraction of systems that belong to the Gaussian component. We see that λ peak = 0 (pure power law) is disfavored. It was envisioned (Talbot & Thrane 2018) that the power-law component of the Power Law + Peak model would terminate in the vicinity of this peak to create a high-mass gap. However, in order to  \n100  \nFigure 3. Astrophysical primary BH mass distribution for the Truncated , Broken Power Law , Power Law + Peak and Multi Peak models. The solid curve is the posterior population distribution (averaging over model uncertainty) while the shaded region shows the 90% credible interval. While the median rate is always inside the credible region, the solid curve represents the mean , which can be outside the credible region. Top (navy) is the Truncated model, second from the top (green) is the Broken Power Law model, third from the top (blue) is for the Power Law + Peak model, and bottom (red) is for the Multi Peak model. The Truncated model is disfavored compared to the three latter models that predict a feature at ∼ 40 M glyph[circledot] : a break in the mass spectrum in the Broken Power Law model or additional Gaussian peaks in the Power Law + Peak and Multi Peak models. The vertical gray bands show 90% credible bounds on the locations of these additional features.  \n<!-- image -->  \n/circledot  \nα  \n<!-- image -->  \n/circledot  \nFigure 4. Observed primary BH mass distributions predicted by each mass model. For each model, we average over the uncertainty in the hyper-parameter posterior. The observed distribution describes the events successfully detected by LIGO-Virgo, preferentially favoring more massive systems relative to the astrophysical distribution due to selection effects.Figure 5. Constraints on the power-law indices governing the primary mass distribution within the Broken Power Law model. The parameter α 1 is the power-law index below the break, which is found to be m break = 39 . 7 +20 . 3 -9 . 1 while α 2 is the index above the break. The dashed and solid contours mark the central 50% and 90% posterior credible regions, respectively, under a flat prior on α 1 , α 2 in the range ( -4 , 12) . We rule out with high confidence the hypothesis that α 1 = α 2 , indicated by the dashed diagonal line, finding α 2 > α 1 with 98 % credibility.  \n<!-- image -->  \naccommodate the most massive binaries in GWTC-2, the power-law extends to values of m max = 86 +12 -13 M glyph[circledot] .  \nWhile the mass spectrum must extend to these high masses, we find that 99% of primary BH masses lie below m 99% ∼ 60 M glyph[circledot] ; see Table 3. The astrophysical rate density at ∼ 80 M glyph[circledot] (the primary mass of GW190521) is two orders of magnitude lower than the rate density at ∼ 40 M glyph[circledot] . However, because of selection effects, the posterior predictive distribution skews to much higher  \nmasses, as seen in Fig. 4, so that the probability of detecting at least one event with m 1 ≥ 80 M glyph[circledot] after observing 44 BBH events drawn from the Power Law + Peak posterior predictive distribution of Fig. 4 is high: 32% .  \nFigure 6. Posterior distribution on the fraction of binaries ( λ peak ) in the Gaussian component of Power Law + Peak model, under a flat prior on λ peak ; see Appendix B.2. We find that λ peak = 0 (which corresponds to no Gaussian peak) is disfavored, supporting the hypothesis that there is a feature in the BH primary mass spectrum.  \n<!-- image -->  \nWe cannot determine whether the high-mass events of GWTC-2 belong to a distinct subpopulation rather than a high-mass tail of the normal BBH population. An additional subpopulation may be expected if highmass BHs have a different origin from low mass ones; for example, if they are the products of hierarchical mergers (Fishbach et al. 2017; Gerosa & Berti 2017;  \nChatziioannou et al. 2019; Doctor et al. 2019; Kimball et al. 2020a). Using the Multi Peak model (bottom panel in Fig. 3), which allows for a second high-mass Gaussian component at m 1 > 50 M glyph[circledot] in addition to the Gaussian component at m 1 glyph[lessorsimilar] 40 M glyph[circledot] , we find that the addition of a second Gaussian peak is not preferred by the data. The Multi Peak model is mildly disfavored compared to Power Law + Peak , with a log 10 B = -0 . 3 (or Bayes factor B = 0 . 5 ) for Multi Peak relative to Power Law + Peak . Motivated by the hypothesis of hierarchical mergers, we consider a variation of the Multi Peak model, in which the location of the second peak is required to be at twice the value of the first peak; that is, µ m, 2 = 2 µ m, 1 (see Appendix B.4. Studying the ( µ m, 1 , µ m, 2 ) panel of the corner plot in Fig. 18, we see that the data mildly prefer the second peak at ≈ 70 M glyph[circledot] , consistent with twice the value of the first peak at ≈ 35 M glyph[circledot] . This ' Modified MultiPeak ' is mildly preferred over the original version by a Bayes factor of ∼ 2 . A similar conclusion is found using the Multi Spin model; as discussed in Section 5.2, we find hints, but no significant evidence for subpopulations with distinct spin distributions. Additional evidence for hierarchical mergers is presented in Kimball et al. (2020b).  \nWithin the framework of the Broken Power Law , Power Law + Peak and Multi Peak models, the most massive event, GW190521, appears to be a normal member of the BBH population in the context of the other GWTC-2 events (see Appendix C.2). The event GW190521 is an outlier if we consider it in the context of GWTC-1 with the Truncated model, but we interpret this as a limitation of the Truncated model (see Fig. 2; Abbott et al. 2020f).  \nThe GWTC-1 detections showed that BHs more massive than ∼ 45 M glyph[circledot] merge relatively rarely, based on simple extrapolations from below 45 M glyph[circledot] . With GWTC2, we are beginning to resolve the shape of the primary BH mass spectrum above 45 M glyph[circledot] . The implications are not yet clear, but there are intriguing possibilities. One hypothesis is that the events with m 1 > 45 M glyph[circledot] are simply the high-mass tail of the ordinary BBH population, and do not form through a distinct channel. For example, if the lower edge of the PPSN gap may be modeled as a smooth tapering rather than a sharp cutoff, the feature at m break = 39 . 7 +20 . 3 -9 . 1 may represent the onset of pair-instability. This explanation may pose challenges to our understanding of stellar evolution since the pairinstability cutoff of BH masses at ∼ 40 M glyph[circledot] is thought to be relatively abrupt (Woosley et al. 2002; Woosley 2017; Farmer et al. 2019), even though its precise location is uncertain (Mapelli et al. 2020; Spera & Mapelli 2017; Giacobbo et al. 2017; van Son et al. 2020; Farmer et al.  \n2020; Croon et al. 2020; Marchant & Moriya 2020). If the PPSN cutoff is indeed sharp and all observed BBH systems lie below the PPSN gap, the cutoff must occur at relatively high masses; in the Truncated model, m max = 78 . 5 +14 . 1 -9 . 4 M glyph[circledot] (or, excluding the most massive event, GW190521, m max = 57 . 0 +11 . 9 -6 . 6 M glyph[circledot] ). This may have significant implications for nuclear (Farmer et al. 2020) and particle (Croon et al. 2020; Ziegler & Freese 2020) physics.  \nAnother hypothesis is that the events with m 1 > 45 M glyph[circledot] constitute a distinct population, created, for example, from hierarchical mergers of lower mass binaries in globular clusters or galactic nuclei (Miller & Hamilton 2002a; Antonini & Rasio 2016; Kimpson et al. 2016; Rodriguez et al. 2018; Yang et al. 2019; Arca Sedda 2020a). Alternatively, the high-mass gap might be populated from low-metallicity stellar mergers in young star clusters, the remnants of which can merge dynamically (Carlo et al. 2019a,b), BH growth through accretion (Roupas & Kazanas 2019; Rice & Zhang 2020; Natarajan 2020; Safarzadeh & Haiman 2020), or Population III stellar remnants with large hydrogen envelopes (Tanikawa et al. 2020; Farrell et al. 2020; Liu & Bromm 2020).  \nThe BBH primary mass distribution exhibits a global maximum between ∼ 4 and ∼ 10 M glyph[circledot] . Figure 7 shows the joint posterior for the m min and δ m parameters inferred using the Power Law + Peak and Broken Power Law mass models, including only BBH events with m 2 > 3 M glyph[circledot] . Recall that, while the Truncated model has a sharp cutoff at m min , the remaining models implement a smooth turn-on of width δ m above m min , causing the mass spectrum to peak and turn over between m min , and δ m + m min (Talbot & Thrane 2018). Using both the Broken Power Law model and the Power Law + Peak model, we find that the primary mass spectrum does not decrease monotonically from 3 M glyph[circledot] . Rather, it turns over at 7 . 8 +1 . 8 -2 . 0 M glyph[circledot] ( Power Law + Peak model) or 6 . 02 +0 . 78 -1 . 96 M glyph[circledot] ( Broken Power Law model). In other words, the mass distribution must turn over at m 1 > 3 M glyph[circledot] , with 99 . 9% credibility (assuming the Power Law + Peak model) or 98 . 5% credibility ( Broken Power Law model). As seen in Fig. 7, if the BH low-mass cutoff is sharp ( δ m = 0 ), then m min glyph[greaterorsimilar] 4 M glyph[circledot] . Conversely, if the BH mass spectrum extends below m min glyph[lessorsimilar] 4 M glyph[circledot] , an extended turn-on δ m glyph[greaterorsimilar] 3 M glyph[circledot] is required. These results support the existence of a low-mass gap (or dip) between ∼ 2 . 6 M glyph[circledot] (the secondary mass of GW190814; the most massive component mass observed below 3 M glyph[circledot] ) and ∼ 6 M glyph[circledot] , strengthening results from Fishbach et al. (2020a), although we cannot determine whether the low-mass gap is empty.  \nFigure 7. Posterior distribution for population parameters m min , the minimum BH mass, and δ m , which controls the sharpness of the low-mass cut-off. A sharp cut-off corresponds to δ m = 0 . Analyzing the 44 BBH events (with the exclusion of GW190814), both models exclude ( m min = 3M glyph[circledot] , δ m = 0) with > 99% credibility, indicating that the rate drops off at low masses. To varying degrees, both models allow for m min ≤ 3 M glyph[circledot] , δ m > 0 , suggesting that the low-mass gap may not be empty. See Appendices B.3 and B.2 for additional details about the δ m , m min parameters.  \n<!-- image -->  \n/circledot  \n/circledot  \nSince our models do not permit additional features beyond a smooth turn on to a power law at low masses, they struggle to accommodate GW190814 (with secondary mass at m 2 = 2 . 59 +0 . 08 -0 . 09 ). We can see this by comparing the mass distribution inferred from the events with m 1 ≥ m 2 > 3 M glyph[circledot] , discussed above, to the distribution inferred with GW190814. If GW190814 is a BBH system, the minimum BH mass must extend to m min = 2 . 18 +0 . 27 -0 . 16 M glyph[circledot] (see the dashed histograms in Fig. 8a). In Fig. 8b we show how the inclusion/exclusion of GW190814 affects the shape of the primary mass distribution below glyph[lessorsimilar] 5 M glyph[circledot] . We see that the two distributions are inconsistent at the low mass end, suggesting that there is a feature in the mass distribution between ∼ 2 . 6 M glyph[circledot] and ∼ 6 M glyph[circledot] that our models cannot capture. This effect can also be seen in the m 1% values inferred with/without GW190814, shown in Table 3. Assuming the Truncated and Power Law + Peak models, the m min posteriors inferred with GW190814 lie at the 0 . 25 +0 . 44 -0 . 20 and 0 . 78 +2 . 22 -0 . 69 percentiles, respectively, of the m min posteriors obtained without GW190814. Even using the Broken Power Law model, which admits greater overlap in the 1-dimensional m min posteriors inferred with/without GW190814, the addition of GW190814 significantly shifts the two-dimensional ( δ m , m min ) pos-  \nand the inferred mass spectrum. Assuming the Broken Power Law model with GW190814, the mass at which the mass spectrum turns over is shifted down to 2 . 59 +0 . 78 -0 . 39 M glyph[circledot] , which is inconsistent with the turnover mass (the low-mass local maximum) inferred without GW190814, 6 . 02 +0 . 78 -1 . 96 M glyph[circledot] . This indicates a failure of our models to fit GW190814 together with the BBH systems of GWTC-2. This finding is supported by additional studies described in Appendix C.3. Because GW190814 is a population outlier with respect to the BBH events of GWTC-2 and our choice of models, we exclude it from the analyses here unless otherwise indicated.  \nThe distribution of mass ratios is broad. The GWTC-1 events are all individually consistent with q = 1 . Describing the conditional mass-ratio distribution as a power law p ( q | m 1 ) ∝ q β q , a population analysis of GWTC-1 allowed β q = 12 (our maximum prior bound), consistent with a mass-ratio distribution sharply peaked at equal-mass pairings (Abbott et al. 2019a; Roulet & Zaldarriaga 2019; Fishbach & Holz 2020b). GWTC-2 saw the first detections of confidently asymmetric systems: GW190412 (Abbott et al. 2020a) and GW190814 (Abbott et al. 2020b). Excluding GW190814, our reconstruction of the mass ratio distribution is consistent with the results published in Abbott et al. (2020a): β q = 1 . 3 +2 . 4 -1 . 5 for the Power Law + Peak model and β q = 1 . 4 +2 . 5 -1 . 5 for the Broken Power Law model. We rule out distributions that are sharply peaked around q = 1 , with β q < 2 . 9 ( Power Law + Peak ) and β q < 3 . 1 ( Broken Power Law ) at 90% credibility. However, we also disfavor distributions that prefer unequal-mass pairings, with β q > 0 at 92 % ( Power Law + Peak ) and 94 % ( Broken Power Law ) credibility. We find that 90% of systems in the underlying population have mass ratios q > 0 . 26 +0 . 14 -0 . 08 .", '5.2. Spin Distribution': "In this subsection, we highlight the results from the Gaussian , Default , and Multi Spin models. We fix the redshift distribution to a Non-Evolving merger rate. The Gaussian and Multi Spin models assume the mass distributions described in Appendix D.2 and D.3, respectively. For the Default spin model, we employ the Power Law + Peak mass model, simultaneously fitting the mass and spin distribution as in the previous subsection.  \nWe observe spin-induced general relativistic precession of the orbital plane. As two BHs merge, the morphology of the resulting gravitational waveform depends on their spins. The spin-dependence of a gravitational waveform is determined in part by two phenomenological parameters. First, the effective inspiral spin pa-  \nFigure 8. (a) Posterior distribution for the minimum mass parameter in the Truncated (navy), Broken Power Law (green) and Power Law + Peak (blue) models. The solid bars show the posterior for m min when fitting the models to the confident BBH events excluding GW190814, while the dashed lines show the fits to the BBH events including GW190814. The Truncated model is disfavored. (b) Distribution of primary masses inferred using the Power Law + Peak model when including (navy, dashed line) and excluding (light blue, solid line) GW190814. Only the secondary mass of GW190814 is below 3 M glyph[circledot] . However, the primary and secondary mass distributions share a common m min parameter in all models we consider. The effect of GW190814 on the mass spectrum at low masses inferred using the Broken Power Law model is similar.  \n<!-- image -->  \n/circledot  \n(a) Minimum primary mass distributions.  \n<!-- image -->  \n/circledot  \n(b) Power Law + Peak mass distribution fit.  \nrameter χ eff quantifies the spin components aligned with the orbital angular momentum (Damour 2001):  \nχ eff = χ 1 cos θ 1 + q χ 2 cos θ 2 1 + q . (5)  \nHere, χ 1 and χ 2 are the dimensionless component spins, defined by χ i = | cS i / ( Gm 2 i ) | where S i is the spin angular momentum of component i , and θ 1 and θ 2 are the misalignment angles between the component spins and the orbital angular momentum. Second, spins with components perpendicular to the orbital angular momentum drive relativistic precession of the orbital plane (Apostolatos et al. 1994). The effect is quantified by the effective precession spin parameter (Schmidt et al. 2012; Hannam et al. 2014; Schmidt et al. 2015)  \nχ p = max [ χ 1 sin θ 1 , ( 4 q +3 4 + 3 q ) q χ 2 sin θ 2 ] . (6)  \nA non-zero value of χ p indicates the presence of relativistic spin-induced precession of the orbital plane. Although the component spin tilts θ 1 and θ 2 appearing in Eqs. (5) and (6) generically evolve over the course of a binary inspiral, χ eff and χ p are themselves approximately conserved quantities (Kidder 1995; Schmidt et al. 2015).  \nThe first unambiguous measurements of BH spin in gravitational-wave astronomy came from analyses of the BBH event GW151226. This system had χ eff > 0 at 99% credibility, with at least one of its components having spin magnitude χ > 0 . 2 , and spin misalignment  \nangles consistent with θ 1 = θ 2 = 0 (Abbott et al. 2016b). While analyses of GWTC-1 found no clear evidence for spin in the other events in GWTC-1 (Miller et al. 2020, but also see Zackay et al. 2019 and Huang et al. 2020), GWTC-1 is collectively inconsistent with a population of non-spinning BHs, if one allows for both spinning and non-spinning subpopulations (Kimball et al. 2020a). Moreover, population analyses of GWTC-1 mildly disfavor the scenario in which all spins are perfectly aligned ( θ 1 = θ 2 = 0 ), although the degree of misalignment is degenerate with the spin magnitude distribution (Farr et al. 2017; Farr et al. 2018; Tiwari et al. 2018; Abbott et al. 2019a; Wysocki et al. 2019b).  \nIn GWTC-2, additional BBH events are observed with confidently positive effective inspiral spin parameter. No individual event is observed with confidently negative χ eff (Abbott et al. 2020c). Several events, including GW190521 (Abbott et al. 2020e,f) and GW190412 (Abbott et al. 2020a), show moderate evidence for non-zero χ p , but no single event unambiguously exhibits spininduced precession (Abbott et al. 2020c).  \nUsing the Default model described in Sec. 3 (see also Appendix D.1), we obtain evidence for non-vanishing spin-orbit misalignment among the population of BBH events in GWTC-2. The Default model describes the distribution of spin-orbit misalignments as a mixture between two components: a component with isotropicallyoriented spins and a preferentially-aligned component  \n<!-- image -->  \nFigure 9. Left : Joint posterior on the fraction ζ of BBHs with preferentially-aligned spins (versus isotropic spins) and the spread σ t of misalignment angles among this population obtained using the Default spin model (see Appendix D.1 for additional details). We rule out a population with perfectly aligned spins corresponding to ζ = 1 and σ t = 0 . The gray shaded region represents the region of parameter space inaccessible to our analysis. This region is artificially excluded due to sampling uncertainties even when analyzing uninformative samples drawn from the spin tilt prior. The dashed and solid contours mark the central 50% and 90% posterior credible regions, respectively, assuming a flat prior on ζ and σ t . Right : Population predictive distributions for the effective precession spin parameter χ p of BBH systems obtained using the Gaussian (blue) and Default (orange) spin models. Shaded regions show the central 90% credible bounds on p ( χ p ) at a given spin value, while the solid lines show the median posterior prediction. The inset shows draws of the Gaussian χ p distributions implied by the posterior on µ p and σ p . Broadly, we see support for two possible morphologies, indicated schematically by the dashed black curves. GWTC-2 is compatible with a χ p distribution that is either broad, or one that is narrow and centered at µ p ∼ 0 . 2 .  \n<!-- image -->  \nwith z = cos θ values centered at z = 0 (perfect alignment) with a Gaussian spread of width σ t . In the left side of Fig. 9, we show the joint posterior on σ t and the fraction ζ of events in the preferentially-aligned subpopulation. Perfect alignment corresponds to ζ = 1 and σ t = 0 . We see that this case is ruled out at > 99% credibility. Thus, either a non-zero fraction of BBH events exhibit isotropically-oriented spins, or BBH spins are preferentially aligned to their orbits but with a nonvanishing spread. Either case constitutes an observation of in-plane spin components among the BBH population. The shaded gray tiles in the left side of Fig. 9 show the values of ζ and σ t that are artificially excluded by the prior and/or finite sampling effects; the true measurement using GWTC-2 lies well away from this artificial exclusion region. 7  \nIn Fig. 10, discussed further below, we plot the range of component spin magnitude and tilt angle distributions recovered using the Default model. Although the data are consistent with tilt angle distributions that  \nfavor alignment, distributions that are highly peaked at cos θ 1 , 2 = 1 are ruled out.  \nA similar conclusion regarding the presence of in-plane spin components may be drawn using Gaussian spin model, which imposes an entirely different parameterization for the BH spin distribution and makes different assumptions regarding their masses. In particular, when measuring the mean µ p and standard deviation σ p of the χ p distribution, the case µ p = σ p = 0 is ruled out at > 99% credibility; fewer than 1% of posterior samples occur at µ p ≤ 0 . 05 and σ p ≤ 0 . 05 . Since any non-zero µ p or σ p implies the existence of spin-induced precession, this result supports the observation of spin misalignment seen in the Default model. In the right side of Fig. 9, the dark blue curve and shaded blue region mark the median and 90% credible bound, respectively, on p ( χ p ) as inferred by the Gaussian model. While the blue region in this figure suggests a χ p distribution that peaks at ∼ 0 . 2 , there are in fact two morphologies preferred by the data according to the Gaussian model: the recovered χ p distribution is either broad-or narrowly peaked at χ p ≈ 0 . 2 . This is illustrated by the inset, in which we plot an ensemble of distributions corresponding to individual draws from the ( µ p , σ p ) posterior; the dashed black curves highlight traces representative of the two permitted morphologies.  \nFor comparison, the orange curve in the right side of Fig. 9 shows the χ p distribution implied by the Default results discussed above. There are several potentially meaningful differences between the results from the Gaussian and Default models. In particular, the Default model predicts χ p distributions that are generally broader and peaked at lower values. This is due to additional physical constraints imposed by the Default spin model; component spins are presumed to preferentially cluster about θ = 0 , an assumption that preferentially favors smaller χ p values. Nevertheless, the two models agree well within statistical uncertainties, 8 indicating that the identification of spin-induced precession is robust to the systematic modeling choices and prior uncertainties.  \nAs mentioned above, GW190521 and GW190412 individually show mild evidence of precession, with χ p posteriors shifted away from their respective priors (Abbott et al. 2020a,e,f). To verify that our populationlevel conclusions are not driven primarily by these two events, we have repeated the Gaussian analysis excluding GW190521 and GW190412. Our results again exclude µ p = σ p = 0 at a similar level of confidence ( > 99% credibility). This implies that the signature of precession observed here is due to the combined influence of many systems with only weakly measured χ p , consistent with expectations from simulation studies (Fairhurst et al. 2019; Wysocki et al. 2019b).  \nThe injection sets used to quantify search selection effects (see Appendix A) contain only events whose component spins are perfectly aligned with their orbital angular momenta. The results in Fig. 9 therefore do not account for systematics possibly affecting our ability to detect events with misaligned spins. The matched filter template banks adopted by the GstLAL and PyCBC search pipelines, for instance, are composed of purely alignedspin waveforms, and so may have reduced sensitivity to events with high χ p (Harry et al. 2016; Calderón Bustillo et al. 2017). Selection effects can, however, only decrease the efficiency with which events with large in-plane spins are detected; incorporating such effects would further shift the posterior in Fig. 9 away from ζ = 1 and σ t = 0 and/or more strongly rule out a delta function at χ p = 0 . Thus, the presence of in-plane spin components is robust to selection effects. The specific preference for µ p ≈ 0 . 2 , though, may not be. In the future, accurately characterizing the effects of in-plane spins on detection efficiency  \nwill be crucial in order to robustly determine the shape of the cos θ and χ p distributions.  \nWe observe anti-aligned spin, which may suggest the presence of more than one binary formation channel. Using the Gaussian model described in Sec. 3, we infer the presence of systems with negative effective inspiral spin parameter: χ eff < 0 . Thus, there exist BBH systems with at least one component spin tilted by θ > 90 · relative to the orbital angular momenta. Figure 11 shows posteriors for the mean µ eff and standard deviation σ eff of the χ eff distribution, marginalized over µ p , σ p , and the covariance between the effective inspiral spin parameter and the effective precession spin parameter. With a peak at µ eff = 0 . 06 +0 . 05 -0 . 05 , we find that most systems have small but positive χ eff , in agreement with the inference from GWTC-1 (Miller et al. 2020; Roulet & Zaldarriaga 2019). With GWTC-2, we can now also constrain the width of the χ eff distribution. The result, σ eff = 0 . 12 +0 . 06 -0 . 04 , requires that a nonzero fraction of BBH systems have χ eff < 0 . Unlike the constraints on the χ p distribution presented above, the results for the presence of negative effective inspiral spin parameter do incorporate selection effects via the prescription described in Sec. 4.  \nAnalysis with the Default spin model is also suggestive of an anisotropic distribution of spin orientations. In Fig. 10, we plot the population distribution of cos θ 1 , 2 reconstructed using the Default model. While the cos θ 1 , 2 distribution shows a preference for primarily aligned spins, with cos θ 1 , 2 > 0 , it also exhibits non-vanishing posterior support for cos θ 1 , 2 < 0 , indicating the presence of component spins misaligned by more than 90 · . The χ eff distribution inferred with the Default model closely matches the distribution inferred using the Gaussian model; compare the orange and blue bands in the right panel of Fig. 11. The two models therefore agree on the fraction of systems with anti-aligned component spins.  \nTo further verify that the apparent presence of events with negative χ eff is physical and not an artifact of our choice of models, we repeat our inference of the Gaussian χ eff distribution, this time permitting the minimum allowed effective inspiral spin parameter χ min eff (until now fixed to χ min eff = -1 ) to vary as an additional hyper-parameter to be inferred from the data. When fitting for χ min eff alongside µ eff and σ eff , we find that χ min eff is less than zero at 99% credibility (see Fig. 27 in the Appendix), confirming that the evidence for anti-aligned spin is not an artifact of our parameterization. Allowing χ min eff to vary yields similar results for the implied χ eff distribution, and in particular, the fraction of systems with negative χ eff .  \n<!-- image -->  \nFigure 10. Reconstructions of the BH spin magnitude and tilt distributions. Left : The distribution of dimensionless spin magnitude χ as inferred using the Default spin model (see Appendix D.1). Light traces show individual draws from the Default posterior, while the solid black curve shows the posterior population distribution for χ . Dashed lines mark the central 90% quantiles. Right : the reconstructed distribution of tilt angle cos θ 1 , 2 of BH component spins relative to the orbital angular momenta. An isotropic spin orientation, which corresponds to a uniform distribution in cos θ 1 , 2 , is disfavored but not ruled out. The data do, however, rule out a highly peaked distribution at cos θ 1 , 2 = 1 . Rather, the data are consistent with a gently peaked distribution, with a modest preference for aligned spin ( cos θ 1 , 2 > 0 ).  \n<!-- image -->  \n<!-- image -->  \nFigure 11. Left : Posterior for the mean µ eff and standard deviation σ eff of the BBH χ eff distribution, obtained using the Gaussian model described in Appendix D.2. We marginalize over the parameters governing the distribution of the effective inspiral spin parameter and the effective precession spin parameter. While we infer a χ eff distribution that is peaked at positive values, its measured width implies that a non-zero fraction of BBH systems have negative χ eff , implying component spins misaligned by t 1 , 2 > 90 · relative to the orbital angular momentum. Right : Population predictive distributions for the effective inspiral spin parameter χ eff obtained with both the Gaussian and Default spin models. Shaded regions show the central 90% credible bounds on p ( χ eff ) and the solid lines show the median posterior prediction for the χ eff distribution.  \n<!-- image -->  \nThe presence of BBH systems with negative effective inspiral spin parameter carries implications for the formation channels that give rise to stellar-mass BBH mergers. BBHs born in the field from isolated stellar progenitors are predicted to contain components whose spins are nearly aligned with their orbital angular momenta, although sufficiently strong supernova kicks might produce  \nmodest misalignment (O'Shaughnessy et al. 2017; Stevenson et al. 2017; Gerosa et al. 2018; Rodriguez et al. 2016; Bavera et al. 2019). In contrast, binaries assembled dynamically in dense stellar environments are expected to have randomly oriented component spins, yielding positive or negative χ eff with equal probabilities (Kalogera 2000; Mandel & O'Shaughnessy 2010; Rodriguez et al.  \n2016; Zevin et al. 2017; Rodriguez et al. 2018; Doctor et al. 2019).  \nUsing the posteriors for µ eff and σ eff from Fig. 11, in Fig. 12 we show posteriors for the implied fractions of BBH systems with negative ( χ eff < -0 . 01 ) and positive effective inspiral spin parameter ( χ eff > 0 . 01 ). Motivated by recent work suggesting that BHs are born with natal spins as small as χ ≈ 10 -2 (Qin et al. 2018; Fuller et al. 2019; Fuller & Ma 2019; Bavera et al. 2019), as well as the tendency of vanishingly small spins to confound efforts to distinguish between positive and negative χ eff (Farr et al. 2018; Abbott et al. 2019a), we include a third bin containing vanishingly small spins between -0 . 01 ≤ χ eff ≤ 0 . 01 . At 90% credibility, we find that fractions f p = 0 . 67 +0 . 16 -0 . 16 , f n = 0 . 27 +0 . 17 -0 . 15 , and f v = 0 . 05 +0 . 02 -0 . 01 of BBH systems have positive, negative, and vanishing χ eff , respectively. All three posterior distributions are peaked away from zero. In particular, f n > 7% at 99% credibility. This result is in contrast to results obtained using GWTC-1 alone, which did not exhibit a confidently non-zero fraction of events with negative χ eff (Abbott et al. 2019a; Miller et al. 2020). Additionally, the relatively small fraction f v of binaries with vanishing spins may provide clues about how BHs gain angular momentum, given recent studies suggesting that most BHs are born slowly rotating (Fuller & Ma 2019). While we define the vanishing bin to be -0 . 01 ≤ χ eff ≤ 0 . 01 , the exact choice of width for the vanishing bin does not strongly affect the values of f p and f n , relative to one another. We obtain nearly identical results, albeit with a slightly weaker lower bound on f n , when we additionally allow the minimum effective inspiral spin χ min eff to vary as described above; in this case f p = 0 . 65 +0 . 17 -0 . 16 and f n = 0 . 26 +0 . 17 -0 . 21 .  \nAs mentioned above, dynamical formation in dense clusters is not the only astrophysical explanation of negative effective inspiral spin parameter. If stellar progenitors experience both strong natal kicks in supernovae and inefficient spin realignment, glyph[lessorsimilar] 10% of BBH systems formed through isolated binary evolution may have χ eff < 0 (Rodriguez et al. 2016; O'Shaughnessy et al. 2017; Stevenson et al. 2017; Wysocki et al. 2018), although these results depend on the poorly understood physics of natal kicks and binary interaction via torques and mass transfer. Moreover, we have so far neglected the possibility of other formation channels that may operate in both the field and dynamical regimes. Isolated hierarchical triples, for example, may produce binary mergers with preferentially in-plane component spins (Rodriguez & Antonini 2018; Antonini et al. 2018; Liu et al. 2019; Fragione & Kocsis 2020). Mergers in the disks of active galactic nuclei, meanwhile, yield component spins that are preferentially parallel or anti-parallel to a binary's  \norbital angular momentum (McKernan et al. 2018; Yang et al. 2019; McKernan et al. 2020).  \nWith these qualifications in mind, if we interpret negative χ eff as indicative of dynamical formation in stellar clusters, then our constraints on f n can be used to infer the fraction of dynamically assembled binaries. We assume that dynamical assembly in dense stellar environments yields a χ eff distribution that is symmetric about zero, while isolated binary evolution produces only positive χ eff . Among the binaries with non-negligible spin (excluding those in the 'vanishing' category above), the fractions f d and f i of binaries arising from dynamical and isolated channels are  \nf d = 2 f n f p + f n , (7)  \nf i = f p -f n f p + f n . (8)  \nWe find 0 . 25 ≤ f d ≤ 0 . 93 at 90% credibility, suggesting that both the field and the dynamical cluster scenarios contribute to the BBH mergers observed in GWTC-2. Because the relative values of f n and f p are not sensitive to the width of the vanishing χ eff bin, this conclusion does not depend strongly on the definition of vanishing spin.  \nAt present, we are unable to include a systematic investigation of waveform error in our analysis of anti-aligned spin and orbital precession. However, preliminary studies suggest that waveform error is unlikely to significantly affect this and other results in this paper. As described in Abbott et al. (2020c), there is good agreement regarding the parameters of the GWTC-2 events when inferred with different waveforms. There is a caveat: our studies do not include eccentric waveforms. Romero-Shaw et al. (2020a); Gayathri et al. (2020) suggest that GW190521 may have been eccentric and Calderón Bustillo et al. (2020) point out that eccentricity can be confused with precession for high-mass events. Currently the only event likely to be affected is GW190521 (Calderón Bustillo et al. 2020). It is not clear how our results would change if we accounted for eccentricity, but we note that eccentricity can be a signature of dynamical assembly (Samsing et al. 2014; Samsing & Ramirez-Ruiz 2017; Samsing 2018; Romero-Shaw et al. 2019; Lower et al. 2018; Fragione & Kocsis 2019; Zevin et al. 2019). We also note that Fishbach & Holz (2020a); Nitz & Capano (2020) find that GW190521 may be an intermediate mass ratio inspiral, which could potentially alter our conclusions if true. For future work, it would be worthwhile to estimate the systematic error using different waveform approximants.  \nNo strong evidence for variation of the spin distribution with mass. BHs born in hierarchical mergers  \nFigure 12. Posterior distribution for the fraction of BBH events with positive or negative χ eff (corresponding to alignment or anti-alignment of the BH spin with the orbital angular momentum; see Eq. 5). We also include the fraction of events that is consistent with vanishingly small spins | χ eff | < 0 . 01 . We confidently infer that a nonzero fraction of events have positive χ eff , which requires at least one component to have spin tilt less than 90 · . A smaller fraction of events have negative χ eff , which requires at least one component with spin tilt > 90 · . A nonzero fraction of events have vanishingly small spins. The prior distributions on the parameters are marked with the non-filled histograms.  \n<!-- image -->  \np  \ninherit the orbital angular momenta of their progenitor systems, leading to significant spin magnitudes χ ≈ 0 . 7 for nearly equal-mass systems (Pretorius 2005; Berti & Volonteri 2008; Fishbach et al. 2017; Gerosa et al. 2018; Doctor et al. 2019; Rodriguez et al. 2018, 2019; Kimball et al. 2020a). If hierarchical mergers are present in GWTC-2, then one may expect correlations between the spins and masses of BBH systems, with more massive hierarchical mergers also possessing larger spins. We use the Multi Spin model to explore possible trends in the BBH spin distribution with mass, allowing for a distinct low-mass and high-mass subpopulation (with primary mass distributions parameterized by a power-law and Gaussian, respectively; see Sec. D.3), each with a distinct spin distribution. The low-mass power-law has a weak preference for smaller spins, as compared to the highmass Gaussian. Both subpopulations disfavor perfectly aligned systems, though the low-mass subpopulation has more support for small misalignments. In spite of these differences, the uncertainties on both of these subpopulations are broad enough that the two are fully consistent  \nwith each other, and we cannot confidently claim to detect a mass-dependence to the spin distribution at this stage. This is demonstrated in Fig. 13, which shows the posteriors for the spin distribution hyper-parameters associated with each mass subpopulation. These findings support the results of previous studies on GWTC-1, which could neither exclude nor confidently detect variation of the spin distribution with mass (Safarzadeh et al. 2020; Tiwari 2020).", '5.3. Merger rate and redshift evolution': "In this subsection we use the Power Law + Peak and Broken Power Law mass models with the Default spin model, and infer the merger rate using the NonEvolving redshift model and the Power-law evolution redshift model.  \nWe better constrain the binary black hole merger rate. Assuming a log-uniform prior, we find a BBH merger rate of R BBH = 23 . 9 +14 . 3 -8 . 6 Gpc -3 yr -1 using the Power Law + Peak mass distribution and the assumption of a non-evolving merger rate den-  \nFigure 13. Posterior distribution for Multi Spin spin hyper-parameters. Each subplot corresponds to a different parameter controlling the shape of the spin distribution (see Appendix D.3 for details). Each color corresponds to a different subpopulation (power-law or Gaussian) or component (primary mass or secondary mass) of the binary. The power-law subpopulation is slightly better measured than the Gaussian component, as a large number of detections are assigned to it. The results hint at a potential correlation between mass and spin, but the large measurement errors mean the spin distributions are consistent between the two subpopulations.  \n<!-- image -->  \nsity. We find that estimates of the BBH merger rate are robust to our choice of mass model, with excellent agreement between the Power Law + Peak , Broken Power Law , and Multi Peak models. The Truncated model yields a higher merger rate than the other models, but the results agree within statistical uncertainties: R BBH = 33 +22 -12 Gpc -3 yr -1 . Our BBH merger rate estimate is consistent with the hypothesis that a significant fraction of the merger rate is due to dynamically assembled binaries in globular clusters (Portegies Zwart & McMillan 2000; O'Leary et al. 2006; Moody & Sigurdsson 2009; Downing et al. 2011; Rodriguez et al. 2015; Park et al. 2017; Fragione & Kocsis 2018; Antonini & Gieles 2020), young/open star clusters (Banerjee et al. 2010; Ziosi et al. 2014), nuclear star clusters (O'Leary et al. 2009; Miller & Lauburg 2009; Antonini & Rasio 2016), or active galactic nuclei discs (McKernan et al. 2018; Tagawa et al. 2020; Gröbner et al. 2020).  \nIn Table 4, we show the merger rate in different primary mass bins, as inferred with the Broken Power Law , Power Law + Peak and Multi Peak models. Taking the range between the lowest 5% and highest 95% estimate across these three models in each mass bin, we find the merger rate to be ∼ 4 -14 Gpc -3 yr -1 in the range of 10 -20 M glyph[circledot] , ∼ 1 . 3 -5 . 3Gpc -3 yr -1 in the range of 20 -30 M glyph[circledot] , and ∼ 1 . 3 -5 . 2Gpc -3 yr -1 in the range of 30 -40 M glyph[circledot] . In Fig. 7, we showed that the primary mass spectrum turns over between 4 -10 M glyph[circledot] . We estimate the merger rate in this range to be ∼ 3 -21 Gpc -3 yr -1 .  \nOur estimate of the BBH merger rate include only systems with m 1 ≥ m 2 > 3 M glyph[circledot] , which notably excludes GW190814. If we calculate the merger rate for all systems down to m 2 ≥ 2 M glyph[circledot] using our models, thereby including GW190814, we infer a higher merger rate:  \nR BBH = 52 +52 -26 Gpc -3 yr -1 for the Power Law + Peak model. The reason for this change is that including GW190814 increases the low-mass rate (see Fig. 8b). However, because our mass distribution models do not extrapolate well to m 2 < 3 M glyph[circledot] (see Section 5.1), the fit with GW190814 likely overestimates the rate of systems with masses between ∼ 2 . 6 M glyph[circledot] and ∼ 6 M glyph[circledot] . Because of the uncertainty regarding the nature of GW190814 and the low significance of GW190426\\_152155 (the other NSBH candidate in GWTC-2), we do not attempt to model the NSBH mass distribution, and do not calculate an NSBH merger rate. An estimate of the merger rate for GW190814-like systems can be found in Abbott et al. (2020b).  \nWe update the binary neutron star merger rate. We give an update to the BNS rate based on the two confident BNS detections in GWTC-2, GW170817 and GW190425. We estimate a rate R BNS = 320 +490 -240 Gpc -3 yr -1 by Monte-Carlo sampling from a non-spinning distribution of systems, with component masses uniformly distributed between 1 M glyph[circledot] and 2 . 5 M glyph[circledot] , 9 and uniformly distributed in comoving volume and detector-frame time with isotropic orientations and sky locations. The detectability of each simulated system was approximated by requiring a network signal-to-noise ratio above 10 with signal-to-noise ratios above 5 in at least two detectors. Because of the longer observing time and the lack of additional detections, we find a slightly smaller value for the BNS rate than previously reported: R BNS = 320 +490 -240 Gpc -3 yr -1 . Assuming that there are 0 . 01 Milky Way equivalent galaxies (MWEG)  \nTable 4. Merger rate estimate in different primary mass bins. These results assume a non-evolving merger rate density, and exclude GW190814 from the analysis.  \nin 1 Mpc 3 (Kopparapu et al. 2008), this implies a rate of R BNS = 32 +49 -24 MWEG -1 Myr -1 .  \nThe BBH merger rate probably increases with redshift, but slower than the star-formation rate. Figure 14 shows the merger rate as a function of redshift using the Power law evolution model (see Appendix E for additional details, and Fig. 30 for a posterior predictive check). When we allow the merger rate to evolve with redshift according to (1+ z ) κ , we find that the z = 0 merger rate is R ( z = 0) = 19 . 3 +15 . 1 -9 . 0 Gpc -3 yr -1 . The posterior for the rate evolution parameter κ is shown in Fig. 15. Since GWTC-2 includes events with greater redshifts than the events in GWTC-1, we obtain a much tighter constraint on the evolution of the merger rate; compare our updated constraints of κ = 1 . 3 +2 . 1 -2 . 1 ( Power Law + Peak model) and κ = 1 . 8 +2 . 1 -2 . 2 ( Broken Power Law model) to the GWTC-1 result of κ = 8 . 4 +9 . 6 -9 . 5 . We find that the merger rate is consistent with a non-evolving distribution ( κ = 0 ), but is more likely to increase with increasing redshift, with κ > 0 at 85% credibility ( Power Law + Peak model) or 91% ( Broken Power Law model).  \nLocally ( z ≈ 0 ), the Madau-Dickinson star-formation rate (Madau & Dickinson 2014) corresponds to κ = 2 . 7 in our Power-Law Redshift parameterization. We infer κ < 2 . 7 at 86 % credibility with the Power Law + Peak mass model ( 77 % with Broken Power Law ). Another way of comparing our inferred merger rate to the star-formation rate is by looking at the ratio between the rate at z = 1 and z = 0 , R BBH ( z = 1) / R BBH ( z = 0) . For the star-formation rate R SFR ( z = 1) / R SFR ( z = 0) ≈ 6 , while for BBH systems, we infer R BBH ( z = 1) / R BBH ( z = 0) = 2 . 5 +8 . 0 -1 . 9 ( Power Law + Peak model). These results are consistent with most astrophysical formation channels, which predict a factor of ∼ 2 increase between the merger rate at z = 0 and z = 1 (Santoliquido et al. 2020; Dominik et al. 2013; Neijssel et al. 2019; Eldridge et al. 2019; Mapelli et al. 2017; Baibhav et al. 2019; Rodriguez & Loeb 2018).", '6. CONCLUSIONS': "The publication of the second LIGO-Virgo gravitationalwave transient catalog has increased the population of BBH events by a factor of more than four. The new  \nFigure 14. Merger rate density as a function of redshift, fit to the Power-law Evolution model. The solid curve shows the median rate density, while the dark (light) shaded region shows 50% (90%) credible intervals. The dashed curve shows the shape of the SFR. The data exhibit a mild preference for the merger rate to increase with redshift, but are consistent with a flat distribution as well as one that tracks the SFR.  \n<!-- image -->  \nFigure 15. Posterior for the redshift evolution parameter κ from the Power-law Evolution model, which assumes that rate density scales like (1 + z ) κ . We assume the Power Law + Peak and Broken Power Law mass models, and take a flat prior on κ .  \n<!-- image -->  \ncatalog has highlighted the limitations of some early population models while yielding remarkable new signatures:  \n- 1. We find that the BBH primary mass spectrum is not well-described as a simple power-law with an abrupt cut-off; there is a strong statistical preference for other models with non-trivial features such as a peak or a tapering. These features occur at ≈ 37M glyph[circledot] , where one might expect pair instability supernovae (and pulsational pair-instability supernovae) to shape the mass distribution of BHs.\n- 2. At the opposite end of the spectrum, we observe a dearth of systems between NS and BH masses, suggesting that the BH mass spectrum likely turns over at ∼ 7 . 8 +1 . 8 -2 . 0 M glyph[circledot] . We constrain the minimum mass of BHs in BBH systems to be m min glyph[lessorsimilar] 5 . 7 M glyph[circledot] at 90% credibility. This is greater than the mass of BH candidates in Galactic binaries, e.g., Thompson et al. (2019). These results hold only when we restrict our analysis to events with both component masses above 3 M glyph[circledot] .\n- 3. Meanwhile, we find that our models fail to fit GW190814 together with the BBH systems with both components above 3 M glyph[circledot] . This may indicate that GW190814 belongs to a distinct population, or that there are additional features at the low mass end of the BBH mass spectrum that are missing in our models. This is perhaps unsurprising as the combination of mass ratio, merger rate, and secondary mass inferred from this system pose a challenge to our current understanding of compact binary formation (Abbott et al. 2020b; Arca Sedda 2020b; Zevin et al. 2020; Safarzadeh 2020).\n- 4. We detect clear evidence of spin-induced, general relativistic precession of the orbital plane. We determine that this signature is not due to a single precessing merger, but from the overall preference of the data for precessing waveforms.\n- 5. We observe that some fraction of the BHs in GWTC-2 are spinning with an orientation that is anti-aligned with respect to the orbital angular momentum of the binary. If we plausibly assume that all binaries with anti-aligned spins are assembled dynamically, this may imply that LIGO-Virgo events merge both dynamically and in the field. Based on the inferred mass and spin distributions, we find no clear evidence for or against hierarchical mergers in GWTC-2.\n- 6. We compute the rate of compact binary mergers, finding R BNS = 320 +490 -240 Gpc -3 yr -1 and R BBH = 23 . 9 +14 . 3 -8 . 6 Gpc -3 yr -1 . The data are consistent with both a merger rate that is constant in time and one that tracks the SFR in the local universe, though the data prefer a merger rate that is somewhere in between. We find that the merger rate at z = 1 differs from the merger rate at z = 0 by a factor of R BBH ( z = 1) / R BBH ( z = 0) = 2 . 5 +8 . 0 -1 . 9 , to be compared with the SFR, R SFR ( z = 1) / R SFR ( z = 0) ∼ 6 .  \nWhile a clearer picture is emerging of the population properties of compact binaries, key questions remain. How do we best characterize the deviations from powerlaw in the primary BH mass spectrum, and what is the physical origin of these new features? What is the origin of BBH mergers in the high-mass gap: hierarchical mergers, stars producing remnants heavier than expected from pair instability supernovae theory, or something else? What is the shape of the mass spectrum between NS and BH masses, and does the current dearth of systems between ∼ 3 M glyph[circledot] and ∼ 6 M glyph[circledot] represent an empty low-mass gap? If so, do systems like the secondary mass of GW190814 belong to the NS or BH side of the gap? Is the observation of anti-aligned spins indicative of dynamically assembled binaries? As the sensitivity of LIGO, Virgo, and KAGRA improves, and as more gravitational-wave transients are detected, we expect to begin to answer these questions. As future observations subject our models to increasing scrutiny, it is inevitable that refinements will be required to fit newly resolved features. This cycle of refining models to account for new data will reveal new questions while providing an evolving understanding of the conclusions presented here.  \nThe authors gratefully acknowledge the support of the United States National Science Foundation (NSF) for the construction and operation of the LIGO Laboratory and Advanced LIGO as well as the Science and Technology Facilities Council (STFC) of the United Kingdom, the Max-Planck-Society (MPS), and the State of Niedersachsen/Germany for support of the construction of Advanced LIGO and construction and operation of the GEO600 detector. Additional support for Advanced LIGO was provided by the Australian Research Council. The authors gratefully acknowledge the Italian Istituto Nazionale di Fisica Nucleare (INFN), the French Centre National de la Recherche Scientifique (CNRS) and the Netherlands Organization for Scientific Research, for the construction and operation of the Virgo detector and the creation and support of the EGO consortium. The authors also gratefully acknowledge research support  \nfrom these agencies as well as by the Council of Scientific and Industrial Research of India, the Department of Science and Technology, India, the Science & Engineering Research Board (SERB), India, the Ministry of Human Resource Development, India, the Spanish Agencia Estatal de Investigación, the Vicepresidència i Conselleria d'Innovació, Recerca i Turisme and the Conselleria d'Educació i Universitat del Govern de les Illes Balears, the Conselleria d'Innovació, Universitats, Ciència i Societat Digital de la Generalitat Valenciana and the CERCA Programme Generalitat de Catalunya, Spain, the National Science Centre of Poland, the Swiss National Science Foundation (SNSF), the Russian Foundation for Basic Research, the Russian Science Foundation, the European Commission, the European Regional Development Funds (ERDF), the Royal Society, the Scottish Funding Council, the Scottish Universities Physics Alliance, the Hungarian Scientific Research Fund (OTKA), the French Lyon Institute of Origins (LIO), the Belgian Fonds de la Recherche Scientifique (FRS-FNRS), Actions de Recherche Concertées (ARC) and Fonds Wetenschappelijk Onderzoek - Vlaanderen (FWO), Belgium,  \nthe Paris Île-de-France Region, the National Research, Development and Innovation Office Hungary (NKFIH), the National Research Foundation of Korea, Industry Canada and the Province of Ontario through the Ministry of Economic Development and Innovation, the Natural Science and Engineering Research Council Canada, the Canadian Institute for Advanced Research, the Brazilian Ministry of Science, Technology, Innovations, and Communications, the International Center for Theoretical Physics South American Institute for Fundamental Research (ICTP-SAIFR), the Research Grants Council of Hong Kong, the National Natural Science Foundation of China (NSFC), the Leverhulme Trust, the Research Corporation, the Ministry of Science and Technology (MOST), Taiwan and the Kavli Foundation. The authors gratefully acknowledge the support of the NSF, STFC, INFN and CNRS for provision of computational resources.  \nWe would like to thank all of the essential workers who put their health at risk during the COVID-19 pandemic, without whom we would not have been able to complete this work.  \nThis is LIGO document number LIGO-P2000077.", 'REFERENCES': 'Aasi, J., et al. 2015, CQGra, 32, 074001 Abbott, B., Abbott, R., Abbott, T., et al. 2013, Prospects for Observing and Localizing Gravitational-Wave Transients with Advanced LIGO, Advanced Virgo and KAGRA, arXiv:1304.0670  \n- -. 2016a, PhRvD, 94, 064035\n- -. 2016b, PhRvL, 116, 241103\n- -. 2016c, PhRvL, 116, 241102\n- -. 2019a, ApJL, 882, L24\n- -. 2019b, PhRvX, 9, 031040\n- -. 2019c, arXiv:1912.11716  \nAbbott, R., Abbott, T., Abraham, S., et al. 2020a, arXiv:2004.08342  \n- -. 2020b, ApJL, 896, L44\n- -. 2020c, arXiv e-prints, arXiv:2010.14527  \nAbbott, R., et al. 2020d  \n- -. 2020e, PhRvL, 125, 101102\n- -. 2020f, ApJL, 900, L13  \nAcernese, F., et al. 2015, CQGra, 32, 024001  \nAdams, T., Buskulic, D., Germain, V., et al. 2016, Classical and Quantum Gravity, 33, 175012  \nAllen, B. 2005, PhRvD, 71, "062001"  \nAllen, B., Anderson, W. G., Brady, P. R., Brown, D. A., & Creighton, J. D. E. 2012, PhRvD, 85, "122006"  \nAntonini, F., & Gieles, M. 2020, arXiv e-prints, arXiv:2009.01861 Antonini, F., Murray, N., & Mikkola, S. 2014, ApJ, 781, 45 Antonini, F., & Rasio, F. A. 2016, ApJ, 831, 187 Antonini, F., Rodriguez, C. L., Petrovich, C., & Fischer, C. L. 2018, Mon. Not. R. Ast. Soc. Lett., 480, L58 Antonini, F., Toonen, S., & Hamers, A. S. 2017, ApJ, 841, 77 Apostolatos, T. A., Cutler, C., Sussman, G. J., & Thorne, K. S. 1994, PhRvD, 49, 6274 Arca Sedda, M. 2020a, ApJ, 891, 47 -. 2020b, Communications Physics, 3, 43 Arca Sedda, M., Mapelli, M., Spera, M., Benacquista, M., & Giacobbo, N. 2020, ApJ, 894, 133 Ashton, G., & Thrane, E. 2020, Mon. Not. Roy. Astron. Soc., 498, 1905 Ashton, G., Hübner, M., Lasky, P. D., et al. 2019, ApJS, 241, 27 Babak, S., Taracchini, A., & Buonanno, A. 2017, Phys. Rev. D, 95, 024010 Baibhav, V., Berti, E., Gerosa, D., et al. 2019, PhRvD, 100, 064060 Bailyn, C., Jain, R. K., Coppi, P., & Orosz, J. A. 1998, ApJ, 499, 367 Banerjee, S., Baumgardt, H., & Kroupa, P. 2010, MNRAS, 402, 371', 'A. ESTIMATING THE DETECTION FRACTION': 'A key ingredient in Eqs. (1) and (2) is the detection fraction ξ (Λ) , the fraction of systems within some prior volume (redshift z < 2 . 3 ) that we expect to successfully detect. The detection fraction quantifies selection biases, and so it is critical to accurately characterize. For a population described by hyper-parameters Λ , the detection fraction is  \nξ (Λ) = ∫ P det ( θ ) π ( θ | Λ) dθ. (A1)  \nHere, P det ( θ ) is the detection probability: the probability that an event with parameters θ is detectable. The detection probability depends primarily on the masses and redshift of a system, and, to a lesser degree, on the spins.  \nWe calculate ξ (Λ) using injections. We simulate compact binary signals from a reference population and record which ones are successfully detected by the PyCBC and GstLAL search pipelines; see Abbott et al. (2020c). Following Tiwari (2018); Farr (2019); Vitale (2020); Loredo (2004), the point estimate for Eq. A1 is calculated using a Monte Carlo integral over found injections:  \nˆ ξ (Λ) = 1 N inj N found ∑ j =1 π ( θ j | Λ) p draw ( θ j ) , (A2)  \nwhere N inj is the total number of injections, N found are the injections that are successfully detected, and p draw is the probability distribution from which the injections are drawn; see LIGO-Virgo (2020) for additional details. When sampling the population likelihood, we marginalize over the uncertainty in ˆ ξ (Λ) following Farr (2019), and ensure that the effective number of found injections remaining after population re-weighting is sufficiently high ( N eff > 4 N det ).  \nFor the O3a observing period, we use the injection campaign described in Abbott et al. (2020c) and characterize the found injections as those recovered with a FAR below our threshold of 1 yr -1 in either PyCBC or GstLAL . For the O1 and O2 observing period, we supplement the O3a pipeline injections with mock injections drawn from the same distribution p draw above. For the mock injections, we calculate P det ( m 1 , m 2 , z, χ 1 ,z , χ 2 ,z ) according to the semi-analytic approximation described in Abbott et al. (2019a), based on a single-detector signal-to-noise ratio threshold ρ = 8 and the Advanced LIGO Early-High Noise PSD (Abbott et al. 2013). We combine O1, O2 and O3 injection sets ensuring a constant rate of injections across the total observing time, yielding N inj ≈ 7 . 7 × 10 7 injections for O3a and N inj ≈ 7 . 1 × 10 7 for O1 and O2. To control computational costs, not all of the injections are performed in real data. Before injecting, the expected network signal-to-noise ratio of the injections is computed, and the hopeless injections with signal-to-noise ratio < 6 are removed.  \nDue to the finite number of injections, we approximate Eq. (A1) with a fixed spin distribution instead of the distribution implied by Λ . When combining the Truncated , Power Law + Peak , Broken Power Law and Multi Peak mass models together with the Default spin distribution, we assume that the aligned spin components χ 1 ,z , χ 2 ,z are independently drawn from a uniform distribution U ( -0 . 5 , 0 . 5) . By making this approximation, we are in effect ignoring selection effects due to spin. Nevertheless, we expect this approximation to have a negligible impact on the inferred spin distribution compared to the statistical uncertainties. For aligned spin components in the range ( -0 . 5 , 0 . 5) , the detection probability varies by no more than a factor of 2 (Ng et al. 2018a). Furthermore, our main conclusions regarding the spin distribution inferred from the Default model are supported by the Gaussian model, which requires no approximations for spin selection effects. The Multi Spin model calculates Eq. (A1) by calibrating a semi-analytic approximation to the list of found injections (Wysocki 2020).  \nThe transfer function between the observed strain and astrophysical strain is subject to a systematic calibration uncertainty. We neglect this calibration uncertainty in our estimates of the search sensitivity above. For the O3a observing run, the amplitude uncertainty was glyph[lessorsimilar] 3% (Sun et al. 2020), which leads to a glyph[lessorsimilar] 10% systematic uncertainty in the sensitive spacetime volume and the inferred merger rate. This systematic uncertainty is subdominant to our uncertainties from Poisson counting error.', 'B. DETAILS OF MASS POPULATION MODELS': 'In this section we provide details about the population models described above in Section 3; see also Fig. 1. Each subsection includes a table with a summary of the parameters for that model and the prior distribution used for each parameter. The prior distributions are indicated using abbreviations: for example, U(0 , 1) translates to uniform on the interval (0 , 1) and LU( 10 -6 , 10 5 ) translates to log-uniform on the interval 10 -6 , 10 5 .  \nTable 5. Summary of Truncated parameters.', 'B.1. Truncated mass model': "This model is equivalent to 'Model B' in Abbott et al. (2019a). The primary mass distribution for this model follows a power-law with spectral index α , and with a sharp cut-off at the lower end m min and the upper end of the distribution m max :  \nπ ( m 1 | α, m min , m max ) ∝    m -α 1 m min < m 1 < m max 0 otherwise , (B3)  \nMeanwhile, the mass ratio q ≡ m 2 /m 1 follows a power-law distribution with spectral index β q  \nπ ( q | β q , m min , m 1 ) ∝    q β q m min < m 2 < m 1 0 otherwise . (B4)  \nThe hyper-parameters for this model are summarized in Table 5.", 'B.2. Power Law + Peak mass model': "This is equivalent to 'Model C' from Abbott et al. (2019a). It is motivated by the idea that the mass loss undergone by pulsational pair-instability supernovae could lead to a pile-up of BBH events before the pair-instability gap (Talbot & Thrane 2018). The primary mass distribution is an extension of Truncated with the addition of tapering at the lower mass end of the distribution and a Gaussian component:  \nπ ( m 1 | λ peak , α, m min , δ m , m max , µ m , σ m ) = [ (1 -λ peak ) P ( m 1 | -α, m max ) + λ peak G ( m 1 | µ m , σ m ) ] S ( m 1 | m min , δ m ) . (B5)  \nHere, P ( m 1 | -α, m max ) is a normalized power-law distribution with spectral index -α and high-mass cut-off m max . Meanwhile, G ( m 1 | µ m , σ m ) is a normalized Gaussian distribution with mean µ m and width σ m . The parameter λ peak is a mixing fraction determining the relative prevalence of mergers in P and G . Finally, S ( m 1 , m min , δ m ) is a smoothing function, which rises from 0 to 1 over the interval ( m min , m min + δ m ) :  \nS ( m | m min , δ m ) =        0 ( m<m min ) [ f ( m -m min , δ m ) + 1] -1 ( m min ≤ m<m min + δ m ) 1 ( m ≥ m min + δ m ) (B6)  \nwith  \nThe conditional mass ratio distribution in this model also includes the smoothing term:  \nf ( m ' , δ m ) = exp ( δ m m ' + δ m m ' -δ m ) . (B7)  \nπ ( q | β, m 1 , m min , δ m ) ∝ q β q S ( qm 1 | m min , δ m ) . (B8)  \nThe hyper-parameters for this model are summarized in Table 6.  \nIn Fig. 16, we provide a corner plot representation of the posterior distribution for the Power Law + Peak hyper-parameters. The ( µ m , λ peak ) panel describes the Gaussian component: µ m is the center of the Gaussian while  \nTable 6. The parameters that describe the BBH mass distribution for Model Power Law + Peak .  \nλ peak is the fraction of mergers taking place in the Gaussian (as opposed to the power-law distribution). Judging from this panel, it appears at first that λ peak peaks close to 0 (corresponding to no Gaussian peak). However, if we zoom in as in Fig. 6, we see that the posterior for λ peak is peaked clearly away from zero at ∼ 0 . 02 . This is consistent with the large Bayes factor indicating preference for Power Law + Peak over Truncated .", 'B.3. Broken Power Law mass model': 'This model is an extension of Truncated . The primary mass distribution consists of a broken power-law. This is motivated by the potential tapering of the primary mass distribution at high masses. Also, the model employs a smoothing function to prevent a sharp cut-off at low masses.  \nπ ( m 1 | α 1 , α 2 , m min , m max ) ∝        m -α 1 1 S ( m 1 | m min , δ m ) m min < m 1 < m break m -α 2 1 S ( m 1 | m min , δ m ) m break < m 1 < m max 0 otherwise . (B9)  \nHere,  \nm break = m min + b ( m max -m min ) , (B10)  \nis the mass where there is a break in the spectral index and b is the fraction of the way between m min and m max at which the primary mass distribution undergoes a break. Meanwhile, S ( m 1 , m min , δ m ) is a smoothing function as in Eq. (B6). The conditional mass ratio distribution is the same as in the Power Law + Peak model; see Eq. B8. The hyper-parameters for this model are summarized in Table 7. In Fig. 17 we provide a corner plot for Broken Power Law . In the limit of no low-mass smoothing ( δ m = 0 ), and in the limit of a second power-law with a steep slope that mimics a sharp cutoff ( m break = m max ), this model reduces to Truncated . Above, we noted that the Broken Power Law model prefers a break in the primary mass spectrum near 40 M glyph[circledot] . On the other hand, if we believe that the feature represented by m break should be closer to a sharp cutoff, then the cut-off must occur at higher masses approaching the maximum mass of Truncated at m max = 74 . 6 +15 . 4 -8 . 6 M glyph[circledot] . This can be seen by the correlation between b and α 2 in Fig. 17.', 'B.4. Multi Peak mass model': 'This model in an extension of Power Law + Peak , where there is an additional Gaussian component at the upper end of the mass distribution motivated by a possible subpopulation of objects in the upper mass gap:  \nπ ( m 1 | λ, α, m min , δ m , m max , µ m , σ m ) = [ (1 -λ ) P ( m 1 | -α, m max ) + λλ 1 G ( m 1 | µ m, 1 , σ m, 1 ) + λ (1 -λ 1 ) G ( m 1 | µ m, 2 , σ m, 2 ) ] S ( m 1 | m min , δ m ) . (B11)  \nHere, the parameters λ and λ 1 correspond to the fraction of binaries in any Gaussian component and the fraction of binaries in the lower-mass Gaussian of the Gaussian components, respectively. The distribution G ( m 1 | µ m, 1 , σ m, 1 ) is a  \n10  \nFigure 16. Posterior distribution for mass hyper-parameters for Power Law + Peak . The fit excludes GW190814. The contours represent 50% and 90% credible bounds.  \n<!-- image -->  \nnormalized Gaussian distribution for the lower-mass peak with mean µ m, 1 and width σ m, 2 and G ( m 1 | µ m, 2 , σ m, 1 ) is a normalized Gaussian distribution for the upper-mass peak with mean µ m, 2 and width σ m, 2 . The hyper-parameters for this model are summarized in Table 8. In Fig. 18 we provide a corner plot for Multi Peak for parameters corresponding to the two Gaussian peaks. The mean of the upper-mass peak µ m, 2 = 68 +18 -14 M glyph[circledot] is located at approximately twice the mean of the lower-mass peak µ m, 1 = 33 . 4 +4 . 4 -4 . 9 M glyph[circledot] . The remaining parameters are: α = 2 . 9 +1 . 9 -1 . 4 , β q = 0 . 9 +1 . 9 -1 . 3 , m min = 4 . 6 +1 . 3 -1 . 8 M glyph[circledot] and m max = 65 +31 -30 M glyph[circledot] . In Fig. 19, we provide a posterior predictive check for all of the mass models used in this analysis.  \nTable 7. Summary of Broken Power Law parameters.  \nTable 8. Parameters for the BBH mass distribution for Model Multi Peak .', 'C. MASS MODEL CHECKING': 'Section 5.1 describes the inferred mass distribution obtained with the Truncated , Broken Power Law , Power Law + Peak , and Multi Peak models, and compares the different models by calculating their Bayes factors. Here we assess the goodness-of-fit of the models using posterior predictive checks, comparing predicted and empirical catalogs of observed m 1 distributions in Fig. 19. The light colored bands show the cumulative distribution of m 1 as predicted by the model, while the darker bands show the empirical distribution based on the actual events observed in GWTC-2. The bands represent a family of curves, where each curve corresponds to a different draw from the population hyper-posterior. Each draw from the hyper-posterior updates both the predicted distribution (in the lighter color) and the empirical distribution (in the darker color), as the individual event posteriors are updated according to the inferred population distribution (Fishbach et al. 2020b; Galaudage et al. 2019; Miller et al. 2020). If the model is a good fit to the data, the dark colored bands should overlap with the light colored bands. Figure 19 shows the relatively poor fit for the Truncated model, which cannot capture the excess of events at ∼ 30 -40 M glyph[circledot] compared to glyph[greaterorsimilar] 40 M glyph[circledot] . The remaining panels show the improved fits with the Power Law + Peak , Broken Power Law and Multi Peak models. These results are consistent with the Bayes factors in Table 2, which conclude that the Truncated model is disfavored by a Bayes factor of 10-80 relative to the other models.', 'C.1. On GW190412': 'Other than GW190814, we find that GW190412 (Abbott et al. 2020a), when analyzed with a population informed prior, remains the only system for which we can confidently bound the mass ratio away from unity, yielding q < 0 . 53 at 99% credibility (using the Power Law + Peak mass model). All other events, when analyzed with a population-informed prior, are consistent with q = 1 at 99% credibility. Repeating the analysis in Abbott et al. (2020a), we perform a leave-one-out analysis without GW190412 and find β q = 4 . 0 +6 . 4 -3 . 2 ( β q = 4 . 5 +5 . 9 -3 . 5 ) for the Power Law + Peak ( Broken Power Law ) model. The β q posterior inferred with the inclusion of GW190412 ( β q = 1 . 3 +2 . 4 -1 . 5 for the Power Law +  \nFigure 17. Posterior distribution for mass hyper-parameters for Broken Power Law . The fit excludes GW190814. The contours represent 50% , 90% credible bounds.  \n<!-- image -->  \nPeak model; β q = 1 . 4 +2 . 5 -1 . 5 for the Broken Power Law model) has moderate ( ∼ 50% ) overlap with the leave-one-out β q posterior, indicating that, consistent with the conclusion in Abbott et al. (2020a), GW190412 likely belongs to the low mass ratio tail of the distribution rather than a distinct subpopulation of asymmetric systems.', 'C.2. On GW190521': 'As discussed in Section 5.1, the most massive event, GW190521 (Abbott et al. 2020e), is an outlier with respect to the Truncated model (see Fig. 2), but fits well within the mass distributions inferred from the other models. In Fig. 20, we show the effect of GW190521 on the primary mass distribution. This event shifts the best-fit mass distribution, but this shift is within the statistical uncertainties. Thus, we find no evidence that GW190521 is an outlier within the framework of the Power Law + Peak and Broken Power Law mass models. This finding is supported by the posterior predictive check in subsection C.4. In Fig. 21, we show how the primary mass posterior distribution for GW190521 changes when we use the Power Law + Peak model to inform our prior. While the population-informed posterior on the primary mass prefers smaller masses (Fishbach et al. 2020b), the conclusion that the primary mass of  \nFigure 18. Posterior distribution for mass hyper-parameters for Multi Peak . The fit excludes GW190814. The contours represent 50% and 90% credible bounds.  \n<!-- image -->  \nGW190521 is above 67 M glyph[circledot] (99% credibility) is robust to the choice of prior, consistent with the claims in Abbott et al. (2020e).', 'C.3. On GW190814': 'On the other hand, we see clear indication that GW190814 is an outlier with respect to the BBH population within the framework of the Power Law + Peak and Broken Power Law mass models, as discussed in Section 5.1. As an additional posterior predictive check, following the analysis described in Fishbach et al. (2020b) and Abbott et al. (2020a), we use the posterior predictive distribution, inferred without GW190814, to construct a distribution for the  \nFigure 20. A comparison of the primary BH mass distribution for the population with and without GW190521. The data are fit using the Power Law + Peak model; the Broken Power Law model produces similar results. The solid curves are the posterior predictive distributions while the shaded regions show the 90% credible interval. The inclusion/exclusion of GW190521 does not have a significant effect on the fit. GW190814 is excluded from this analysis.  \n<!-- image -->  \n/circledot  \n/circledot  \nFigure 19. A posterior predictive check: the cumulative density function (CDF) of the observed primary mass distribution for the Truncated , Power Law + Peak , Broken Power Law and Multi Peak models. The observed event distribution is shown in the darker colors. The thickness of the bands indicates the 90% credibility range. The Power Law + Peak , Broken Power Law and Multi Peak models are a better fit than the Truncated model; the dark band overlaps entirely with the light colored band. This is due to the Power Law + Peak , Broken Power Law and Multi Peak models having more flexibility to fit the relative excess of binaries in the 30 M glyph[circledot] -40 M glyph[circledot] region compared to the glyph[greaterorsimilar] 40 M glyph[circledot] region. GW190814 is excluded from this analysis.  \n<!-- image -->  \n/circledot  \nFigure 21. The posterior probability density for the primary mass of GW190521 using the original default prior (blue; flat in redshifted masses) and a reweighted version (green) obtained by using the Power Law + Peak mass model. Results using the Broken Power Law mass model are similar. The reweighted version shifts the posterior support to lower masses, with m 1 < 83 M glyph[circledot] (99% credibility).  \n<!-- image -->  \n/circledot  \nminimum m 2 detected in a sample of 45 events. When using both the Power Law + Peak and Broken Power Law models, we find that the observation of a system with a secondary mass equal to or smaller than the that of GW190814 ( 2 . 59 +0 . 08 -0 . 09 M glyph[circledot] ) is highly improbable, with probability < 0 . 02% for both Power Law + Peak and Broken Power Law ; see Fig. 22b for the distribution of the minimum observed secondary mass in a sample of 45 events predicted by the Power Law + Peak model. The distribution for Broken Power Law is qualitatively similar. The mass ratio of GW190814 is also somewhat unusual according to this posterior predictive check; see Fig. 22a. Observing an event with the mass ratio of GW190814 or smaller, based on the fit to the other 44 BBH events, has probability < 0 . 02% in both the Power Law + Peak and Broken Power Law models. These posterior predictive checks suggest that GW190814 is not a typical BBH, and support the conclusion that there may be a dearth of systems between ∼ 2 . 6 M glyph[circledot] and ∼ 6 M glyph[circledot] . Future observations will reveal the precise shape of the mass distribution at low masses and extreme mass ratios, and better determine the nature of GW190814.', 'C.4. Mass and distance checks with a burst analysis': 'The earlier posterior predictive checks in this section compared simulated sets of BBH masses to the observed set of catalog events. As a complimentary posterior predictive check, we can simulate the gravitational-wave signals from these predicted events, run it through our search pipelines, and compare the synthetic data, as detected by the pipeline, to the observed data. As a proof of principle, we carry out this posterior predictive check with the Coherent WaveBurst ( cWB ) pipeline (Klimenko & Mitselmakher 2004; Klimenko et al. 2016), which is designed to detect unmodeled gravitational-wave transients. The cWB analysis resulted in the detection of 22 BBH events in GWTC-2. We investigate whether the set of cWB observations is consistent with the model predictions. We focus on assessing possible outliers at high masses and high redshifts, where cWB is especially sensitive to BBH signals. In particular, we examine whether GW190521 ,which was recovered with higher significance by the cWB saerch than the templated searches, is an outlier in the context of our mass and redshift models. Following Klimenko et al. (2016), we calculate the expected distribution of the central frequency f (which depends on the redshifted mass of the BBH) and coherent signal-to-noise ratio ρ (which, for a given redshifted mass, depends on the distance of the BBH) for two different population models: Power Law + Peak and Broken Power Law , using the Non-Evolving redshift model. We then compare the empirical distribution of ( f, ρ ) , as recovered by the cWB pipeline to the distribution predicted by the population model. We quantify the comparison by calculating a p -value for each event i , which measures how unusual its observed ( f i , ρ i ) is, given the distribution of predicted ( f, ρ ) . The central frequency is a proxy for the redshifted mass while the signal-to-noise ratio is a proxy for the distance.  \nTo compute the predicted distribution of ( f, ρ ) , we inject simulated waveforms drawn from the Power Law + Peak and Broken Power Law distributions into the O1, O2 and O3a data and compile injections recovered by cWB  \nFigure 22. Left: the posterior population distribution for primary and secondary mass with and without GW190814. Shown here are the 99% credible intervals. Dark blue is with GW190814 and light blue is without. The shaded regions show the astrophysical , while the colored contours show the observed distribution (as it appears in the catalog due to selection effects). The median values of the posterior of GW190814 is marked with a star. Right: Distribution of the minimum secondary mass detected out of 45 detections, predicted from the fit to the Power Law + Peak model to the BBH population excluding GW190814 (light blue) and the BBH population including GW190814 (dark blue). The dashed line and shaded region (gray) denote the median and 90% symmetric credible interval on the secondary mass of GW190814. This distribution is qualitatively similar to the distribution predicted from the fit to the Broken Power Law model.  \n<!-- image -->  \n/circledot  \n- (a) The model-averaged astrophysical mass distribution p ( m 1 , m 2 ) .  \n<!-- image -->  \n/circledot  \n(b) Observed minimum secondary mass distribution.  \nwith a FAR < 1 yr -1 . The central frequencies and coherent signal-to-noise ratios of the recovered injections generated according to Power Law + Peak are plotted in Fig. 23. The results for Broken Power Law are similar. The locations of the 22 detections on this plane are visually consistent with the model predictions, indicating that the model is a reasonably good fit to the data. The event with the lowest p -value (least consistent with predictions) is GW190521, with p -values of 0.053 and 0.077 for the Broken Power Law model and the Power Law + Peak model respectively. These p -values indicate that GW190521 is a moderately unusual detection, but it is consistent with the population models.', 'D. DETAILS OF SPIN POPULATION MODELS': "D.1. Default spin model  \nThis model was introduced in Abbott et al. (2019a). Following Wysocki et al. (2019a), the dimensionless spin magnitude distribution is taken to be a Beta distribution,  \nπ ( χ 1 , 2 | α χ , β χ ) = Beta ( α χ , β χ ) , (D12)  \nwhere α χ and β χ are the standard shape parameters that determine the distribution's mean and variance. The Beta distribution is convenient because it is bounded on (0,1). The distributions for χ 1 and χ 2 are assumed to be the same. Following Talbot & Thrane (2017), we define z = cos θ 1 , 2 as the cosine of the tilt angle between component spin and a binary's orbital angular momentum, and assume that z is distributed as a mixture of two populations:  \nπ ( z | ζ, σ t ) = ζ G t ( z | σ t ) + (1 -ζ ) I ( z ) . (D13)  \nHere, I ( z ) is an isotropic distribution, while G t ( z | σ t ) is a truncated Gaussian, peaking at z = 0 (perfect alignment) with width σ t . The mixing parameter ζ controls the relative fraction of mergers drawn from the isotropic distribution and Gaussian subpopulations. The isotropic subpopulation is intended to accommodate dynamically assembled binaries, while G t is a model for field mergers. The hyper-parameters for this model and their priors are summarized in Table 9. Additional constraints to the priors on µ χ and σ 2 χ are applied by setting α χ , β χ > 1 .  \nIn Fig. 24 we provide a corner plot for the Default spin model. This model prefers modest spin magnitudes. It favors the hypothesis that binaries are preferentially aligned ( ζ → 1 and σ t glyph[lessorsimilar] 2 ), albeit with potentially large misalignment  \n<!-- image -->  \n˜{√√∐⌈ ˜√˜√√˜{̂glyph[arrowbt] ⋃∮︷⊎  \nFigure 23. The coherent signal-to-noise ratio and central frequency for 22 BBH events detected by one of our detection pipelines, cWB , in O1, O2 and O3a (violet dots) compared to simulated BBH events from the Power Law + Peak model detected by cWB . The consistency between the distributions of simulated and observed triggers shows that the Power Law + Peak mass model coupled with the non-evolving redshift distribution is a good fit to the data. Results for the Broken Power Law model are similar.Table 9. Summary of Default spin parameters.  \nangles ( σ t > 0 ). The case of perfect alignment, which would correspond to ζ = 1 and σ = 0 , is disfavored, lying outside the 99% credible bound on ζ and σ t . Within the main text, Fig. 10 shows the implied distributions of component spin magnitudes and tilt angles. The implied distributions of the effective precession spin parameter ( χ p ) and the effective inspiral spin parameter ( χ eff ) are shown in Figs. 9 and 11, respectively; these distributions are in good agreement with the results obtained using the Gaussian model described below. In particular, both models predict the existence of systems with anti-aligned spins (negative χ eff ), and in-plane spin components (non-zero χ p ).", 'D.2. Gaussian spin model': "The Gaussian spin model offers an alternative description of BBH spins. It is convenient to measure the distribution of the effective inspiral spin parameter ( χ eff ) and the effective precession spin parameter ( χ p ), which are better constrained than individual component spin magnitudes or tilts. We parameterize the distributions of χ eff and χ p , using a bivariate Gaussian:  \nπ ( χ eff , χ p | µ eff , σ eff , µ p , σ p , ρ ) ∝ G ( χ eff , χ p | µµµ, ΣΣΣ) . (D14)  \nThe distribution has a mean µµµ = ( µ eff , µ p ) and a covariance matrix  \nΣΣΣ = ( σ 2 eff ρσ eff σ p ρσ eff σ p σ 2 p ) . (D15)  \nThe population parameters appearing in Eqs. (D14) and (D15) and their associated priors are summarized in Table 10. We truncate and normalize Eq. (D14) based on the allowed regions of the effective inspiral spin parameter: χ eff ∈ ( -1 , 1) and χ p ∈ (0 , 1) . The results from the Gaussian model are obtained assuming a Truncated mass model with α = -2 . 2 , β q = 1 . 3 , m min = 5 M glyph[circledot] , and m max = 75 M glyph[circledot] , consistent with the median values obtained when fitting the Truncated  \n√  \n}√{  \nFigure 24. Posterior distribution for spin hyper-parameters for Default , assuming the Power Law + Peak mass model and Non-Evolving redshift distribution. The fit excludes GW190814. The contours represent 50% , 90% credible bounds. A perfectly aligned spin distribution ( σ t = 0 , ζ = 1 ) is ruled out at > 99% credibility, consistent with the results of the Gaussian model, but the data disfavor a purely isotropic distribution ( ζ = 0 or σ t glyph[greaterorsimilar] 2 ).  \n<!-- image -->  \nTable 10. Summary of Gaussian spin parameters.  \n-  \nmodel to GWTC-2. We additionally assume a comoving merger rate density that grows as (1 + z ) 2 . 7 . Although the Truncated model is disfavored relative to the more complex mass models discussed above, it is sufficient for purposes of constructing an informed mass ratio distribution, the primary confounding factor in efforts to measure χ eff and χ p (Ng et al. 2018a).  \nThe full posterior on parameters of the Gaussian model is shown in Fig. 25. We find no correlation between the parameters of the χ eff and χ p distributions, nor do we obtain any information regarding the degree of correlation ρ between the effective inspiral spin parameter and the effective precession spin parameter. As discussed in Sec. 5.2,  \nρ  \nFigure 25. Posterior distribution for spin hyper-parameters under the Gaussian model. The fit again excludes GW190814, and contours represent 50% and 90% credible bounds. Consistent with the results of the Default model, which results out a perfectly aligned spin distribution at 99% credibility, here we find that a vanishing χ p distribution ( µ p = σ p = 0 ) is disfavored.  \n<!-- image -->  \nanalysis with the Gaussian spin model is consistent with the identification of spin-orbit misalignment using the Default spin model; with µ p = σ p = 0 disfavored. For the Default spin model, we verified that the signature of spin-orbit misalignment was not a spurious prior artifact, finding that our posterior lies safely outside the artificial exclusion region in Fig. 9. A similar exclusion region also exists for the Gaussain model around µ p = σ p = 0 , but our estimate of its exact size is subject to sampling uncertainties driven by the relatively small number of prior samples close to χ p = 0 for each event.  \nWith the Gaussian model, we also find evidence that at least some BHs have anti-aligned spins, with θ > 90 · , such that χ eff < 0 . To further evaluate the robustness of our Gaussian model fits, in Fig. 26 we show posterior predictive comparisons between predicted and empirical catalogs of χ eff and χ p measurements. The light blue bands mark 90% credible bounds on the predicted cumulative distribution of observed effective inspiral spin parameter values, given our posterior on the Gaussian model parameters. The dark shaded regions, meanwhile, show 90% credible bounds on the true distribution observed within GWTC-2, achieved by reweighting single event χ eff and χ p by repeated random draws from the Gaussian hyper-posterior. A similar predictive comparison between observation and an alternative strictly positive model, in which the effective inspiral spin distribution is truncated on the interval 0 ≤ χ eff ≤ 1 , reveals possible tension; when asserting that all effective inspiral spins are positive, the resulting population model underpredicts the number of observations with χ eff < 0 . 1 approximately 75% of the time.  \n<!-- image -->  \nFigure 26. Population predictive checks for the effective inspiral spin parameter χ eff ( left ) and the effective precession spin parameter χ p ( right ) of BBH mergers using the Gaussian spin model. The light shaded regions show the central 90% credible bounds on the posterior predictive distributions. According to the model, we expect the observed distributions on χ eff and χ p to lie within the light shaded region 90% of the time. The dark shaded regions show the 90% credible bounds on the observed distributions in GWTC-2, found using the population-informed posteriors of the confident BBH events in GWTC-2. The overlap between the dark and light regions shows that the model passes the posterior predictive check. The results for the Default model are similar, indicating that both models are a good fit to the data.  \n<!-- image -->  \nAs described in the main text, we can leverage the assumption that BBH with negative values of χ eff are formed dynamically to infer the fraction of binaries formed via dynamical ( f d ) and isolated ( f i ) channels; see Eq. (7). In Fig. 28 we show the corresponding posterior distributions for f i and f d .  \nThere are several sources of possible bias that might influence our Gaussian model conclusions. One possible source of bias is the mass model presumed for the Gaussian spin analysis. As noted above, measurements of a binary's χ eff and mass ratio q are generally anti-correlated (Ng et al. 2018a). Therefore, our particular choice of β q = 1 . 3 could conceivably affect conclusions regarding the χ eff distribution. We have directly verified that the results in Fig. 27 remain robust under different fiducial choices of β q between -1 . 5 and 2 .  \nAnother source of bias may be the Gaussian functional form we impose on the χ eff and χ p distributions, enforcing a unimodal distribution with smooth tails. As discussed in Sec. 5.2, though, the Default spin model yields near-identical χ eff and χ p distributions, despite its different parametrization and different physical assumptions. As an additional check, we repeat the Gaussian spin analysis on our data, truncating the χ eff distribution not on ( -1 , 1) , but on ( χ min eff , 1) , where χ min eff is inferred from the data. Figure 27 shows the marginal posterior for χ min eff . We find that χ min eff is constrained to be negative at 99% credibility, confirming that support for negative effective inspiral spin parameter is a feature of the data and not simply an artifact of the Gaussian model. In contrast, when we repeat the measurement of χ min eff using simulated catalogs drawn from ( i ) a Gaussian population truncated to strictly positive values, and ( ii ) a pair of delta functions at χ eff = 0 and 0 . 1 , we correctly observe consistency with χ min eff = 0 . Examining χ min eff is therefore a useful safeguard even when the true population is poorly fit by a Gaussian. However, if the χ eff distribution deviates too strongly from a Gaussian functional form, then it may remain possible to spuriously conclude that χ min eff < 0 . In the case of significant tidal torques, for example, it is predicted that the χ eff distribution is strongly bimodal, with effective inspiral spins clustered near χ eff = 0 or 1 , with no support at χ eff < 0 (Kushnir et al. 2016; Zaldarriaga et al. 2018; Bavera et al. 2019; Belczynski et al. 2020). To illustrate how results can be biased from model misspecification, we analyze mock observations drawn from a similar bimodal distribution. Using this intentionally misspecified model, we incorrectly conclude that χ min eff < 0 . However, in cases of such extreme model mismatch, we expect that our data would fail the predictive check of Fig. 26. Indeed, Fig. 29 shows the result of such a predictive check in the case of the bimodal, tidally-torqued χ eff distribution. When fitting this population with a Gaussian model, we find that the model overpredicts the fraction of mock observations with negative χ eff as well as the range 0 . 3 glyph[lessorsimilar] χ eff glyph[lessorsimilar] 0 . 8 , showing a clear deviation in the predicted and observed cumulative χ eff distributions.  \n<!-- image -->  \neff  \nFigure 27. Posterior distribution for χ min eff , below which we truncate the Gaussian χ eff distribution. While the results shown in the main text presume χ min eff = -1 , in Sec. D.2 we elevate χ min eff to a free hyper-parameter to be determined by the data; the resulting marginalized posterior distribution is shown here. In this case, we exclude χ min eff ≥ 0 at 99% credibility. This finding affirms that the signatures of anti-aligned BH spins are present in our BBH catalog, and not a bias due to our choice of parameterized spin model.  \n<!-- image -->  \nFigure 28. Posterior distributions for f i ( left ) and f d ( right ) of BBH mergers as defined Eq. 7. Assuming all binaries with χ eff < 0 are dynamically formed, f d corresponds to the fraction of dynamically assembled BBH systems. f i corresponds to the fraction of BBH systems formed through the isolated channel.  \n<!-- image -->", 'D.3. Multi Spin model': 'This model is an extension of the Truncated mass model with an additional Gaussian component. It is similar to the Power Law + Peak model, but there are several differences. First, the high-mass subpopulation in Multi Spin is described by a Gaussian in both m 1 and m 2 (up to the m 1 ≥ m 2 truncation) while Power Law + Peak only models m 1 as a Gaussian, and assumes that all BBH systems are described by a power-law distribution in mass ratio q . Most importantly, as its name suggests, Multi Spin allows each subpopulation to have its own independent spin distribution, each of which follows the Default model, with ζ = 1 . This allows us to probe whether the spin distribution varies with mass. The parameters for Multi Spin are summarized in Table 11.  \nFigure 29. Example of a failed posterior predictive check for the effective inspiral spin parameter χ eff distribution. We fit the Gaussian spin model to a mock catalog drawn from the strongly bimodal χ eff distributions predicted in the presence of tidal torques (Kushnir et al. 2016; Zaldarriaga et al. 2018; Bavera et al. 2019; Belczynski et al. 2020). The dark shaded region shows the central 90% credible bounds on the mock observed cumulative χ eff distribution, while the dark region corresponds to that predicted by the model. In this case, our Gaussian model is an extremely poor representation of the underlying χ eff distribution, and so significant tension is seen between the observed and predicted distributions, with the model overpredicting the fraction of observations with χ eff < 0 as well the fraction of observations with 0 . 3 glyph[lessorsimilar] χ eff glyph[lessorsimilar] 0 . 8 .  \n<!-- image -->  \nFigure 30. A posterior predictive check: the cumulative density function (CDF) for the Power Law Evolution model. The model is a good fit to the data. GW190814 is excluded from this analysis.  \n<!-- image -->', 'E. REDSHIFT EVOLUTION MODELS': 'The power-law redshift evolution model parameterizes the merger rate density as  \nR ( z ) = R 0 (1 + z ) κ , (E16)  \nwhere R 0 denotes the merger rate density at z = 0 . This implies that the redshift distribution is  \ndN dz = C dV c dz (1 + z ) κ -1 , (E17)  \nwhere dV c /dz is the differential comoving volume, and C is related to R 0 by  \nR 0 = C dV c dz [∫ z max 0 dV c dz (1 + z ) κ -1 ] -1 . (E18)  \nTable 11. Summary of Multi Spin parameters.  \nWe take z max = 2 . 3 in the analysis, as this is a conservative upper bound on the redshift at which we could detect BBH systems during O3a within the mass range considered here. When fitting this model, we employ a uniform prior on κ centered at κ = 0 . The value κ = 0 corresponds to no evolution; i.e., a merger rate that is uniform-in-comoving volume and source-frame time. We take a sufficiently wide prior so that the likelihood is entirely within the prior range, κ ∈ ( -6 , 6) .', 'F. GRAVITATIONAL-WAVE LENSING': 'It has been suggested that gravitational-wave lensing could bias the estimate of binary masses (Dai et al. 2017; Ng et al. 2018b; Li et al. 2018; Oguri 2018; Hannuksela et al. 2019), which could lead to a biased population inference. However, based on the predictions on the number of expected gravitational-wave sources and the distribution of galaxy lenses in the Universe, Li et al. (2018); Oguri (2018) predict that only around one in a thousand observed events are lensed, although this estimate can vary depending on the redshift evolution of the merger-rate density. The lensing rate  \nis expected to be similarly rare for galaxy clusters (Smith et al. 2018). As the expected lensing rate is low compared to the number of events in GWTC-2, we assume that all events in our sample are unlensed.'}