PT Journal AU ATLAS Collaboration (Aad, Gea Amos, KR Aparisi Pozo, JA Bailey, AJ Bouchhar, N Cabrera Urban, S Cantero, J Cardillo, F Castillo Gimenez, V Chitishvili, M Costa, MJ Didenko Escobar, C Fiorini, L Fullana Torregrosa, E Fuster, J Garcia, C Garcia Navarro, JE Gomez Delegido, AJ Gonzalez de la Hoz, S Gonzalvo Rodriguez, GR Guerrero Rojas, JGR Higon-Rodriguez, E Lacasta, C Lozano Bahilo, JJ Marti-Garcia, S Martinez Agullo, P Miralles Lopez, M Mitsou, VA Monsonis Romero, L Moreno Llacer, M Munoz Perez, D Navarro-Gonzalez, J Poveda, J Prades Ibañez, A Rubio Jimenez, A Ruiz-Martinez, A Sabatini, P Salt, J Sanchez Sebastian, V Sayago Galvan, I Senthilkumar, V Soldevila, U Sanchez, J Torro Pastor, E Valero, A Valls Ferrer, JA Varriale, L Villaplana Perez, M Vos, M TI Search for Higgs boson pair production in association with a vector boson in pp collisions at √s=13 TeV with the ATLAS detector SO European Physical Journal C JI Eur. Phys. J. C PY 2023 BP 519 - 41pp VL 83 IS 6 DI 10.1140/epjc/s10052-023-11559-y LA English AB This paper reports a search for Higgs boson pair (hh) production in association with a vector boson ( W or Z) using 139 fb(-1) of proton-proton collision data at root s = 13 TeV recorded with the ATLAS detector at the Large Hadron Collider. The search is performed in final states in which the vector boson decays leptonically (W -> l nu, Z -> ll, nu nu with l = e, mu) and the Higgs bosons each decay into a pair of b-quarks. It targets Vhh signals from both non-resonant hh production, present in the Standard Model (SM), and resonant hh production, as predicted in some SM extensions. A 95% confidence-level upper limit of 183 (87) times the SM cross-section is observed (expected) for non-resonant Vhh production when assuming the kinematics are as expected in the SM. Constraints are also placed on Higgs boson coupling modifiers. For the resonant search, upper limits on the production cross-sections are derived for two specific models: one is the production of a vector boson along with a neutral heavy scalar resonance H, in the mass range 260-1000GeV, that decays into hh, and the other is the production of a heavier neutral pseudoscalar resonance A that decays into a Z boson and H boson, where the A boson mass is 360-800GeV and the H boson mass is 260-400GeV. Constraints are also derived in the parameter space of two-Higgs-doublet models. ER