Photoionization-induced dynamics of ammonia: ab initio potential energy surfaces and time-dependent wave packet calculations for the ammonia cation

An analytical anharmonic six-dimensional three-sheeted potential energy surface for the ground and first excited states of the ammonia cation has been developed which is tailored to model the ultrafast photoinduced dynamics. Selected ab initio cuts, obtained by multireference configuration interaction calculations, have been used to determine the parameters of a diabatic representation for this Jahn-Teller and pseudo-Jahn-Teller system. The model includes higher-order coupling terms both for the Jahn-Teller and for the pseudo-Jahn-Teller matrix elements. The relaxation to the ground state is possible via dynamical pseudo-Jahn-Teller couplings involving the asymmetric bending and stretching coordinates. The photoelectron spectrum of NH3 and the internal conversion dynamics of NH3+ have been determined by wave packet propagation calculations employing the multiconfigurational time-dependent Hartree method. Three different time scales are found in the dynamics calculations for the second absorption band. The ultrafast Jahn-Teller dynamics of the two excited states occurs on a 5 fs time scale. The major part of the internal conversion to the ground state takes place within a short time scale of 20 fs. This fast internal conversion is, however, incomplete and the remaining excited state population does not decay completely even within 100 fs.     

Ultrafast anisotropy dynamics of water molecules dissolved in acetonitrile

Infrared pump-probe experiments are performed on isolated H(2)O molecules diluted in acetonitrile in the spectral region of the OH stretching vibration. The large separation between water molecules excludes intermolecular interactions, while acetonitrile as a solvent provides substantial hydrogen bonding. Intramolecular coupling between symmetric and asymmetric modes results in the anisotropy decay to the frequency-dependent values of approximately 0-0.2 with a 0.2 ps time constant. The experimental data are consistent with a theoretical model that includes intramolecular coupling, anharmonicity, and environmental fluctuations. Our results demonstrate that intramolecular processes are essential for the H(2)O stretching mode relaxation and therefore can compete with the intermolecular energy transfer in bulk water.     

Controlling the isomerization dynamics of iodide acetonitrile dimer complex by optimally designed electromagnetic field: A wave packet based approach

Both symmetric and asymmetric forms of the iodide acetonitrile dimer complex are known to exist with a higher stability shown by the symmetric isomer. Dynamics of this isomerization has been investigated on reduced dimensional surface in absence/presence of external electromagnetic field using Gaussian wave packet as the starting wave function. When the symmetric isomer has been taken as the initial structure, isomerization does not takes place without the external field due to high energy barrier, whereas with the asymmetric isomer as the starting one, isomerization takes place spontaneously with an equal population of both isomers. Our target in this study has been to get a dominant population of asymmetric isomer which is not possible without the application of electromagnetic field. However, it has also been seen that a simple monochromatic field cannot achieve this goal. Only by application of optimally designed polychromatic field the dominating population of asymmetric isomer can be obtained.     

Infrared activity of symmetric bridging ligand modes in pyrazine-bridged hexaruthenium mixed-valence clusters

A fully symmetric (A(g)) vibrational mode of pyrazine is observed in the infrared spectrum of four pyrazine-bridged hexaruthenium mixed-valence complexes with varying degrees of electronic coupling between clusters. Deuteration of the bridging pyrazine ligand and the accompanying shift in frequency confirm the assignment of this mode. Previous observation of infrared line coalescence in the carbonyl stretching region assigns all of these complexes to Robin-Day class II (partial localization of charge) on the picosecond time scale. The infrared activity of the fully symmetric bridging ligand mode could provide a complementary assignment of these complexes to class II on a faster, femtosecond time scale. However, the extinction coefficient for this band is much greater than that observed in similar asymmetric, non-mixed-valence complexes and suggests that its strong IR activity is due to vibronic enhancement rather than electronic asymmetry.     

Ultrafast dynamics of halogens in rare gas solids

