Excited Potential Energy Surfaces Research Articles

Exploring how molluscan purple has survived throughout the Ages: The excited state dynamics of 6,6′-dibromoindigotin

The color purple has long been used as a symbol of royalty and power in art, fashion, and architecture. This correlation with aristocracy can be traced back to the use of molluscan purple, a rich purple dye that is still presently more expensive than gold, by royalty as early as the 18th century BCE. While the molluscan purple dye is composed of a mixture of various proteins, its color is derived from indigotin, indirubin and their brominated derivatives, including the molecule, 6,6′-dibromoindigotin, an analogue of the indigotin dye which is also widely used in various modern industries. Both dyes can exhibit high stability, in particular against photodamage from UV and Visible radiation. In this study, we present the gas phase absorption spectrum and excited-state lifetimes of 6,6′-dibromoindigotin combined with static calculations of the excited and ground-state potential energy surfaces. The lifetime measurements reveal that molluscan purple has nearly the same relaxation rate as indigotin providing new insights in the possible relaxation mechanisms of the indigotin family of dye molecules.

Quantum Chemical Studies on Possible Molecular Devices Based on Electric Field-Induced Intramolecular Charge Transfer.

We report quantum chemical studies on possible molecular devices working based on electric field-induced intramolecular charge transfer (EFIMCT). In the case of donor-acceptor (DA)-type molecular systems, intramolecular charge transfer (IMCT) can be induced by applying the external electric field to molecular systems along the charge transport direction, providing a possible switching mechanism which does not depend upon the electron-phonon coupling effect and is different from the negative differential resistance mechanism observed in the well-known NO2-substituted phenylene ethynylene oligomers. When the EFIMCT proceeds, the molecular systems have strong static electron correlation effects, where the standard nonequilibrium Green's function-density functional theory (DFT) approach cannot be applied to the molecular junction. As a first step toward practical switching devices, we do quantum chemical studies on the EFIMCT in such molecular systems as an isolated molecule, instead of using the electrode-junction-electrode open quantum system model. A prototype molecule P1 is designed as a tentative candidate molecule where the EFIMCT can proceed. The complete active space self-consistent field (CASSCF) molecular orbital calculations on P1 indicate that the EFIMCT can proceed at the external electric field intensity of 0.003 au, corresponding to about 2.25 V bias voltage. This calculated result strongly suggests that the development of this type of switching devices working at practically low bias voltage is feasible if the molecular system is properly designed. Broken symmetry unrestricted Hartree-Fock and spin-polarized Kohn-Sham DFT calculations also qualitatively reproduce the CASSCF results on P1, to some extent, indicating that these approaches can be employed for rough estimations on the EFIMCT such as the first screening of a large quantity of candidate molecules for this type of molecular devices. The possibility of molecular memory devices based on the EFIMCT is also discussed by analyzing the ground and excited potential energy surface model. Remaining challenges to develop practical molecular devices are discussed.

Semiclassical Dynamics on Machine-Learned Coupled Multireference Potential Energy Surfaces: Application to the Photodissociation of the Simplest Criegee Intermediate.

Determination of high-dimensional potential energy surfaces (PESs) and nonadiabatic couplings have always been quite challenging. To this end, machine learning (ML) models, trained with a finite set of ab initio data, allow accurate prediction of such properties. To express the PESs in terms of atomic contributions is the cornerstone of any ML based technique because it can be easily scaled to large systems. In this work, we have constructed high fidelity PESs and nonadiabatic coupling terms at the CASSCF level of ab initio data using a machine learning technique, namely, kernel-ridge regression. Additional MRCI-level calculations were carried out to assess the quality of the PESs. We use these machine-learned PESs and nonadiabatic couplings to simulate excited-state molecular dynamics based on Tully's fewest-switches surface hopping method (FSSH). FSSH is a semiclassical method in which nuclei move on the PESs due to the electrons according to the laws of classical mechanics. Nonadiabatic effects are taken into account in terms of transitions between PESs. We apply this scheme to study the O-O photodissociation of the simplest Criegee intermediate (CH2OO). The FSSH trajectories were initiated on the lowest optically bright singlet excited state (S2) and propagated along the three most important internal coordinates, namely, O-O and C-O bond distances and the COO bond angle. Some of the trajectories end up on energetically lower PESs as a result of radiationless transfer through conical intersections. All of the trajectories lead to the dissociation of the O-O bond due to the dissociative nature of the excited PESs through one of the two dissociative channels. The simulation reveals that there is about 88.4% probability of dissociation through the lower channel leading to the H2CO (X1A1) and O (1D) products, whereas there is only 11.6% probability of dissociation through the upper channel leading to H2CO (a3A″) and O (3P) products.

