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.     

Momentum representation of the solute/bath interaction in the dynamic theory of chemical processes in condensed phase

For the system consisting of the chemically reactive solute immersed in the oscillator bath, we consider an approach based on the solute/medium interaction expressed in terms of momenta rather than coordinates. In the adiabatic representation the medium reorganization effects are suppressed, being hidden in the solute renormalized potential and new spectral density function. The advantage proposed by the bilinear interaction in momentum representation is its spatial uniformity important for approximate dynamical treatments. The procedure of explicit transforming a standard spectral density (coordinate representation of interaction) into the spectral density in adiabatic representation (momentum representation of interaction) is the main new result of the present study. Illustrative calculations for several types of spectral functions are performed. Special discussion is devoted to clarifying the nature of the slow diffusion coordinate, to which the present approach is mainly addressed.     

Density functional theory study of 1,2‐dioxetanone decomposition in condensed phase

The decomposition of 1,2‐dioxetanone into a CO₂ molecule and into an excited state formaldehyde molecule was studied in condensed phase, using a density functional theory approach. Singlet and triplet ground and excited states were all included in the calculations. The calculations revealed a novel mechanism for the chemiluminescence of this compound. The triplet excitation can be explained by two intersystem crossings (ISCs) with the ground state, while the singlet excitation can be accounted by an ISC with the triplet state. The experimentally verified small excitation yield can then be explained by the presence of an energy barrier present in the potential energy surface of the triplet excited state, which will govern both triplet and singlet excitation. It was also found that the triplet ground state interacts with both the triplet excited and singlet ground states. A MPWB1K/mPWKCIS approach provided results in agreement with the existent literature. © 2012 Wiley Periodicals, Inc.     

Ultrafast coherent dynamics of nonadiabatically coupled quasi-degenerate excited states in molecules: Population and vibrational coherence transfers

Results of a theoretical study of ultrafast coherent dynamics of nonadiabatically coupled quasi-degenerate π-electronic excited states of molecules were presented. Analytical expressions for temporal behaviors of population and vibrational coherence were derived using a simplified model to clarify the quantum mechanical interferences between the two coherently excited electronic states, which appeared in the nuclear wavepacket simulations [M. Kanno, H. Kono, Y. Fujimura, S.H. Lin, Phys. Rev. Lett 104 (2010) 108302]. The photon-polarization direction of the linearly polarized laser, which controls the populations of the two quasi-degenerate electronic states, determines constructive or destructive interference. Features of the vibrational coherence transfer between the two coupled quasi-electronic states through nonadiabatic couplings are also presented. Information on both the transition frequency and nonadiabatic coupling matrix element between the two states can be obtained by analyzing signals of two kinds of quantum beats before and after transfer through nonadiabatic coupling.     

Transition rates between selected vibrational states in methyl halides; influence of anharmonicities

A semiclassical model of collision induced vibrational relaxation is discussed in terms of an effective collision mass for different values of vibrational energy release. Selected one, two and three quantum transitions of the methyl halides upon collision with rare gases are evaluated in the presence of resonant and nonresonant anharmonic couplings. It is found, that due to the anharmonic coupling the rates between CH stretching modes and the overtones of the CH bending modes become as large as the transfer rates between two CH stretching modes. This is in qualitative agreement with experiments. Without the anharmonic coupling they differ by two orders of magnitude.     

A study of molecular vibrational relaxation mechanism in condensed phase based upon mixed quantum-classical molecular dynamics. II. Noncollisional mechanism for the relaxation of a polar solute in supercritical water

Mixed quantum-classical molecular dynamics method has been applied to vibrational relaxation of a hydrophilic model NO in supercritical water at various densities along an isotherm above the critical temperature. The relaxation rate was determined based on Fermi's golden rule at each state point and showed an inverse S-shaped curve as a function of bulk density. The hydration number was also calculated as a function of bulk density based on the calculated radial distribution function, which showed a good correlation with the relaxation rate. Change of the survival probability of the solute vibrational state was analyzed as a function of time together with the trajectory of the solvent water and the interaction with it. We will show that the solvent molecule resides near the solute molecule for a while and the solvent contributes to the relaxation by the random-noiselike Coulombic interaction only when it stays near the solute. After the solvent leaves the solute, it shows no contribution to the relaxation. The relaxation mechanism for this system is significantly different from the collisional one found for a nonpolar solute in nonpolar solvent in Paper I. Then, the relaxation rate is determined, on average, by the hydration number or local density of the solvent. Thus, the density dependence of the relaxation rate for the polar solute in supercritical water is apparently similar to that found for the nonpolar solute in nonpolar solvent, although the molecular process is quite different from each other.     