We perform time resolved pump-probe spectroscopy on small halogen molecules ClF, Cl2, Br2, and I2 embedded in rare gas solids (RGS). We find that dissociation, angular depolarization, and the decoherence of the molecule is strongly influenced by the cage structure. The well ordered crystalline environment facilitates the modelling of the experimental angular distribution of the molecular axis after the collision with the rare gas cage. The observation of many subsequent vibrational wave packet oscillations allows the construction of anharmonic potentials and indicate a long vibrational coherence time. We control the vibrational wave packet revivals, thereby gaining information about the vibrational decoherence. The coherence times are remarkable larger when compared to the liquid or high pressure gas phase. This fact is attributed to the highly symmetric molecular environment of the RGS. The decoherence and energy relaxation data agree well with a perturbative model for moderate vibrational excitation and follow a classical model in the strong excitation limit. Furthermore, a wave packet interferometry scheme is applied to deduce electronic coherence times. The positions of those cage atoms, excited by the molecular electronic transitions are modulated by long living coherent phonons of the RGS, which we can probe via the molecular charge transfer states.     

Photodecay dynamics of octaethylporphine in the condensed phase explored via resonance Raman spectroscopy and density functional theory calculation

Resonance Raman spectra of free-base octaethylporphine (OEP) were obtained with 368.9 nm, 397.9 nm and 416.0 nm excitation wavelengths, and density functional calculations were done to help the elucidation of Soret (B(x) and B(y)-band) electronic transitions and the corresponding photo relaxation dynamics of OEP. The RRs indicate that the Franck-Condon region photo relaxation dynamics upon S(0)→S(8) electronic transition is predominantly along the totally symmetric C(m)C(α) stretch, the C(β)C(β) stretch, and simultaneously along the asymmetric δ(pyr deformation),γ(CH(2)) vibrational relaxation processes. The excited state structural dynamics of OEP determined from resonance Raman spectra show that the internal conversion between B(y) and B(x) electronic states occurs in tens of femtoseconds and the electronic relaxation dynamics were firstly interpreted with account of the time-dependent wave packet theory and Herzberg-Teller (vibronic coupling) contributions.     

Time-resolved intraband electronic relaxation dynamics of Hg(n) (-) clusters (n=7-13,15,18) excited at 1.0 eV

Time-resolved photoelectron imaging has been used to study the relaxation dynamics of small Hg(n) (-) clusters (n=7-13,15,18) following intraband electronic excitation at 1250 nm (1.0 eV). This study furthers our previous investigation of single electron, intraband relaxation dynamics in Hg(n) (-) clusters at 790 nm by exploring the dynamics of smaller clusters (n=7-10), as well as those of larger clusters (n=11-13,15,18) at a lower excitation energy. We measure relaxation time scales of 2-9 ps, two to three times faster than seen previously after 790 nm excitation of Hg(n) (-), n=11-18. These results, along with size-dependent trends in the absorption cross-section and photoelectron angular distribution anisotropy, suggest significant evolution of the cluster anion electronic structure in the size range studied here. Furthermore, the smallest clusters studied here exhibit 35-45 cm(-1) oscillations in pump-probe signal at earliest temporal delays that are interpreted as early coherent nuclear motion on the excited potential energy surfaces of these clusters. Evidence for evaporation of one or two Hg atoms is seen on a time scale of tens of picoseconds.     

Dynamic light scattering in sols of a poly(N-vinylcaprolactam) and polyacrylamide mixture

The properties of semidilute solutions of a mixture of polyacrylamide and thermoresponsive polyvinylcaprolactam are studied via dynamic light scattering. It is found that the relaxation time (correlation length) distribution contains several modes. The obtained experimental data are interpreted in terms of the theory describing the dynamics of semidilute and concentrated solutions with allowance for viscoelastic interactions. The investigation of mode dispersion (the dependence of relaxation time on the wave vector of scattering) shows that the fastest and slowest modes are of the diffusion type. The fast mode is the mode of collective diffusion; it may be attributed to the fluctuation motion of polymer units inside blobs. The slowest mode corresponds to the diffusion motion of macromolecule clusters. The rate of relaxation of medium modes nonlinearly depends on the squared wave vector, and these modes are apparently viscoelastic.     