An intramolecular vibrationally excited intermolecular potential energy surface and predicted 2OH overtone spectroscopy of H2O-Kr.

A new five-dimensional potential energy surface (PES) for H2O-Kr which explicitly includes the intramolecular 2OH overtone state of the H2O monomer is presented. The intermolecular potential energies were evaluated using explicitly correlated coupled cluster theory [CCSD(T)-F12] with a large basis set. Four vibrationally averaged analytical intermolecular PESs for H2O-Kr with H2O molecules in its |00+〉, |02+〉, |02-〉, and |11+〉 states are obtained by fitting to the multi-dimensional Morse/Long-Range potential function form. Each vibrationally averaged PES fitted to 578 points has root-mean-square (RMS) deviations smaller than 0.14 cm-1 and requires only 58 parameters. The combined radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm were employed to calculate the rovibrational energy levels for |00+〉, |02+〉, |02-〉, and |11+〉 states of the H2O-Kr complexes. The calculated |02-〉Πf/e(101) ← |00+〉Σe(000) and |02+〉Πf/e(110) ← |00+〉Σe(101) infrared transitions are in excellent agreement with the experimental values with RMS discrepancies being only 0.007 and 0.016 cm-1, respectively. These analytical PESs can be used to provide reliable theoretical guidance for future infrared overtone spectroscopy of H2O-Kr.

Clusters of Fullerenes: Structures and Dynamics.

The geometric structures and reaction dynamics of clusters of carbon fullerene molecules are reviewed. The topics on structure cover the elementary building blocks, the interatomic and intermolecular potentials, and the geometric structures of the aggregates. The dynamics part describes the time development after excitation with laser light, in bimolecular collisions, and in collisions with high energy atomic ions. These reactions produce singly or multiply charged clusters, fragmented or intact fullerene monomers, and fusion products. The duration of the interaction time and the transferred energy are found to play determining roles for the nature of the products. Short interaction times and high energy excitation favor intramolecular reactions on excited potential energy surfaces, producing fused products, often with a strong release of the combined excess energy and heat of fusion. Clusters excited at low energies tend to disintegrate into intact fullerene molecules.

Charge transfer in ultrafast isomerization of acetylene ions

First-principle calculations are employed to investigate the ultrafast isomerization of the acetylene cation and dication. We use the time-dependent density functional theory together with the Ehrenfest dynamics to track the coupled electron-nuclear dynamics. For both the acetylene cation and the dication, we observe nonadiabatic behaviors during the isomerization. We find that the charge transfer not only alters the electronic structure through nonadiabatic transitions, but also plays a key role in the subsequent hydrogen migration. We show that nonadiabatic transitions affect the structural modification of the excited potential energy surface, and also facilitate the ultrafast isomerization through the creation of a channel of increased negative charge that facilitates the proton movement. For the acetylene cation, we find a timescale for hydrogen isomerization of $66\ifmmode\pm\else\textpm\fi{}4$ fs, which is consistent with previous pump-probe experiments and on-the-fly calculations. For the dication, we find nonadiabatic transitions occur before the isomerization and identify a similar channel for the proton. Moreover, we find the formation of vinylidene-like structures is always accompanied by a characteristic charge separation on the carbon skeleton. These heuristics will be useful in identifying tautomers and motivating the ultrafast charge-transfer detection methods for future experiments.

Exploring the Photophysics of a Zn2+ Fluorescence Sensor and Its Sensing Mechanism.

