Alberto Castro, ChemPhysChem 17, 1601 (2016)

The combination of nonadiabatic Ehrenfest-path molecular dynamics (EMD) based on time-dependent density functional theory (TDDFT) and quantum optimal control formalism (QOCT) was used to optimize the shape of ultra-short laser pulses to achieve photodissociation of a hydrogen molecule and the trihydrogen cation H3 + . This work completes a previous one [A. Castro, ChemPhysChem, 2013, 14, 1488–1495], in which the same objective was achieved by demonstrating the combination of QOCT and TDDFT for many-electron systems on static nuclear potentials. The optimization model, therefore, did not include the nuclear movement and the obtained dissociation mechanism could only be sequential: fast laser-assisted electronic excitation to nonbonding states (during which the nuclei are considered to be static), followed by field-free dissociation. Here, in contrast, the optimization was performed with the QOCT constructed on top of the full dynamic model comprised of both electrons and nuclei, as described within EMD based on TDDFT. This is the first numerical demonstration of an optimal control formalism for a hybrid quantum–classical model, that is, a molecular dynamics method.