Resonance Raman spectroscopy and density functional theory calculation study of photodecay dynamics of tetra(4-carboxyphenyl) porphyrin

The 397.9 nm, 416.0 nm and 435.7 nm resonance Raman spectra were acquired for meso-tetrakis(4-carboxyphenyl)porphyrin (TCPP) in tetrahydrofuran solution, and the Raman effect of relaxation dynamics was analyzed according to Herzberg-Teller (vibronic coupling) contributions. Density functional calculations were done to help the elucidation of the Soret (B(x) and B(y)-band) electronic transitions and the corresponding photo relaxation dynamics of TCPP. The spectra indicate that the Franck-Condon region photo relaxation dynamics upon S(0) → S(4) electronic transition are predominantly along the totally symmetric C(m)-ph stretch and Porphin ring breath stretch, and simultaneously along the asymmetric ν(C(m)-Phenyl) + δ(N-H) and ν(C(α)-C(m)-C(α))(as) + def (pyr) vibrational relaxation processes. The excited state structural dynamics of TCPP determined from the resonance Raman spectra show that the internal conversion between the B(y) and B(x) electronic states occurs in tens of femtoseconds, and the electronic relaxation dynamics were firstly interpreted taking into account the time-dependent wave packet theory and Herzberg-Teller (vibronic coupling) contributions.     

Relaxation of bimolecular layer films on water surfaces

Ferric stearate, a three-tailed amphiphile, forms bimolecular layers on water surfaces. Molecules in the lower layer are in an "asymmetric" configuration, Fe-containing heads touching water and three hydrocarbon tails in air, while molecules in the upper layer are in a "symmetric" configuration, in pairs of "Y and inverted Y" disposition of tails about the Fe-bearing head. Pressure relaxation at constant area (pi- t curves) and area relaxation at constant pressure ( A- t curves) of this bimolecular layer can be modeled as a sum of three exponential decay terms with distinct time constants and weight factors. Relating the long-term decay with desorption of the total film thus indicates a remarkable long-term stability of the bimolecular layer film. An X-ray reflectivity study of the bimolecular films deposited horizontally on Si(001) at various conditions of relaxation shows no further growth along the vertical of any other layer. Under pressure relaxation molecules are transferred from the upper layer to the lower layer with a change from symmetric to asymmetric configuration, while under area relaxation the transfer is from the lower layer to the upper layer with a configurational change from symmetric to asymmetric.     

Ultrafast laser studies of the photothermal properties of gold nanocages

Au nanocages were synthesized via a galvanic replacement reaction. The extinction peak of these hollow structured particles is shifted into the near-IR compared with the Ag nanocube templates. Energy transfer from the Au nanocages into the surrounding environment (water) as well as the coherently excited vibrational modes of the nanocages were studied by femtosecond pump-probe spectroscopy. The time scale for energy relaxation was found to increase with the size of the particles, with the relaxation time being independent of the laser intensity. The time scales for relaxation are comparable to those for solid spherical gold particles and are consistent with energy relaxation being controlled by heat dissipation in the solvent. The period of the coherently excited vibrational mode is proportional to the dimensions of the nanocages. Intensity-dependent measurements show that in solution the nanocages maintain their integrity up to lattice temperatures of 1100 +/- 100 K.     

Quantum dynamics study of fulvene double bond photoisomerization: the role of intramolecular vibrational energy redistribution and excitation energy

