@article{Tremblay2011, author = {Tremblay, Jean Christophe}, title = {Laser control of molecular excitations in stochastic dissipative media}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {134}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {17}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.3587093}, pages = {16}, year = {2011}, abstract = {In the present work, ideas for controlling photochemical reactions in dissipative environments using shaped laser pulses are presented. New time-local control algorithms for the stochastic Schrodinger equation are introduced and compared to their reduced density matrix analog. The numerical schemes rely on time-dependent targets for guiding the reaction along a preferred path. The methods are tested on the vibrational control of adsorbates at metallic surfaces and on the ultrafast electron dynamics in a strong dissipative medium. The selective excitation of the specific states is achieved with improved yield when using the new algorithms. Both methods exhibit similar convergence behavior and results compare well with those obtained using local optimal control for the reduced density matrix. The favorable scaling of the methods allows to tackle larger systems and to control photochemical reactions in dissipative media of molecules with many more degrees of freedom.}, language = {en} } @misc{ToepferTremblay2016, author = {T{\"o}pfer, Kai and Tremblay, Jean Christophe}, title = {How surface reparation prevents catalytic oxidation of carbon monoxide on atomic gold at defective magnesium oxide surfaces}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-394978}, pages = {8}, year = {2016}, abstract = {In this contribution, we study using first principles the co-adsorption and catalytic behaviors of CO and O2 on a single gold atom deposited at defective magnesium oxide surfaces. Using cluster models and point charge embedding within a density functional theory framework, we simulate the CO oxidation reaction for Au1 on differently charged oxygen vacancies of MgO(001) to rationalize its experimentally observed lack of catalytic activity. Our results show that: (1) co-adsorption is weakly supported at F0 and F2+ defects but not at F1+ sites, (2) electron redistribution from the F0 vacancy via the Au1 cluster to the adsorbed molecular oxygen weakens the O2 bond, as required for a sustainable catalytic cycle, (3) a metastable carbonate intermediate can form on defects of the F0 type, (4) only a small activation barrier exists for the highly favorable dissociation of CO2 from F0, and (5) the moderate adsorption energy of the gold atom on the F0 defect cannot prevent insertion of molecular oxygen inside the defect. Due to the lack of protection of the color centers, the surface becomes invariably repaired by the surrounding oxygen and the catalytic cycle is irreversibly broken in the first oxidation step.}, language = {en} } @article{ScholzLindnerLoncaricetal.2019, author = {Scholz, Robert and Lindner, Steven and Loncaric, Ivor and Tremblay, Jean Christophe and Juaristi, J. and Alducin, Maite and Saalfrank, Peter}, title = {Vibrational response and motion of carbon monoxide on Cu(100) driven by femtosecond laser pulses: Molecular dynamics with electronic friction}, series = {Physical review : B, Condensed matter and materials physics}, volume = {100}, journal = {Physical review : B, Condensed matter and materials physics}, number = {24}, publisher = {American Physical Society}, address = {College Park}, issn = {2469-9950}, doi = {10.1103/PhysRevB.100.245431}, pages = {20}, year = {2019}, abstract = {Carbon monoxide on copper surfaces continues to be a fascinating, rich microlab for many questions evolving in surface science. Recently, hot-electron mediated, femtosecond-laser pulse induced dynamics of CO molecules on Cu(100) were the focus of experiments [Inoue et al., Phys. Rev. Lett. 117, 186101 (2016)] and theory [Novko et al., Phys. Rev. Lett. 122, 016806 (2019)], unraveling details of the vibrational nonequilibrium dynamics on ultrashort (subpicoseconds) timescales. In the present work, full-dimensional time-resolved hot-electron driven dynamics are studied by molecular dynamics with electronic friction (MDEF). Dissipation is included by a friction term in a Langevin equation which describes the coupling of molecular degrees of freedom to electron-hole pairs in the copper surface, calculated from gradient-corrected density functional theory (DFT) via a local density friction approximation (LDFA). Relaxation due to surface phonons is included by a generalized Langevin oscillator model. The hot-electron induced excitation is described via a time-dependent electronic temperature, the latter derived from an improved two-temperature model. Our parameter-free simulations on a precomputed potential energy surface allow for excellent statistics, and the observed trends are confirmed by on-the-fly ab initio molecular dynamics with electronic friction (AIMDEF) calculations. By computing time-resolved frequency maps for selected molecular vibrations, instantaneous frequencies, probability distributions, and correlation functions, we gain microscopic insight into hot-electron driven dynamics and we can relate the time evolution of vibrational internal CO stretch-mode frequencies to measured data, notably an observed redshift. Quantitatively, the latter is found to be larger in MDEF than in experiment and possible reasons are discussed for this observation. In our model, in addition we observe the excitation and time evolution of large-amplitude low-frequency modes, lateral CO surface diffusion, and molecular desorption. Effects of surface atom motion and of the laser fluence are also discussed.}, language = {en} } @misc{TremblayBlancoRey2015, author = {Tremblay, Jean Christophe and Blanco-Rey, Maria}, title = {Manipulating interfacial hydrogens at palladium via STM}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-99511}, pages = {11}, year = {2015}, abstract = {In this contribution, we provide a detailed dynamical analysis of the interfacial hydrogen migration mediated by scanning tunneling microscopy (STM). Contributions from the STM-current and from the non-adiabatic couplings are taken into account using only first principle models. The slight asymmetry of the tunnelling rates with respect to the potential bias sign inferred from experimental observations is reproduced by weighting the contributions of the metal acceptor-donor states for the propagation of the impinging electrons. The quasi-thermal inelastic collision mechanism is treated perturbatively. The influence of hydrogen pre-coverage is also investigated using new potential energy surfaces obtained from periodic density functional theory calculations. Fully quantum dynamical simulations of the system evolution are performed by solving the Pauli master equation, providing insight into the reaction mechanism of STM manipulation of subsurface hydrogens. It is observed that the hydrogen impurity favors resurfacing over occupation of the bulk and subsurface sites whenever possible. The present simulations give strong indication that the experimentally observed protuberances after STM-excitation are due to hydrogen accumulating in the vicinity of the surface.}, language = {en} } @article{TremblaySaalfrank2009, author = {Tremblay, Jean Christophe and Saalfrank, Peter}, title = {Selective subsurface absorption of hydrogen in palladium using laser distillation}, issn = {0021-9606}, doi = {10.1063/1.3212695}, year = {2009}, abstract = {A theoretical model for the selective subsurface absorption of atomic hydrogen in a Pd(111) surface by infrared (IR) laser pulses is presented. The dynamics of the adsorbate is studied within the reduced density matrix approach. Energy and phase relaxation of the hydrogen atom are treated using the semigroup formalism. The vibrational excitation leading to subsurface absorption is performed using rationally designed pulses as well as IR laser pulses optimized on- the-fly. It is shown that dissipation can be used as a tool to transfer population to an otherwise inaccessible state via a mechanism known as "laser distillation." We demonstrate that when the reaction path is generalized from a reduced one-dimensional to full three-dimensional treatment of the system, the laser control strategy can prove very different.}, language = {en} } @article{TremblayMonturetSaalfrank2011, author = {Tremblay, Jean Christophe and Monturet, Serge and Saalfrank, Peter}, title = {The Effects of electron-hole pair coupling on the infrared laser-controlled vibrational excitation of NO on Au(111)}, series = {The journal of physical chemistry : A, Molecules, spectroscopy, kinetics, environment \& general theory}, volume = {115}, journal = {The journal of physical chemistry : A, Molecules, spectroscopy, kinetics, environment \& general theory}, number = {39}, publisher = {American Chemical Society}, address = {Washington}, issn = {1089-5639}, doi = {10.1021/jp205902k}, pages = {10698 -- 10707}, year = {2011}, abstract = {In this work, we present theoretical simulations of laser-driven vibrational control of NO adsorbed on a gold surface. Our goal is to tailor laser pulses to selectively excite specific modes and vibrational eigenstates, as well as to favor photodesorption of the adsorbed molecule. To this end, various control schemes and algorithms are applied. For adsorbates at metallic surfaces, the creation of electron hole pairs in the substrate is known to play a dominant role in the transfer of energy from the system to the surroundings. These nonadiabatic couplings are included perturbatively in our reduced density matrix simulations using a generalization of the state-resolved position-dependent anharmonic rate model we recently introduced. An extension of the reduced density matrix is also proposed to provide a sound model for photodesorption in dissipative systems.}, language = {en} } @article{TremblayKlinkuschKlamrothetal.2011, author = {Tremblay, Jean Christophe and Klinkusch, Stefan and Klamroth, Tillmann and Saalfrank, Peter}, title = {Dissipative many-electron dynamics of ionizing systems}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {134}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {4}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.3532410}, pages = {10}, year = {2011}, abstract = {In this paper, we perform many-electron dynamics using the time-dependent configuration-interaction method in its reduced density matrix formulation (rho-TDCI). Dissipation is treated implicitly using the Lindblad formalism. To include the effect of ionization on the state-resolved dynamics, we extend a recently introduced heuristic model for ionizing states to the rho-TDCI method, which leads to a reduced density matrix evolution that is not norm-preserving. We apply the new method to the laser-driven excitation of H(2) in a strongly dissipative environment, for which the state-resolve lifetimes are tuned to a few femtoseconds, typical for dynamics of adsorbate at metallic surfaces. Further testing is made on the laser-induced intramolecular charge transfer in a quinone derivative as a model for a molecular switch. A modified scheme to treat ionizing states is proposed to reduce the computational burden associated with the density matrix propagation, and it is thoroughly tested and compared to the results obtained with the former model. The new approach scales favorably (similar to N(2)) with the number of configurations N used to represent the reduced density matrix in the rho-TDCI method, as compared to a N(3) scaling for the model in its original form.}, language = {en} } @article{TremblayKrauseKlamrothetal.2010, author = {Tremblay, Jean Christophe and Krause, Pascal and Klamroth, Tillmann and Saalfrank, Peter}, title = {Time-dependent response of dissipative electron systems}, issn = {1050-2947}, doi = {10.1103/Physreva.81.063420}, year = {2010}, abstract = {We present a systematic study of the influence of energy and phase relaxation on dynamic polarizability simulations in the linear response regime. The nonperturbative approach is based on explicit electron dynamics using short laser pulses of low intensities. To include environmental effects on the property calculation, we use the time- dependent configuration-interaction method in its reduced density matrix formulation. Both energy dissipation and nonlocal pure dephasing are included. The explicit treatment of time-resolved electron dynamics gives access to the phase shift between the electric field and the induced dipole moment, which can be used to define a useful uncertainty measure for the dynamic polarizability. The nonperturbative treatment is compared to perturbation theory expressions, as applied to a simple model system, the rigid H-2 molecule. It is shown that both approaches are equivalent for low field intensities, but the time-dependent treatment provides complementary information on the phase of the induced dipole moment, which allows for the definition of an uncertainty associated with the computation of the dynamic polarizability in the linear response regime.}, language = {en} } @article{FuechselTremblayKlamrothetal.2012, author = {F{\"u}chsel, Gernot and Tremblay, Jean Christophe and Klamroth, Tillmann and Saalfrank, Peter}, title = {Selective excitation of molecule-surface vibrations in H2 and D2 dissociatively adsorbed on Ru(0001)}, series = {Israel journal of chemistry}, volume = {52}, journal = {Israel journal of chemistry}, number = {5}, publisher = {Wiley-VCH}, address = {Weinheim}, issn = {0021-2148}, doi = {10.1002/ijch.201100097}, pages = {438 -- 451}, year = {2012}, abstract = {In this contribution we report about the selective vibrational excitation of H2 and D2 on Ru(0001) as an example for nonadiabatic coupling of an open quantum system to a dissipative environment. We investigate the possibility of achieving state-selective vibrational excitations of H2 and D2 adsorbed on a Ru(0001) surface using picosecond infrared laser pulses. The systems behavior is explored using pulses that are rationally designed and others that are optimized using a time-local variant of Optimal Control Theory. The effects of dissipation on the laser-driven dynamics are studied using the reduced-density matrix formalism. The non-adiabatic couplings between adsorbate and surface are computed perturbatively, for which our recently introduced state-resolved anharmonic rate model is used. It is shown that mode- and state-selective excitation can be achieved in the absence of dissipation when using optimized laser pulses. The inclusion of dissipation in the model reduces the state selectivity and the population transfer yield to highly excited states. In this case, mode activation is most effectively realized by a rational pulse of carefully chosen duration rather than by a locally optimized pulse.}, language = {en} } @article{TremblayFuechselSaalfrank2012, author = {Tremblay, Jean Christophe and F{\"u}chsel, Gernot and Saalfrank, Peter}, title = {Excitation, relaxation, and quantum diffusion of CO on copper}, series = {Physical review : B, Condensed matter and materials physics}, volume = {86}, journal = {Physical review : B, Condensed matter and materials physics}, number = {4}, publisher = {American Physical Society}, address = {College Park}, issn = {1098-0121}, doi = {10.1103/PhysRevB.86.045438}, pages = {13}, year = {2012}, abstract = {We investigate the effect of intermode coupling and anharmonicity on the excitation and relaxation dynamics of CO on Cu(100). The nonadiabatic coupling of the adsorbate to the surface is treated perturbatively using a position-dependent state-resolved transition rate model. Using the potential energy surface of Marquardt et al. [J. Chem. Phys. 132, 074108 (2010)], which provides an accurate description of intermode interactions, we propose a four-dimensional model that represents simultaneously the diffusion and the desorption of the adsorbate. The system is driven by both rational and optimized infrared laser pulses to favor either selective mode and state excitations or lateral displacement along the diffusion coordinate. The dissipative dynamics is simulated using the reduced density matrix in its Lindblad form. We show that coupling between the degrees of freedom, mediated by the creation and annihilation of electron-hole pairs in the metal substrate, significantly affects the system excitation and relaxation dynamics. In particular, the angular degrees of freedom appear to play an important role in the energy redistribution among the molecule-surface vibrations. We also show that coherent excitation using simple IR pulses can achieve population transfer to a specific target to some extent but does not allow enforcement of the directionality to the diffusion motion.}, language = {en} } @article{FuechselTremblayKlamrothetal.2012, author = {F{\"u}chsel, Gernot and Tremblay, Jean Christophe and Klamroth, Tillmann and Saalfrank, Peter and Frischkorn, C.}, title = {Concept of a single temperature for highly nonequilibrium laser-induced hydrogen desorption from a ruthenium surface}, series = {Physical review letters}, volume = {109}, journal = {Physical review letters}, number = {9}, publisher = {American Physical Society}, address = {College Park}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.109.098303}, pages = {5}, year = {2012}, abstract = {Laser-induced condensed phase reactions are often interpreted as nonequilibrium phenomena that go beyond conventional thermodynamics. Here, we show by Langevin dynamics and for the example of femtosecond-laser desorption of hydrogen from a ruthenium surface that light adsorbates thermalize rapidly due to ultrafast energy redistribution after laser excitation. Despite the complex reaction mechanism involving hot electrons in the surface region, all desorption product properties are characterized by equilibrium distributions associated with a single, unique temperature. This represents an example of ultrahot chemistry on the subpicosecond time scale.}, language = {en} } @article{FuechselKlamrothTremblayetal.2010, author = {F{\"u}chsel, Gernot and Klamroth, Tillmann and Tremblay, Jean Christophe and Saalfrank, Peter}, title = {Stochastic approach to laser-induced ultrafast dynamics : the desorption of H-2/D-2 from Ru(0001)}, issn = {1463-9076}, doi = {10.1039/C0cp00895h}, year = {2010}, abstract = {The desorption of molecular hydrogen and deuterium induced by femtosecond-laser pulses is studied theoretically for the so-called DIMET (Desorption Induced by Multiple Electronic Transitions) process. These investigations are based on nonadiabatic classical Monte Carlo trajectory (CMCT) simulations on a ground and an excited state potential energy surface, including up to all six adsorbate degrees of freedom. The focus is on the hot-electron mediated energy transfer from the surface to the molecule and back, and the energy partitioning between the different degrees of freedom of the desorbing molecules. We first validate for a two-mode model comprising the desorption mode and the internal vibrational coordinate, the classical Monte Carlo trajectory method by comparing with Monte Carlo wavepacket (MCWP) calculations arising from a fully quantum mechanical open-system density matrix treatment. We then proceed by extending the CMCT calculations to include all six nuclear degrees of freedom of the desorbing molecule. This allows for a detailed comparison between theory and experiment concerning isotope effects, energy partitioning (translational, vibrational, and rotational energies and their distributions), and the dependence of these properties on the laser fluence. The most important findings are as follows. (i) CMCT agrees qualitative with the MCWP scheme. (ii) The basic experimental features such as the large isotope effect, the non-linear increase of yield with laser fluence, translationally hot products (in the order of several 1000 K) and non-equipartitioning of translational and internal energies (E-trans > E- vib > E-rot) are well reproduced. (iii) Predictions concerning a strong angular dependence of translational energies at large observation angles are also made.}, language = {en} } @article{FuechselTremblayKlamrothetal.2013, author = {F{\"u}chsel, Gernot and Tremblay, Jean Christophe and Klamroth, Tillmann and Saalfrank, Peter}, title = {Quantum dynamical simulations of the femtosecond-laser-induced ultrafast desorption of H2 and D2 from Ru(0001)}, series = {ChemPhysChem : a European journal of chemical physics and physical chemistry}, volume = {14}, journal = {ChemPhysChem : a European journal of chemical physics and physical chemistry}, number = {7}, publisher = {Wiley-VCH}, address = {Weinheim}, issn = {1439-4235}, doi = {10.1002/cphc.201200940}, pages = {1471 -- 1478}, year = {2013}, abstract = {We investigate the recombinative desorption of hydrogen and deuterium from a Ru(0001) surface initiated by femtosecond laser pulses. We adopt a quantum mechanical two-state model including three molecular degrees of freedom to describe the dynamics within the desorption induced by electronic transition (DIET) limit. The energy distributions as well as the state-resolved and ensemble properties of the desorbed molecules are analyzed in detail by using the time-energy method. Our results shed light on the experimentally observed 1) large isotopic effects regarding desorption yields and translational energies and 2) the nonequal energy partitioning into internal and translational modes. In particular, it is shown that a single temperature is sufficient to characterize the energy distributions for all degrees of freedom. Further, we confirm that quantization effects play an important role in the determination of the energy partitioning.}, language = {en} } @article{FuechselTremblaySaalfrank2014, author = {F{\"u}chsel, Gernot and Tremblay, Jean Christophe and Saalfrank, Peter}, title = {A six-dimensional potential energy surface for Ru(0001)(2x2):CO}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {141}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {9}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.4894083}, pages = {9}, year = {2014}, abstract = {We present a new global ground state potential energy surface (PES) for carbon monoxide at a coverage of 1/4, on a rigid Ru(0001) surface [Ru(0001)(2x2):CO]. All six adsorbate degrees of freedom are considered. For constructing the PES, we make use of more than 90 000 points calculated with periodic density functional theory using the RPBE exchange-correlation functional and an empirical van der Waals correction. These points are used for interpolation, utilizing a symmetry-adapted corrugation reducing procedure (CRP). Three different interpolation schemes with increasing accuracy have been realized, giving rise to three flavours of the CRP PES. The CRP PES yields in agreement with the DFT reference and experiments, the atop position of CO to be the most stable adsorption geometry, for the most accurate interpolation with an adsorption energy of 1.69 eV. The CRP PES shows that diffusion parallel to the surface is hindered by a barrier of 430 meV, and that dissociation is facilitated but still activated. As a first "real" application and further test of the new potential, the six-dimensional vibrational Schrodinger equation is solved variationally to arrive at fully coupled, anharmonic frequencies and vibrational wavefunctions for the vibrating, adsorbed CO molecule. Good agreement with experiment is found also here. Being analytical, the new PES opens an efficient way towards multidimensional dynamics. (C) 2014 AIP Publishing LLC.}, language = {en} } @incollection{SaalfrankBedurkeHeideetal.2020, author = {Saalfrank, Peter and Bedurke, Florian and Heide, Chiara and Klamroth, Tillmann and Klinkusch, Stefan and Krause, Pascal and Nest, Mathias and Tremblay, Jean Christophe}, title = {Molecular attochemistry: correlated electron dynamics driven by light}, series = {Advances in quantum chemistry}, volume = {81}, booktitle = {Advances in quantum chemistry}, editor = {Ruud, Kenneth and Br{\"a}ndas, Erkki J.}, publisher = {Elsevier}, address = {Amsterdam [u.a.]}, isbn = {978-0-12-819757-8}, issn = {0065-3276}, doi = {10.1016/bs.aiq.2020.03.001}, pages = {15 -- 50}, year = {2020}, abstract = {Modern laser technology and ultrafast spectroscopies have pushed the timescales for detecting and manipulating dynamical processes in molecules from the picosecond over femtosecond domains, to the attosecond regime (1 as = 10(-18) s). This way, real-time dynamics of electrons after their photoexcitation can be probed and manipulated. In particular, experiments are moving more and more from atomic and solid state systems to molecules, opening the fields of "molecular electron dynamics" and "attosecond chemistry." Also on the theory side, powerful quantum dynamical tools have been developed to rationalize experiments on ultrafast electron dynamics in molecular species.
In this contribution, we concentrate on theoretical aspects of ultrafast electron dynamics in molecules, mostly driven by lasers. The dynamics will be described with the help of wavefunction-based ab initio methods such as time-dependent configuration interaction (TD-CI) or the multiconfigurational time-dependent Hartree-Fock (MCTDHF) methods. Besides a survey of the methods and their extensions toward, e.g., treatment of ionization, laser pulse optimization, and open quantum systems, two specific examples of applications will be considered: The creation and/or dynamical fate of electronic wavepackets, and the nonlinear optical response to laser pulse excitation in molecules by high harmonic generation (HHG).}, language = {en} } @article{TremblayMonturetSaalfrank2010, author = {Tremblay, Jean Christophe and Monturet, Serge and Saalfrank, Peter}, title = {Electronic damping of anharmonic adsorbate vibrations at metallic surfaces}, issn = {1098-0121}, doi = {10.1103/Physrevb.81.125408}, year = {2010}, abstract = {The nonadiabatic coupling of an adsorbate close to a metallic surface leads to electronic damping of adsorbate vibrations and line broadening in vibrational spectroscopy. Here, a perturbative treatment of the electronic contribution to the lifetime broadening serves as a building block for a new approach, in which anharmonic vibrational transition rates are calculated from a position-dependent coupling function. Different models for the coupling function will be tested, all related to embedding theory. The first two are models based on a scattering approach with (i) a jellium-type and (ii) a density functional theory based embedding density, respectively. In a third variant a further refined model is used for the embedding density, and a semiempirical approach is taken in which a scaling factor is chosen to match harmonic, single-site, first-principles transition rates, obtained from periodic density functional theory. For the example of hydrogen atoms on (adsorption) and below (subsurface absorption) a Pd(111) surface, lifetimes of and transition rates between vibrational levels are computed. The transition rates emerging from different models serve as input for the selective subsurface adsorption of hydrogen in palladium starting from an adsorption site, by using sequences of infrared laser pulses in a laser distillation scheme.}, 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} }