@article{WitzorkyParamonovBouaklineetal.2021, author = {Witzorky, Christoph and Paramonov, Guennaddi and Bouakline, Foudhil and Jaquet, Ralph and Saalfrank, Peter and Klamroth, Tillmann}, title = {Gaussian-type orbital calculations for high harmonic generation in vibrating molecules}, series = {Journal of chemical theory and computation}, volume = {17}, journal = {Journal of chemical theory and computation}, number = {12}, publisher = {American Chemical Society}, address = {Washington}, issn = {1549-9618}, doi = {10.1021/acs.jctc.1c00837}, pages = {7353 -- 7365}, year = {2021}, abstract = {The response of the hydrogen molecular ion, H-2(+), to few-cycle laser pulses of different intensities is simulated. To treat the coupled electron-nuclear motion, we use adiabatic potentials computed with Gaussian-type basis sets together with a heuristic ionization model for the electron and a grid representation for the nuclei. Using this mixed-basis approach, the time-dependent Schrodinger equation is solved, either within the Born-Oppenheimer approximation or with nonadiabatic couplings included. The dipole response spectra are compared to all-grid-based solutions for the three-body problem, which we take as a reference to benchmark the Gaussian-type basis set approaches. Also, calculations employing the fixed-nuclei approximation are performed, to quantify effects due to nuclear motion. For low intensities and small ionization probabilities, we get excellent agreement of the dynamics using Gaussian-type basis sets with the all-grid solutions. Our investigations suggest that high harmonic generation (HHG) and high-frequency response, in general, can be reliably modeled using Gaussian-type basis sets for the electrons for not too high harmonics. Further, nuclear motion destroys electronic coherences in the response spectra even on the time scale of about 30 fs and affects HHG intensities, which reflect the electron dynamics occurring on the attosecond time scale. For the present system, non-Born-Oppenheimer effects are small. The Gaussian-based, nonadiabatically coupled, time-dependent multisurface approach to treat quantum electron-nuclear motion beyond the non-Born-Oppenheimer approximation can be easily extended to approximate wavefunction methods, such as time-dependent configuration interaction singles (TD-CIS), for systems where no benchmarks are available.}, language = {en} } @article{BouaklineTremblay2020, author = {Bouakline, Foudhil and Tremblay, Jean Christophe}, title = {Is it really possible to control aromaticity of benzene with light?}, series = {Physical chemistry, chemical physics : PCCP}, volume = {22}, journal = {Physical chemistry, chemical physics : PCCP}, number = {27}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c9cp06794a}, pages = {15401 -- 15412}, year = {2020}, abstract = {Recent theoretical investigations claim that tailored laser pulses may selectively steer benzene's aromatic ground state to localized non-aromatic excited states. For instance, it has been shown that electronic wavepackets, involving the two lowest electronic eigenstates, exhibit subfemtosecond charge oscillation between equivalent Kekule resonance structures. In this contribution, we show that such dynamical electron-localization in the molecule-fixed frame contravenes the principle of the indistinguishability of identical particles. This breach stems from a total omission of the nuclear degrees of freedom, giving rise to nonsymmetric electronic wavepackets under nuclear permutations. Enforcement of the latter leads to entanglement between the electronic and nuclear states. To obey quantum statistics, the entangled molecular states should involve compensating nuclear-permutation symmetries. This in turn engenders complete quenching of dynamical electron-localization in the molecule-fixed frame. Indeed, for the (six-fold) equilibrium geometry of benzene, group-theoretic analysis reveals that any electronic wavepacket exhibits a (D-6h) totally symmetric electronic density, at all times. Thus, our results clearly show that the six carbon atoms, and the six C-C bonds, always have equal Mulliken charges, and equal bond orders, respectively. However, electronic wavepackets may display dynamical localization of the electronic density in the space-fixed frame, whenever they involve both even and odd space-inversion (parity) or permutation-inversion symmetry. Dynamical spatial-localization can be probed experimentally in the laboratory frame, but it should not be deemed equivalent to charge oscillation between benzene's identical electronic substructures, such as Kekule resonance structures.}, language = {en} } @article{Bouakline2018, author = {Bouakline, Foudhil}, title = {Unambiguous signature of the berry phase in intense laser dissociation of diatomic molecules}, series = {The journal of physical chemistry letters}, volume = {9}, journal = {The journal of physical chemistry letters}, number = {9}, publisher = {American Chemical Society}, address = {Washington}, issn = {1948-7185}, doi = {10.1021/acs.jpclett.8b00607}, pages = {2271 -- 2277}, year = {2018}, abstract = {We report strong evidence of Berry phase effects in intense laser dissociation of D-2(+) molecules, manifested as Aharonov-Bohm-like oscillations in the photofragment angular distribution (PAD). Our calculations show that this interference pattern strongly depends on the parity of the diatom initial rotational state, (-1)(j). Indeed, the PAD local maxima (minima) observed in one case (j odd) correspond to local minima (maxima) in the other case (j even). Using simple topological arguments, we clearly show that such interference conversion is a direct signature of the Berry phase. The sole effect of the latter on the rovibrational wave function is a sign change of the relative phase between two interfering components, which wind in opposite senses around a light-induced conical intersection (LICI). Therefore, encirclement of the LICI leads to constructive (j odd) or destructive (j even) self-interference of the initial nuclear wavepacket in the dissociative limit. To corroborate our theoretical findings, we suggest an experiment of strong-field indirect dissociation of D-2(+) molecules, comparing the PAD of the ortho and para molecular species in directions nearly perpendicular to the laser polarization axis.}, language = {en} } @phdthesis{Bouakline2023, author = {Bouakline, Foudhil}, title = {Manifestations of Quantum-Mechanical Effects in Molecular Reaction Dynamics}, school = {Universit{\"a}t Potsdam}, pages = {316}, year = {2023}, abstract = {This habilitation thesis summarises the research work performed by the author during the last quindecennial period. The dissertation reflects his main research interests, which revolve around quantum dynamics of small-sized molecular systems, including their interactions with electromagnetic radiation or dissipative environments. This covers various dynamical processes that involve bound-bound, bound-free, and free-free molecular transitions. The latter encompass light-triggered rovibrational or rovibronic dynamics in bound molecules, molecular photodissociation induced by weak or strong laser fields, state-to-state reactive and/or inelastic molecular collisions, and phonon-driven vibrational relaxation of adsorbates at solid surfaces. Although the dissertation covers different topics of molecular reaction dynamics, most of these studies focus on nuclear quantum effects and their manifestations in experimental measures. The latter are assessed through comparison between quantum and classical predictions, and/or direct confrontation of theory and experiment. Most well known quantum concepts and effects will be encountered in this work. Yet, almost all these quantum notions find their roots in the central pillar of quantum theory, namely, the quantum superposition principle. Indeed, quantum coherence is the main source of most quantum effects, including interference, entanglement, and even tunneling. Thus, the common and predominant theme of all the investigations of this thesis is quantum coherence, and the survival or quenching of subsequent interference effects in various molecular processes. The lion's share of the dissertation is devoted to two associated quantum concepts, which are usually overlooked in computational molecular dynamics, viz. the Berry phase and identical nuclei symmetry. The importance of the latter in dynamical molecular processes and their direct fingerprints in experimental observables also rely very much on quantum coherence and entanglement. All these quantum phenomena are thoroughly discussed within the four main topics that form the core of this thesis. Each topic is described in a separate chapter, where it is briefly summarised and then illustrated with three peer-reviewed publications. The first topic deals with the relevance of quantum coherence/interference in molecular collisions, with a focus on the hydrogen-exchange reaction, H+H2 --> H2+H, and its isotopologues. For these collision processes, the significance of interference of probability amplitudes arises because of the existence of two main scattering pathways. The latter could be inelastic and reactive scattering, direct and time-delayed scattering, or two encircling reaction paths that loop in opposite senses around a conical intersection (CI) of the H3 molecular system. Our joint theoretical-experimental investigations of these processes reveal strong interference and geometric phase (GP) effects in state-to-state reaction probabilities and differential cross sections. However, these coherent effects completely cancel in integral cross sections and reaction rate constants, due to efficient dephasing of interference between the different scattering amplitudes. As byproducts of these studies, we highlight the discovery of two novel scattering mechanisms, which contradict conventional textbook pictures of molecular reaction dynamics. The second topic concerns the effect of the Berry phase on molecular photodynamics at conical intersections. To understand this effect, we developed a topological approach that separates the total molecular wavefunction of an unbound molecular system into two components, which wind in opposite senses around the conical intersection. This separation reveals that the only effect of the geometric phase is to change the sign of the relative phase of these two components. This in turn leads to a shift in the interference pattern of the molecular system---a phase shift that is reminiscient of the celebrated Aharonov-Bohm effect. This procedure is numerically illustrated with photodynamics at model standard CIs, as well as strong-field dissociation of diatomics at light-induced conical intersections (LICIs). Besides the fundamental aspect of these studies, their findings allow to interpret and predict the effect of the GP on the state-resolved or angle-resolved spectra of pump-probe experimental schemes, particularly the distributions of photofragments in molecular photodissociation experiments. The third topic pertains to the role of the indistinguishability of identical nuclei in molecular reaction dynamics, with an emphasis on dynamical localization in highly symmetric molecules. The main object of these studies is whether nuclear-spin statistics allow dynamical localization of the electronic, vibrational, or even rotational density on a specific molecular substructure or configuration rather than on another one which is identical (indistinguishable). Group-theoretic analysis of the symmetrized molecular wavefunctions of these systems shows that nuclear permutation symmetry engenders quantum entanglement between the eigenstates of the different molecular degrees of freedom. This subsequently leads to complete quenching of dynamical localization over indistinguishable molecular substructures---an observation that is in sharp contradiction with well known textbook views of iconic molecular processes. This is illustrated with various examples of quantum dynamics in symmetric double-well achiral molecules, such as the prototypical umbrella inversion motion of ammonia, electronic Kekul{\´e} dynamics in the benzene molecule, and coupled electron-nuclear dynamics in laser-induced indirect photodissociation of the dihydrogen molecular cation. The last part of the thesis is devoted to the development of approximate wavefunction approaches for phonon-induced vibrational relaxation of adsorbates (system) at surfaces (bath). Due to the so-called 'curse of dimensionality', these system-bath complexes cannot be handled with standard wavefunction methods. To alleviate the exponential scaling of the latter, we developed approximate yet quite accurate numerical schemes that have a polynomial scaling with respect to the bath dimensionality. The corresponding algorithms combine symmetry-based reductions of the full vibrational Hilbert space and iterative Krylov techniques. These approximate wavefunction approaches resemble the 'Bixon-Jortner model' and the more general 'quantum tier model'. This is illustrated with the decay of H-Si (D-Si) vibrations on a fully H(D)-covered silicon surface, which is modelled with a phonon-bath of more than two thousand oscillators. These approximate methods allow reliable estimation of the adsorbate vibrational lifetimes, and provide some insight into vibration-phonon couplings at solid surfaces. Although this topic is mainly computational, the developed wavefunction approaches permit to describe quantum entanglement between the system and bath states, and to embody some coherent effects in the time-evolution of the (sub-)system, which cannot be accounted for with the widely used 'reduced density matrix formalism'.}, language = {en} } @article{BouaklineSaalfrank2021, author = {Bouakline, Foudhil and Saalfrank, Peter}, title = {Seemingly asymmetric atom-localized electronic densities following laser-dissociation of homonuclear diatomics}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistry}, volume = {154}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistry}, number = {23}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/5.0049710}, pages = {10}, year = {2021}, abstract = {Recent experiments on laser-dissociation of aligned homonuclear diatomic molecules show an asymmetric forward-backward (spatial) electron-localization along the laser polarization axis. Most theoretical models attribute this asymmetry to interference effects between gerade and ungerade vibronic states. Presumably due to alignment, these models neglect molecular rotations and hence infer an asymmetric (post-dissociation) charge distribution over the two identical nuclei. In this paper, we question the equivalence that is made between spatial electron-localization, observed in experiments, and atomic electron-localization, alluded by these theoretical models. We show that (seeming) agreement between these models and experiments is due to an unfortunate omission of nuclear permutation symmetry, i.e., quantum statistics. Enforcement of the latter requires mandatory inclusion of the molecular rotational degree of freedom, even for perfectly aligned molecules. Unlike previous interpretations, we ascribe spatial electron-localization to the laser creation of a rovibronic wavepacket that involves field-free molecular eigenstates with opposite space-inversion symmetry i.e., even and odd parity. Space-inversion symmetry breaking would then lead to an asymmetric distribution of the (space-fixed) electronic density over the forward and backward hemisphere. However, owing to the simultaneous coexistence of two indistinguishable molecular orientational isomers, our analytical and computational results show that the post-dissociation electronic density along a specified space-fixed axis is equally shared between the two identical nuclei-a result that is in perfect accordance with the principle of the indistinguishability of identical particles. Published under an exclusive license by AIP Publishing.}, language = {en} }