The double bond photoisomerization of fulvene has been studied with quantum dynamics calculations using the multi-configuration time-dependent Hartree method. Fulvene is a test case to develop optical control strategies based on the knowledge of the excited state decay mechanism. The decay takes place on a time scale of several hundred femtoseconds, and the potential energy surface is centered around a conical intersection seam between the ground and excited state. The competition between unreactive decay and photoisomerization depends on the region of the seam accessed during the decay. The dynamics are carried out on a four-dimensional model surface, parametrized from complete active space self-consistent field calculations, that captures the main features of the seam (energy and locus of the seam and associated branching space vectors). Wave packet propagations initiated by single laser pulses of 5-25 fs duration and 1.85-4 eV excitation energy show the principal characteristics of the first 150 fs of the photodynamics. Initially, the excitation energy is transferred to a bond stretching mode that leads the wave packet to the seam, inducing the regeneration of the reactant. The photoisomerization starts after the vibrational energy has flowed from the bond stretching to the torsional mode. In our propagations, intramolecular energy redistribution (IVR) is accelerated for higher excess energies along the bond stretch mode. Thus, the competition between unreactive decay and isomerization depends on the rate of IVR between the bond stretch and torsion coordinates, which in turn depends on the excitation energy. These results set the ground for the development of future optical control strategies.     

Water as a molecular hinge in amidelike structures

The authors have studied the reorientational dynamics of isolated water molecules in a solution of N,N-dimethylacetamide (DMA). From linear spectra, the authors find that the water in this solution forms double hydrogen bond connections to the DMA molecules, resulting in the formation of DMA-water-DMA complexes. The authors use polarization-resolved mid-infrared pump-probe spectroscopy on the water in these complexes to measure the depolarization of three distinct transition dipole moments, each with a different directionality relative to the molecular frame (OH stretch in HDO, symmetric and asymmetric stretch normal modes in H(2)O). By combining these measurements, the authors find that the system exhibits bimodal rotational dynamics with two distinct time scales: a slow (7+/-1 ps) reorientation of the entire DMA-water complex and a fast (0.5+/-0.2 ps) "hinging" motion of the water molecule around the axis parallel to the connecting hydrogen bonds. Additionally, the authors observe an exchange of energy between the two normal modes of H(2)O at a time scale of 0.8+/-0.1 ps and find that the vibrational excitation decays through the symmetric stretch normal mode with a time constant of 0.8+/-0.2 ps.     

Electron-nuclear correlations for photo-induced dynamics in molecular dimers

Ultrafast photoinduced dynamics of electronic excitation in molecular dimers is drastically affected by the dynamic reorganization of inter- and intra- molecular nuclear configuration modeled by a quantized nuclear degree of freedom. The dynamics of the electronic population and nuclear coherence is analyzed by solving the chain of coupled differential equations for population inversion, electron-vibrational correlation, etc. Intriguing results are obtained in the approximation of a small change of the nuclear equilibrium upon photoexcitation. In the limiting case of resonance between the electronic energy gap and the frequency of the nuclear mode these results are justified by comparison to the exactly solvable Jaynes-Cummings model. It is found that the photoinduced processes in the model dimer are arranged according to their time scales: (i) Fast scale of nuclear motion, (ii) intermediate scale of dynamical redistribution of electronic population between excited states as well as growth and dynamics of electron-nuclear correlation, (iii) slow scale of electronic population approach to the quasi-equilibrium distribution, decay of electron-nuclear correlation, and decrease of the amplitude of mean coordinate oscillation. The latter processes are accompanied by a noticeable growth of the nuclear coordinate dispersion associated with the overall nuclear wave packet width. The demonstrated quantum relaxation features of the photoinduced vibronic dynamics in molecular dimers are obtained by a simple method, applicable to systems with many degrees of freedom.     

The time scale for electronic reorganization upon sudden ionization of the water and water-methanol hydrogen bonded dimers and of the weakly bound NO dimer