Sensor X is a turn-on sensor, which is applied in the fluorescence detection of Zn2+ ions. Its photophysical process is comprehensively investigated to clarify its weak fluorescence. With the aid of density functional theory (DFT) and time-dependent density functional theory (TDDFT), the potential energy surfaces (PES) of X on both ground and first excited states are studied. Excited-state intramolecular proton transfer (EPT) processes as well as molecule twisting motion are observed, which induces several minima on the excited-state PES. Transition states as well as rate constants for these dynamic processes are obtained to evaluate their occurrences. The twisting motion of the sensor is an ultrafast process, which is initiated by a specific EPT process and leads to a nonemissive twisted intramolecular charge transfer (TICT) state. The fluorescence of the sensor is barely observable because of the easily attainable TICT state on the excited PES. This mechanism is trustworthy and intrinsically different from the previously proposed mechanism. After clarifying the photophysical process of the sensor, the Zn2+ sensing mechanism is uncovered. Also, the selectivity against Cd2+ and Hg2+ is fully discussed.

Coherent Control of Molecular Dissociation by Selective Excitation of Nuclear Wave Packets.

We report on pump-probe control schemes to manipulate fragmentation product yields in p-nitrotoluene (PNT) cation. Strong field ionization of PNT prepares the parent cation in the ground electronic state, with coherent vibrational excitation along two normal modes: the C–C–N–O torsional mode at 80 cm−1 and the in-plane ring-stretching mode at 650 cm−1. Both vibrational wave packets are observed as oscillations in parent and fragment ion yields in the mass spectrum upon optical excitation. Excitation with 650 nm selectively fragments the PNT cation into , whereas excitation with 400 nm selectively produces and . In both cases the ion yield oscillations result from torsional wave packet excitation, but 650 and 400 nm excitation produce oscillations with opposite phases. Ab initio calculations of the ground and excited electronic potential energy surfaces of PNT cation along the C–C–N–O dihedral angle reveal that 400 nm excitation accesses an allowed transition from D0 to D6 at 0° dihedral angle, whereas 650 nm excitation accesses a strongly allowed transition from D0 to D4 at a dihedral angle of 90°. This ability to access different electronic excited states at different locations along the potential energy surface accounts for the selective fragmentation observed with different probe wavelengths. The ring-stretching mode, only observed using 800 nm excitation, is attributed to a D0 to D2 transition at a geometry with 90° dihedral angle and elongated C–N bond length. Collectively, these results demonstrate that strong field ionization induces multimode coherent excitation and that the vibrational wave packets can be excited with specific photon energies at different points on their potential energy surfaces to induce selective fragmentation.

Excited-State Energy Surfaces in Molecules Revealed by Impulsive Stimulated Raman Excitation Profiles.

Photophysical and photochemical processes are ruled by the interplay between transient vibrational and electronic degrees of freedom, which are ultimately determined by the multidimensional potential energy surfaces (PESs). Differences between ground and excited PESs are encoded in the relative intensities of resonant Raman bands, but they are experimentally challenging to access, requiring measurements at multiple wavelengths under identical conditions. Herein, we perform a two-color impulsive vibrational scattering experiment to launch nuclear wavepacket motions by an impulsive pump and record their coupling with a targeted excited-state potential by resonant Raman processes with a delayed probe, generating in a single measurement background-free vibrational spectra across the entire sample absorption. Building on the interference between the multiple pathways resonant with the excited-state manifold that generate the Raman signal, we show how to experimentally tune their relative phase by varying the probe chirp, decoding nuclear displacements along different normal modes and revealing the multidimensional PESs. Our results are validated against time-dependent density functional theory.

Trans-cis Photoisomerization of a biomimetic cyclocurcumin analogue rationalized by molecular modelling.