Umbrella sampling of solute vibrational line shifts in mixed quantum-classical molecular dynamics simulations

An umbrella sampling approach for vibrational frequency line shifts is presented. The technique allows for efficient sampling of the solvent configurations corresponding to frequency shifts of a solute in mixed quantum-classical simulations. The approach is generally applicable and can also be used within traditional perturbation theory calculations of frequency shifts. It is particularly useful in the extraction of detailed mechanistic information about the solute-solvent interactions giving rise to the frequency shifts. The method is illustrated by application to the simple I2 in a liquid Xe system, and the advantages are discussed.     

Femtosecond quantum control of molecular dynamics in the condensed phase

We review the progress in controlling quantum dynamical processes in the condensed phase with femtosecond laser pulses. Due to its high particle density the condensed phase has both high relevance and appeal for chemical synthesis. Thus, in recent years different methods have been developed to manipulate the dynamics of condensed-phase systems by changing one or multiple laser pulse parameters. Single-parameter control is often achieved by variation of the excitation pulse's wavelength, its linear chirp or its temporal subpulse separation in case of pulse sequences. Multiparameter control schemes are more flexible and provide a much larger parameter space for an optimal solution. This is realized in adaptive femtosecond quantum control, in which the optimal solution is iteratively obtained through the combination of an experimental feedback signal and an automated learning algorithm. Several experiments are presented that illustrate the different control concepts and highlight their broad applicability. These fascinating achievements show the continuous progress on the way towards the control of complex quantum reactions in the condensed phase.     

Comparison between the Landau-Teller and flux-flux methods for computing vibrational energy relaxation rate constants in the condensed phase

The calculation of vibrational energy relaxation (VER) rate constants in the condensed phase is usually based on the Landau-Teller formula, which puts them in terms of the Fourier transform, at the vibrational frequency, of the autocorrelation function of the force exerted on the relaxing mode by the bath modes. An alternative expression for the VER rate constant puts it in terms of the autocorrelation function of the vibrational energy flux. In this paper, we compare the predictions obtained via those two methods in the case of iodine in liquid xenon. We find that the computational cost underlying both methods is comparable and that they predict similar VER rates. However, while the calculation of the VER rate via the Landau-Teller formula is somewhat more direct, the predictions obtained via the flux-flux formula are in somewhat better agreement with the VER rates obtained from nonequilibrium molecular dynamics simulations.     

eQE: An open‐source density functional embedding theory code for the condensed phase