When the valence molecular orbital is localized sudden ionization can cause the nascent hole to move rapidly even before any relaxation of the geometry occurs. Hydrogen bonded clusters offer suitable test systems where the hole is initially localized on one moiety. Computational studies are reported for the water dimer and water-methanol bimer. The local ionization potential of water is different in the methanol-water and water-methanol conformers and this difference is very clearly reflected in the dynamics of charge migration. For the NO dimer the results are that its structure is symmetric so that the two NO molecules are equivalent and do not exhibit the required localization. The role of symmetry is also evident in the charge propagation for holes created in different orbitals. Localization of the initial hole distribution even if absent in the bare molecule can still be induced by the intense electric field of a sudden photoionization. This effect is computationally studied for the NO dimer in the presence of a static electric field.     

Vibrational relaxation of normal and deuterated liquid nitromethane

Anti-Stokes Raman scattering is used to monitor vibrational energy redistribution in the ambient temperature liquids nitromethane (NM-h3) and perdeuterated nitromethane (NM-d3) after ultrafast IR excitation of either the symmetric or asymmetric CH- or CD-stretch transitions. The instantaneous populations of most of the fifteen NM vibrations are determined with good accuracy, and a global fitting procedure with a master equation is used to fit all the data. The pump pulses excite not only CH- or CD-stretches but also certain combinations of bending and nitro stretching fundamentals. The coupled vibrations that comprise the initial state are revealed via the instantaneous rise of the anti-Stokes transients associated with each vibrational fundamental. In contrast to many other polyatomic liquids studied previously, there is little energy exchange among the CH-stretch (or CD-stretch) excitations, which is attributed to the nearly free rotation of the methyl group in NM. The vibrational cooling process, which is the multistep return to a thermalized state, occurs in three stages in both NM-h3 and NM-d3. In the first stage, the parent CH- or CD-stretch decays in a few picoseconds, exciting all lower-energy vibrations. In the second stage, the midrange vibrations decay in 10-15 ps, exciting the lower-energy vibrations. In the third stage, these lower-energy vibrations decay into the bath in tens of picoseconds. The initial excitations are thermalized in approximately 150 ps in NM-h3 and there is little dependence on which CH-stretch is excited. VC is somewhat faster in NM-d3 with more dependence on the initial CD-stretch, taking approximately 100 ps with symmetric CD-stretch excitation and approximately 120 ps with asymmetric CD-stretch excitation. Comparison is made with earlier nonequilibrium molecular dynamics simulations of VC [Kabadi, V. N.; Rice, B. M. Molecular dynamics simulations of normal mode vibrational energy transfer in liquid nitromethane. J. Phys. Chem. A 2004, 108, 532-540]. The simulations do a good job of reproducing the observed VC process and in addition they predicted the slow interconversion among CH-stretch excitations and the slower relaxation of the asymmetric CH-stretch now observed here.     

Photodissociation of CH3I: a full-dimensional (9D) quantum dynamics study

The photodissociation of methyl iodide in the A band is studied by full-dimensional (9D) wave packet dynamics calculations using the multiconfigurational time-dependent Hartree approach. The potential energy surfaces employed are based on the diabatic potentials of Xie et al. [J. Phys. Chem. A 2000, 104, 1009] and the vertical excitation energy is taken from recent ab initio calculations [Alekseyev et al. J. Chem. Phys.2007, 126, 234102]. The absorption spectrum calculated for exclusively parallel excitation agrees well with the experimental spectrum of the A band. The electronic population dynamics is found to be strongly dependent on the motion in the torsional coordinate related to the H(3)-C-I bend, which presumably is an artifact of the diabatic model employed. The calculated fully product state-selected partial spectra can be interpreted based on the reflection principle and suggests strong coupling between the C-I stretching and the H(3)-C-I bending motions during the dissociation process. The computed rotational and vibrational product distributions typically reproduce the trends seen in the experiment. In agreement with experiment, a small but significant excitation of the total symmetric stretching and the asymmetric bending modes of the methyl fragment can be seen. In contrast, the umbrella mode of the methyl is found to be too highly excited in the calculated distributions.     