Cyclocurcumin is a natural compound extracted from turmeric and showing, in addition to antiinfective, antibacterial, and intinflammatory capabilities, solvent-dependent phtoswitching ability. The solvent-dependent photochemistry of cyclocurcumin has been rationalized on the basis of a competition between π-π* and n-π* states. Recently we have reported the synthesis of a biomimetic analogue showing enhanced photochemical properties and in particular presenting photoswitching capacity in various media. In the present contribution we rely on the use of molecular modeling and simulation, incuding density functional and wavefunction based methods to explore the excited states potential energy surface landscape. We see that the addition of a carbon-carbon double bond to the core of the natural compounds favors the population of the π-π* state, whatever the choice of the solvent, and hence leads to photoisomerisation, with fluorescence reduced to only a minor channel, rationalizing the experimental observations. In addition, the two photon absorption cross section is also strongly increased compared to the parent compound, paving the way to the use in biologically oriented applications.

Ab initio MRCI study on potential energy surfaces for double Cl loss from the palladium tetrachloride anion PdCl42−

We theoretically study double Cl loss from the palladium tetrachloride anion PdCl42− using the ab initio MRCI method. The double Cl loss is probably a significant process in the laser-induced particle formation of Pd in aqueous/EtOH solution which was proposed recently by Saeki et al. (M. Saeki et al., J. Phys. Chem. C 123 (2019) 817). It is predicted that the wave packet created on excited potential energy surfaces (PESs) by the UV irradiation at 266 nm temporally propagates on the PESs and finally spreads over a flat region, where neutral Cl and anionic Cl− are generated.

New insight into the ultrafast excited state deactivation mechanism of guanosine in the gas phase

The deactivation process of guanosine molecule on the first singlet excited state in the gas phase has been explored through the density functional theory (DFT/TDDFT) and the complete active space self-consistent field (CASSCF) methods. To unravel the decay mechanism of this system, the topography of excited potential energy surfaces relevant to the deactivation process were mapped. Two deactivation decay channels characterized by the excited state intramolecular proton transfer process and the puckering motion of ring were explored, and the latter one was proposed to play significant roles in the deactivation process of guanosine for the barrierless energy profile. However, the probability for the excited state intramolecular proton transfer process pathway is lower because of the existence of an energy barrier. This conclusion was further confirmed by the two minimum energy conical intersections (MECIs) located between S1 and S0 states, characterized by the ESIPT process and pyramidalization of the C2 atom, respectively. In addition, by comparing the decay mechanism with the previously reported adenosine and pyrimidine nucleosides, it was found that two competitive decay pathways coexist in the purine nucleosides, which is different from the pyrimidine nucleosides. As the next step in the bottom-up investigation of the DNA photostability properties, the new mechanism of the guanosine in gas phase will be conducive for more comprehensive understanding of the photostability of DNA and RNA molecules and provide theoretical guidance for the investigation of the different decay mechanism of the purine nucleosides and pyrimidine nucleosides.

Prospects for laser cooling of polyatomic molecules with increasing complexity

Optical cycling transitions and direct laser cooling have recently been demonstrated for a number of alkaline-earth dimers and trimer molecules. This is made possible by diagonal Franck-Condon factors between the vibrational modes of the optical transition. Achieving a similar degree of cooling to micro-Kelvin equivalent kinetic energy for larger polyatomic molecules, however, remains challenging. Since polyatomic molecules are characterized by multiple degrees of freedom and have a correspondingly more complex structure, it is far from obvious whether there exist polyatomic molecules that can repeatedly scatter photons. Here, we propose chemical substitution approaches to engineer large polyatomic molecules with optical cycling centers (OCCs) containing alkaline-earth oxide dimers or acetylenic alkaline-earth trimers (i.e. M--C$\equiv$C) connected to CH$_n$ chains or fullerenes. To validate the OCC-character of the selected molecules we performed electronic structure calculations of the equilibrium configuration of both ground and excited potential energy surfaces, evaluated their vibrational and bending modes, and corresponding Franck-Condon factors for all but the most complex molecules. For fullerenes, we have shown that OCCs based on M--C$\equiv$C can perform better than those based on alkaline-earth oxide dimers. In addition, for heavier polyatomic molecules it might be advantageous to attach two OCCs, thereby potentially doubling the photon scattering rate and thus speeding up cooling rates.

