Coupled diabatic potential energy surfaces for studying the nonadiabatic dynamics at conical intersections in angular resolved photodetachment simulations of OHF- --> OHF+e-

16 pages, 14 figures, 1 table.-- PACS nrs.: 82.20.Kh; 82.30.Cf; 33.80.Eh; 82.20.Db; 82.20.Hf.An energy-based method is proposed for the diabatization of the OH(2Π)+F(2P) --> O(3P)+HF(1Σ+) reaction. It is demonstrated that the diabatic representation obtained is regularized, i.e., the residual derivative couplings do not present singularities at the conical intersections appearing along the reaction path. This method only requires the knowledge of the 1,2 3A'' and 1 3A' eigenvalues and does not require any adjustable parameter. Thus, many convergence problems arising in other derivative-based diabatization methods are avoided, and the description of the configuration space along the reaction path is enormously simplified. Three-dimensional coupled diabatic energy surfaces are obtained by an interpolation procedure using ≈ 4000 accurate ab initio points. The angular resolved photodetachment cross sections are obtained in the diabatic and adiabatic representations using a wave packet method. An excellent agreement is obtained with recent experimental data [D. M. Neumark, Phys. Chem. Chem. Phys. 7, 433 (2005)] for high electron kinetic energies where only the triplet electronic states contribute.This work has been supported by DGICYT (Ministerio de Educación y Ciencia, Spain) under Grant Nos. CTQ2004-02415/BQU and CTQ2004-06615/BQU.Peer reviewe