Measuring the extent and width of internal energy deposition in ion activation using nanocalorimetry

The recombination energies resulting from electron capture by a positive ion can be accurately measured using hydrated ion nanocalorimetry in which the internal energy deposition is obtained from the number of water molecules lost from the reduced cluster. The width of the product ion distribution in these experiments is predominantly attributable to the distribution of energy that partitions into the translational and rotational modes of the water molecules that are lost. These results are consistent with a singular value for the recombination energy. For large clusters, the width of the energy distribution is consistent with rapid energy partitioning into internal vibrational modes. For some smaller clusters with high recombination energies, the measured product ion distribution is narrower than that calculated with a statistical model. These results indicate that initial water molecule loss occurs on the time scale of, or faster than energy randomization. This could be due to inherently slow internal conversion or it could be due to a multi-step process, such as initial ion-electron pair formation followed by reduction of the ion in the cluster. These results provide additional evidence for the accuracy with which condensed phase thermochemical values can be deduced from gaseous nanocalorimetry experiments.     

Optical coherence and theoretical study of the excitation dynamics of a highly symmetric cyclophane-linked oligophenylenevinylene dimer

Optoelectronic properties of a polyphenylenevinylene-based oligomer and its paracylophane-linked dimer are studied using a variety of experimental and theoretical techniques. Despite the symmetrical structure and redshifted absorption of the dimer versus the monomer, an exciton picture is not the most appropriate. Electronic structure calculations establish changes in charge density upon optical excitation and show localized excitations that cannot be accounted for by a simple Frenkel exciton model. Visible frequency pump-probe anisotropy measurements suggest that the dimer should be considered as a three-level system with a fast, approximately 130 fs, internal conversion from the higher to lower energy excited electronic state. Signatures of nuclear relaxation processes are compared for electric field-resolved transient grating and two-dimensional photon echo spectra. These measurements reveal that nuclear relaxation occurs on similar time scales for the monomer and dimer. The connection between the spectral phase of four-wave mixing signals and the time dependent width of a nuclear wave packet is discussed. Semiempirical electronic structure and metropolis Monte Carlo calculations show that the dominant line broadening mechanisms for the monomer and dimer are associated with inter-ring torsional coordinates. Together, the theoretical calculations and electric field-resolved four-wave mixing experiments suggest that while the structure of dimer is more rigid than that of monomer, the difference in their rigidities is not sufficient to slow down excited state relaxation of dimer with respect to the monomer.     

Car-Parrinello molecular dynamics simulation of Fe 3+ (aq)

The optimized geometry and energetic properties of Fe(D2O)n 3+ clusters, with n = 4 and 6, have been studied with density-functional theory calculations and the BLYP functional, and the hydration of a single Fe 3+ ion in a periodic box with 32 water molecules at room temperature has been studied with Car-Parrinello molecular dynamics and the same functional. We have compared the results from the CPMD simulation with classical MD simulations, using a flexible SPC-based water model and the same number of water molecules, to evaluate the relative strengths and weaknesses of the two MD methods. The classical MD simulations and the CPMD simulations both give Fe-water distances in good agreement with experiment, but for the intramolecular vibrations, the classical MD yields considerably better absolute frequencies and ion-induced frequency shifts. On the other hand, the CPMD method performs considerably better than the classical MD in describing the intramolecular geometry of the water molecule in the first hydration shell and the average first shell...second shell hydrogen-bond distance. Differences between the two methods are also found with respect to the second-shell water orientations. The effect of the small box size (32 vs 512 water molecules) was evaluated by comparing results from classical simulations using different box sizes; non-negligible effects are found for the ion-water distance and the tilt angles of the water molecules in the second hydration shell and for the O-D stretching vibrational frequencies of the water molecules in the first hydration shell.