Nanosecond simulations of the dynamics of C60 excited by intense near-infrared laser pulses: Impulsive Raman excitation, rearrangement, and fragmentation

Capturing the photo-induced dynamics of nano-molecules by X-ray free electron laser induced Coulomb explosion

We performed reaction dynamics simulations to demonstrate that the vibrational dynamics of C₆₀ induced by infrared (IR) pulses can be traced by triggering Coulomb explosion with intense femtosecond X-ray free electron laser (XFEL) probe pulses. The time series of the angular anisotropy β(t) of fast C⁺ and C²⁺ fragments of C₆₀ ⁶⁰⁺ produced by such an XFEL pulse reflects the instantaneous structure of C₆₀ vibrationally excited by IR pulses. The phases and amplitudes of excited vibrational modes and the coupling between excited modes can be successfully extracted from the expansion of β(t) in terms of vibrational modes. This proof-of-principle simulation clearly demonstrates that various information of the structures and reaction dynamics of large clusters or biomolecules can be retrieved by decomposing the experimentally determined β(t) into vibrational modes.

Diffractive Imaging of C_{60} Structural Deformations Induced by Intense Femtosecond Midinfrared Laser Fields

Theoretical studies indicated that C_{60} exposed to linearly polarized intense infrared pulses undergoes periodic cage structural distortions with typical periods around 100 fs (1  fs=10^{-15}  s). Here, we use the laser-driven self-imaging electron diffraction technique, previously developed for atoms and small molecules, to measure laser-induced deformation of C_{60} in an intense 3.6  μm laser field. A prolate molecular elongation along the laser polarization axis is determined to be (6.1±1.4)% via both angular- and energy-resolved measurements of electrons that are released, driven back, and diffracted from the molecule within the same laser field. The observed deformation is confirmed by density functional theory simulations of nuclear dynamics on time-dependent adiabatic states and indicates a nonadiabatic excitation of the h_{g}(1) prolate-oblate mode. The results demonstrate the applicability of laser-driven electron diffraction methods for studying macromolecular structural dynamics in four dimensions with atomic time and spatial resolutions.

Ultrafast Coulomb explosion of a diiodomethane molecule induced by an X-ray free-electron laser pulse

The Coulomb explosion mechanism of a CH₂I₂ molecule is rather different to that of CH₃I. The kinetic energy of iodine ions is ∼3 times larger due to Coulomb repulsion of the two iodine ions, while that of carbon ions is almost the same for both, as indicated by the red arrows that represent kinetic energies of the atomic ions.

Femtosecond charge and molecular dynamics of I-containing organic molecules induced by intense X-ray free-electron laser pulses

We studied the electronic and nuclear dynamics of I-containing organic molecules induced by intense hard X-ray pulses at the XFEL facility SACLA in Japan. The interaction with the intense XFEL pulse causes absorption of multiple X-ray photons by the iodine atom, which results in the creation of many electronic vacancies (positive charges) via the sequential electronic relaxation in the iodine, followed by intramolecular charge redistribution. In a previous study we investigated the subsequent fragmentation by Coulomb explosion of the simplest I-substituted hydrocarbon, iodomethane (CH₃I). We carried out three-dimensional momentum correlation measurements of the atomic ions created via Coulomb explosion of the molecule and found that a classical Coulomb explosion model including charge evolution (CCE-CE model), which accounts for the concerted dynamics of nuclear motion and charge creation/charge redistribution, reproduces well the observed momentum correlation maps of fragment ions emitted after XFEL irradiation. Then we extended the study to 5-iodouracil (C₄H₃IN₂O₂, 5-IU), which is a more complex molecule of biological relevance, and confirmed that, in both CH₃I and 5-IU, the charge build-up takes about 10 fs, while the charge is redistributed among atoms within only a few fs. We also adopted a self-consistent charge density-functional based tight-binding (SCC-DFTB) method to treat the fragmentations of highly charged 5-IU ions created by XFEL pulses. Our SCC-DFTB modeling reproduces well the experimental and CCE-CE results. We have also investigated the influence of the nuclear dynamics on the charge redistribution (charge transfer) using nonadiabatic quantum-mechanical molecular dynamics (NAQMD) simulation. The time scale of the charge transfer from the iodine atomic site to the uracil ring induced by nuclear motion turned out to be only ∼5 fs, indicating that, besides the molecular Auger decay in which molecular orbitals delocalized over the iodine site and the uracil ring are involved, the nuclear dynamics also play a role for ultrafast charge redistribution. The present study illustrates that the CCE-CE model as well as the SCC-DFTB method can be used for reconstructing the positions of atoms in motion, in combination with the momentum correlation measurement of the atomic ions created via XFEL-induced Coulomb explosion of molecules.