In this work, we present the main features and algorithmic details of a novel implementation of the frozen density embedding (FDE) formulation of subsystem density functional theory (DFT) that is specifically designed to enable ab initio molecular dynamics (AIMD) simulations of large‐scale condensed‐phase systems containing 1000s of atoms. This code (available at http://eqe.rutgers.edu ) has been given the moniker of embedded Quantum ESPRESSO (eQE) as it is a generalization of the open‐source Quantum ESPRESSO (QE) suite of programs. The strengths of eQE reside in a hierarchical parallelization scheme that allows for an efficient and fully self‐consistent treatment of the electronic structure (via the addition of an additional DIIS extrapolation layer) while simultaneously exploiting the inherent symmetries and periodicities in the system (via sampling of subsystem‐specific first Brillouin zones and utilization of subsystem‐specific basis sets). While bulk liquids and molecular crystals are two classes of systems that exemplify the utility of the FDE approach (as these systems can be partitioned into weakly interacting subunits), we show that eQE has significantly extended this regime of applicability by outperforming standard semilocal Kohn–Sham DFT (KS‐DFT) for large‐scale heterogeneous catalysts with quite different layer‐specific electronic structure and intrinsic periodicities. eQE features very favorable strong parallel scaling for a model system of bulk liquid water composed of 256 water molecules, which allows for a significant decrease in the overall time to solution when compared to KS‐DFT. We show that eQE achieves speedups greater than one order of magnitude ( ) when performing AIMD simulations of such large‐scale condensed‐phase systems as: (1) molecular liquids via bulk liquid water represented by 1024 independent water molecules (3072 atoms with a 25.3× speedup over KS‐DFT), (2) polypeptide/biomolecule solvation via (gly )₆ solvated in (H₂O)₃₉₅ (1230 atoms with a 38.6× speedup over KS‐DFT), and (3) molecular crystals via a 3 × 3 × 3 periodic supercell of pentacene (1940 atoms with a 12.0× speedup over KS‐DFT). These results represent a significant improvement over the current state‐of‐the‐art and now enable subsystem DFT‐based AIMD simulations of realistically sized condensed‐phase systems of interest throughout chemistry, physics, and materials science.     

Modeling real dynamics in the coarse-grained representation of condensed phase systems

This work presents a systematic multiscale methodology to provide a more faithful representation of real dynamics in coarse-grained molecular simulation models. The theoretical formalism is based on the recently developed multiscale coarse-graining (MS-CG) method [S. Izvekov and G. A. Voth, J. Phys. Chem. B. 109, 2469 (2005); J. Chem. Phys. 123, 134105 (2005)] and relies on the generalized Langevin equation approach and its simpler Langevin equation limit. The friction coefficients are determined in multiscale fashion from the underlying all-atom molecular dynamics simulations using force-velocity and velocity-velocity correlation functions for the coarse-grained sites. The diffusion properties in the resulting CG Brownian dynamics simulations are shown to be quite accurate. The time dependence of the velocity autocorrelation function is also well-reproduced relative to the all-atom model if sufficient resolution of the CG sites is implemented.     

The theory of atomic-molecular transformations in condensed phase at low temperatures

The probability of the elementary act of a nonadiabatic chemical reaction in condensed medium at low temperature is calculated. General results are illustrated by the formulas for Debye-type and resonance-type absorption spectra of the medium. For the reactions involving the reorganization of the intramolecular structure of the reactants the limiting cases of both weak and strong coupling of the reactants with the medium are considered.     

吲哚分子振动光谱的密度泛函理论研究

用多种密度泛函理论(DFT)方法(BLYP/6-31G^*^*,B3LYP/6-31G^*^*,B3PW91/6-31G^*^*和SVWN/6-31G^*^*)对吲哚分子的平衡几何构型进行了优化。在优化构型的基础上计算了吲哚分子的谐力场、振动基频和红外光谱强度。计算得到的振动频率与实验值比较平均偏差对四种计算方法(BLYP/6-31G^*^*,P3LYP/6-31G^*^*,B3PW91/6-31G^*^*和SVWN/6-31G^*^*)分别为16.3,40.5,45.1和26.4cm^-^1。BLYP/6-31G^*^*理论力场被用于吲哚分子的简正坐标分析计算中。根据振动率的势能分布(PEDs)对此分子的振动基频进行了理论归属。     

On a fokker-planck type approach to molecular dynamics in condensed phase

A fast algorithm for dealing with multiplicative stochastic processes such as dielectric relaxation is developed. This technique allows us to get a straightforward evaluation of the effects of non-linearity, expecially non-Gaussian behaviour. Its use jointly to computer “experiments” on non-Gaussian behaviour is shown to be of basic importance for establishing a new approach to molecular dynamics in terms of non-linear models.     

A study of molecular vibrational relaxation mechanism in condensed phase based upon mixed quantum-classical molecular dynamics. I. A test of IBC model for the relaxation of a nonpolar solute in nonpolar solvent at high density

In order to investigate vibrational relaxation mechanism in condensed phase, a series of mixed quantum-classical molecular dynamics calculations have been executed for nonpolar solute in nonpolar solvent and polar solute in polar solvent. In the first paper (Paper I), relaxation mechanism of I2 in Ar, where Lennard-Jones force is predominant in the interaction, is investigated as a function of density and temperature, focusing our attention on the isolated binary collision (IBC) model. The model was originally established for the relaxation in gas phase. A key question, here, is "can we apply the IBC model to the relaxation in the high-density fluid?" Analyzing the trajectory of solvent molecule as well as its interaction with the solute, we found that collisions between them may be defined clearly even in the high-density fluid. Change of the survival probability of the vibrationally first excited state on collision was traced. The change caused by collisions with a particular solvent molecule was also traced together with the interaction between them. Each collision makes a contribution to the relaxation by a stepwise change in the probability. The analysis clearly shows that the relaxation is caused by collisions even in the high-density fluid. The difference between stepwise relaxation and the continuous one found for the total relaxation in the low-density fluid and in the high-density one, respectively, was clarified to come from just the difference in frequency of the collision. The stronger the intensity of the collision is, the greater the relaxation caused by the collision is. Further, the shorter the collision time is, the greater the resultant relaxation is. The discussion is followed by the succeeding paper (Paper II), where we report that molecular mechanism of the relaxation of a polar molecule in supercritical water is significantly different from that assumed in the IBC model despite that the density dependence of the relaxation rate showed a linear correlation with the local density of water around the solute, the linear correlation being apparently in good accordance with the IBC model. The puzzle will be solved in Paper II.     

Ab initio centroid path integral molecular dynamics: application to vibrational dynamics of diatomic molecular systems

An ab initio centroid molecular dynamics (CMD) method is developed by combining the CMD method with the ab initio molecular orbital method. The ab initio CMD method is applied to vibrational dynamics of diatomic molecules, H2 and HF. For the H2 molecule, the temperature dependence of the peak frequency of the vibrational spectral density is investigated. The results are compared with those obtained by the ab initio classical molecular dynamics method and exact quantum mechanical treatment. It is shown that the vibrational frequency obtained from the ab initio CMD approaches the exact first excitation frequency as the temperature lowers. For the HF molecule, the position autocorrelation function is also analyzed in detail. The present CMD method is shown to well reproduce the exact quantum result for the information on the vibrational properties of the system.     

Pump-probe spectroscopy with phase-locked pulses in the condensed phase: decoherence and control of vibrational wavepackets

Electronic and vibrational coherences of Cl2 embedded in solid Ar are investigated by exciting to the B state with a phase-locked pulse pair from an unbalanced Michelson interferometer, where the chirp difference matches the B state anharmonicity. Recording the A' --> X fluorescence after relaxation is compared to probing to charge transfer states by a third pulse. The three-pulse experiment delivers more details on the decoherence processes. The signal modulation due to phase tuning up to the third vibrational round-trip time indicates that the electronic coherence in the B <-- X transition is preserved for more than 660 fs in the solid Ar environment where many body electronic interactions take place. Vibrational coherence lasts longer than 3 ps according to the observed half revival of the wavepacket. Control of the coupling between wavepacket motion and lattice oscillation is demonstrated by tuning the relative phase between the phase-locked pulses, preparing wavepackets predominantly composed of either zero-phonon lines or phonon side bands.     

Hybrid molecular dynamics-quantum mechanics simulations of solute spectral properties in the condensed phase: evaluation of simulation parameters

We have explored the impact of a number of basic simulation parameters on the results of a recently developed hybrid molecular dynamics-quantum mechanics (MD-QM) method (Mercer et al., J Phys Chem B 1999, 103, 7720). The method utilizes MD simulations to explore the ground-state configuration space of the system and QM evaluation of those structures to yield the time-dependent electronic transition energy, which is transformed into the optical line-broadening function using the second-order cumulant expansion. Both linear and nonlinear optical spectra can then be generated for comparison to experiment. The dependence of the resulting spectra on the length of the MD trajectory, the QM sampling rate, and the QM model chemistry have all been examined. In particular, for the system of oxazine-4 in methanol studied here, at least 20 ps of MD trajectory are needed for qualitative convergence of linear spectral properties, and >100 ps is needed for quantitative convergence. Surprisingly, little difference is found between the 3-21G and 6-31G(d) basis sets, and the CIS and TD-B3LYP methods yield remarkably similar spectra. The semiempirical INDO/s method yields the most accurate results, reproducing the experimental Stokes shift to within 5% and the FWHM to within 20%. Nonlinear 3-pulse photon echo peak shift (3PEPS) decays have also been simulated. Decays are generally poorly reproduced, though the initial peak shift which depends on the overall coupling of motions to the solute transition energy is within 15% of experiment for all model chemistries other than those using the STO-3G basis.     

Density functional theory study of vibrational spectra of acridine and phenazine

Density functional theory (DFT) calculations using Becke's exchange in conjunction with Lee-Yang-Parr's correlation functionals (BLYP), Becke's three-parameter hybrid DFT/HF method using Lee-Yang-Parr's correlation functionals (B3LYP) and ab initio Hartree-Fock (HF) method have been carried out to investigate the structure and vibrational spectra of acridine and phenazine. Structural parameters obtained by B3LYP/6-31G* geometry optimization are in good agreement with available experimental data. The raw BLYP non-CH stretching frequencies approximate the experimental results much better than the HF results with the mean absolute deviation about 16 cm(-1). The scaled B3LYP frequencies are more reliable than that of the BLYP and HF methods with the mean absolute deviation about 17 cm(-1). On the basis of the comparison between calculated and experimental results, assignments of fundamental vibrational modes are examined. Also the structure and vibrational frequencies are compared with those of anthracene, pyridine and benzene to study the similarities and differences.