On the Problem of Multiple Minima on the S2 Excited Potential Energy Surface of Benzene: A Restricted Active Space Self Consistent Field Study

AbstractApplication of the Restricted Active Space Self‐Consistent Field Theory on the S2 excited state potential energy surface of benzene shows the existence of three minima. Of these, the first one is planar with a D6h molecular point group, while the other two have a boat type structure with an approximate C2v symmetry. All three minima have a biradicaloid structure and constitute potential candidates for new photochemical pathways on the S2 surface. A common S2/S1 Conical Intersection, which is accessible from all three S2 minima has been found which decays to the S1 potential energy surface (PES) giving benzene S1 minimum as the only photoproduct.

Effect of isotope on state-to-state dynamics for reactive collision reactions and in ground state 12A″ and first excited 12A′ potential energy surfaces**Project supported by the National Natural Science Foundation of China (Grant No. 11504206) and the Shandong Jiaotong University PhD Research Start-up Fund, China.

We carry out quantum scattering dynamics and quasi-classical trajectory (QCT) calculations for the reactive collision in the ground (12A″) and first excited (12A′) potential energy surface. We calculate the reaction probabilities of and reaction for total angular momentum J = 0. The results calculated by QCT are consistent with those from quantum mechanical wave packet. Using the QCT method, we generate in the center-of-mass frame the product state-resolved integral cross-sections (ICSs); two commonly used generalized polarization-dependent differential cross-sections (PDDCSs), (2π/σ)(dσ00/dωt), (2π/σ)(dσ20/dωt); and three angular distributions of the product rotational vectors, P(θr), P(ϕr), and P(θr, ϕr). We discuss the influence on the scalar and vector properties of the potential energy surface, the collision energy, and the isotope mass. Since there are deep potential wells in these two potential energy surfaces, their kinetic characteristics are similar to each other and the isotopic effect is not obvious. However, the well depths and configurations of the two potential energy surfaces are different, so the effects of isotopic substitution on the integral cross-section and the rotational polarization of product are different.

HCS(A2A″)-based insights into the effect of vibrational excitation on the reactions C+SH (v = 0–20, j = 0) → S+CH, H+CS

The effect of vibrational excitation on reaction C+SH (v = 0–20, j = 0) → S+CH, H+CS is investigated on the excited potential energy surface of HCS(A2A″) by the quasi-classical trajectory method. The obtained reaction probability, total integral cross section (ICS), and the impact parameter show that the influence of vibration excitation presents different characteristics on different reaction channels. The vibrational state-resolved ICSs, differential cross sections as well as two-angle distribution functions P(θr), P(ϕr) of products for different vibrational quantum numbers of reactant are investigated. These results show that (i) the products have obvious forward–backward scattering feature; (ii) for different reactions, the distribution P(θr) varies with vibrational quantum number of reactant; (iii) at high vibrational excitations of the reactant, the insertion mechanism becomes apparent in this reaction, so the product molecules are more positively oriented along the positive direction of the scattering plane.

Electronic spectroscopy of carbon chains (C2n+1, n = 7–10) of astrophysical importance. I. Quantum chemistry

Carbon chains have been predicted to be potential carriers of diffuse interstellar band features in astrophysical observations. Motivated by numerous predictions, we set out to carry out extensive ab initio quantum chemistry calculations to establish the ground and excited electronic potential energy surfaces and their coupling surfaces for carbon chains containing an odd number of carbon atoms (C2n+1, n = 7–10). Vibronic coupling models are established with the aid of the calculated electronic energies to investigate nuclear dynamics from first principles. The latter are reported in Ghosh et al. [J. Chem. Phys. 151, 054304 (2019)]. The mentioned carbon chains possess a linear cumulenic structure at the equilibrium minimum of their electronic ground state, and an electronic excited state of the Σu+1 term appears to be extremely bright optically and absorbs in the visible region of the electromagnetic spectrum. Vertical excitation energy of this state decreases and transition dipole moment increases, and as a result, the oscillator strength of this state linearly increases with an increase of the chain length. There are states belonging to 1Πg, 1Πu, Σg+1, 1Δg, and 1Δu terms, in the immediate vicinity of the Σu+1 state, which are optically dark but can gain intensity through vibronic coupling with the optically bright Σu+1 state. Construction of a coupling scheme considering the Renner-Teller coupling within the degenerate Π states and pseudo-Renner-Teller coupling between the Renner-Teller split component states as well as with the nondegenerate Σ states is another motivation of this work. The coupled-state Hamiltonian is constructed in a diabatic electronic basis in terms of the dimensionless normal coordinates of the vibrational modes of the carbon chains. Both Renner-Teller and pseudo-Renner-Teller types of couplings are included in the Hamiltonian. The theoretical results are discussed in relation to the experimental findings.