Theoretical studies on a carbonaceous molecular bearing: association thermodynamics and dual-mode rolling dynamics

The thermodynamics and dynamics of a carbonaceous molecular bearing comprising a belt-persistent tubular molecule and a fullerene molecule have been investigated using density functional theory (DFT). Among ten representative methods, two DFT methods afforded an association energy that reasonably reproduced the experimental enthalpy of -12.5 kcal mol-1 at the unique curved π-interface. The dynamics of the molecular bearing, which was assembled solely with van der Waals interactions, exhibited small energy barriers with maximum values of 2-3 kcal mol-1 for the rolling motions. The dynamic motions responded sensitively to the steric environment and resulted in two distinct motions, precession and spin, which explained the unique NMR observations that were not clarified in previous experimental studies.

Roles of dynamical symmetry breaking in driving oblate-prolate transitions of atomic clusters

This paper explores the driving mechanisms for structural transitions of atomic clusters between oblate and prolate isomers. We employ the hyperspherical coordinates to investigate structural dynamics of a seven-atom cluster at a coarse-grained level in terms of the dynamics of three gyration radii and three principal axes, which characterize overall mass distributions of the cluster. Dynamics of gyration radii is governed by two kinds of forces. One is the potential force originating from the interactions between atoms. The other is the dynamical forces called the internal centrifugal forces, which originate from twisting and shearing motions of the system. The internal centrifugal force arising from twisting motions has an effect of breaking the symmetry between two gyration radii. As a result, in an oblate isomer, activation of the internal centrifugal force that has the effect of breaking the symmetry between the two largest gyration radii is crucial in triggering structural transitions into prolate isomers. In a prolate isomer, on the other hand, activation of the internal centrifugal force that has the effect of breaking the symmetry between the two smallest gyration radii is crucial in triggering structural transitions into oblate isomers. Activation of a twisting motion that switches the movement patterns of three principal axes is also important for the onset of structural transitions between oblate and prolate isomers. Based on these trigger mechanisms, we finally show that selective activations of specific gyration radii and twisting motions, depending on the isomer of the cluster, can effectively induce structural transitions of the cluster. The results presented here could provide further insights into the control of molecular reactions.

The influence of distribution of hydroxyl groups on vibrational spectra of fullerenol C60(OH)24 isomers: DFT study

The infrared and Raman spectra of C60(OH)24 molecule with uniform and non-uniform distribution of hydroxyl groups have been investigated using first principle DFT calculations at the B3LYP/6-31G(d,p) level of theory. The important features of the obtained geometries have been measured and compared to experimental results. The reference calculations of C60 molecule geometry and vibrational spectra have been made and compared to available experimental data. The striking differences of infrared spectra between C60(OH)24 molecule with uniform and non-uniform distribution of hydroxyl groups have been shown and discussed. The OH modes have been identified as the most sensitive to C60(OH)24 isomer configuration. The C-C stretching modes in the Raman spectra of the C60(OH)24 molecule have been found as a potential sensor of OH groups distribution over fullerene C60 surface.