Vibronic structure of the cyanobutadiyne cation. II. Theoretical exploration of the complex energy landscape of HC5N.

The results of an extensive ab initio study of the cyanobutadiyne cation, initially motivated by threshold-photoelectron spectroscopy experiments [see the study by Gans et al., J. Chem. Phys. 150, 244304 (2019)], are reported in the present paper. Calculations at the internally contracted multireference configuration interaction level of theory have been performed to derive the rovibronic properties of the seven lowest electronic states of HC5N+. Equilibrium geometries, rotational constants, vibrational frequencies, electric dipole moments, and spin-orbit constants have been calculated and compared with experimental data when available. Adiabatic and vertical ionization energies from the neutral ground state as well as transition energies within the cation electronic manifold are predicted, using the convergence to the complete basis set limit. The accurate description of the complex energy landscape up to 32 000 cm-1 above the ionization potential allows us to perform Franck-Condon simulations of the photoionization spectrum to the X+ 2Π, A+ 2Π, B+ 2Σ+, and C+ 2Π states and allows us to simulate the A+ 2Π → X+ 2Π emission spectrum. The vibronic perturbations occurring on the excited potential energy surfaces are revealed and discussed, in particular, for the 3 2Π surface, which presents a double-well topography.

Electrostatic Spectral Tuning Maps for Biological Chromophores.

The mechanism by which the absorption wavelength of a molecule is modified by a protein is known as spectral tuning. Spectral tuning is often achieved by electrostatic interactions that stabilize/destabilize or modify the shape of the excited and ground-state potential energy surfaces of the chromophore. We present a protocol for the construction of three-dimensional "electrostatic spectral tuning maps" that describe how vertical excitation energies in a chromophore are influenced by nearby charges. The maps are built by moving a charge on the van der Waals surface of the chromophore and calculating the change in its excitation energy. The maps are useful guides for protein engineering of color variants, for interpreting spectra of chromophores that act as probes of their environment, and as starting points for further quantum mechanical/molecular mechanical studies. The maps are semiquantitative and can approximate the magnitude of the spectral shift induced by a point charge at a given position with respect to the chromophore. We generate and discuss electrostatic spectral tuning maps for model chromophores of photoreceptor proteins, fluorescent proteins, and aromatic amino acids. Such maps may be extended to other properties such as oscillator strengths, absolute energies (stability), ionization energies, and electron affinities.

Phase-space wavepacket dynamics of internal conversion via conical intersection: Multi-state quantum Fokker-Planck equation approach

We theoretically investigate internal conversion processes of a photoexcited molecule in a condensed phase. The molecular system is described by two-dimensional adiabatic ground and excited potential energy surfaces that are coupled to heat baths. We quantify the role of conical intersection (CI) and avoided crossing (AC) in the PESs in dissipative environments by simulating the time evolution of wavepackets to compute the lifetime of the excited wavepacket, yield of the product, and adiabatic electronic coherence. For this purpose, we employ the multi-state quantum Fokker-Planck equation (MSQFPE) for a two-dimensional Wigner space utilizing the Wigner–Moyal expansion for the potential term and the Brinkman hierarchy expression for the momentum. We find that the calculated results are significantly different between the CI and AC cases due to the transition in the tuning mode and vibrational motion in the coupling mode.