State Potential Energy Surface Research Articles

Strongly bound anions featuring bismuth fluoride building blocks

The stability of polynuclear superhalogen anions composed of BiF5 building blocks was investigated using ab initio electronic structure methods and flexible basis sets. A comprehensive exploration of the ground state potential energy surfaces of (Bi2F11)−, (Bi3F16)− and (Bi4F21)− anions, which can be viewed as comprising BiF5 fragments and an additional fluorine atom, led to the identification of their isomeric structures. It was found that the most stable isomers, predicted to dominate at room temperature, correspond to chain-like extended structures containing BiF6 subunits, with fluorine ligands arranged octahedrally around Bi atoms, sharing F atoms to form Bi–F–Bi bridging linkages. The vertical electron detachment energies of the (BinF5n+1)− anions (n = 1–4) were found to be very high (ranging from 10.91 to 13.36 eV) and increased with the number of bismuth atoms (n) and thus the BiF5 building blocks involved in the structure. Thermodynamic stability of the (BinF5n+1)− anions (i.e., their susceptibility to fragmentation) was also verified and discussed.

Photoisomerization of heptamethine cyanine (Cy7) dyes: A theoretical study

In this study, density functional theory (DFT) combined with time-dependent (TD) DFT calculations were employed to investigate the photoisomerization reaction kinetics of two near infrared (NIR) heptamethine cyanine (Cy7-NH3 and Cy7-SO3) dyes in the ground singlet state and the first excited singlet state. We found that the photoisomerization of the ground state all-trans Cy7 molecules results in at least one mono-cis and one all-cis species that demonstrate redshifted emission, in agreement with recently published transient state excitation modulation spectroscopy and fluorescence correlation spectroscopy measurements. The transition states were estimated for a whole photoisomerization pathway for both the ground singlet and first excited singlet state potential energy surfaces. We have found that all-cis isomers of the studied Cy7 dyes can be achieved through a sequential two-step photoisomerization within the excited singlet state potential energy surface, along the double CC bond adjacent to edge group (leading to mono-cis isomer 1) and along the double CC bond adjacent to the central-chain group (leading to mono-cis isomer 2). Computations show that all-trans→ mono-cis isomer 1→all-cis kinetics is limited by the first trans→ mono-cis isomer 1 stage, while the all-trans→ mono-cis isomer 2→all-cis pathway is limited by the second mono-cis isomer 2→all-cis stage. Accounting for the fact that mono-cis isomer 2 demonstrates red-shifted emission compared to the all-trans form and that this mono-cis isomer 2 is reachable through the energetically favorable all-trans→ mono-cis isomer 2 stage, we concluded that the experimentally observed red-shifted emission by Cy7-NH3 and Cy7-SO3 should be assigned to the formation of mono-cis isomer 2 species. If the all-cis isomer is populated through the double-step photoisomerization it can also be considered as a source of red-shifted emission. However, as follows from our simulations, the all-cis isomer is kinetically intricate to achieve compared to the mono-cis isomer 2.

Kinetic Study and Reaction Mechanism of the Gas-Phase Thermolysis Reaction of Methyl Derivatives of 1,2,4,5-Tetroxane.

Tetroxane derivatives are interesting drugs for antileishmaniasis and antimalaric treatments. The gas-phase thermal decomposition of 3,6,-dimethyl-1,2,4,5-tetroxane (DMT) and 3,3,6,6,-tetramethyl-1,2,4,5-tetroxane (acetone diperoxide (ACDP)) was studied at 493-543 K by direct gas chromatography by means of a flow reactor. The reaction is produced in the injector chamber at different temperatures. The resulting kinetics Arrhenius equations were calculated for both tetroxanes. Including the parent compound of the series 1,2,4,5-tetroxane (formaldehyde diperoxide (FDP)), the activation energy and frequency factors decrease linearly with the number of methyl groups. The reaction mechanisms of ACDP and 3,6,6-trimethyl-1,2,4,5-tetroxane (TMT) decomposition have been studied by means of the DFT method with the BHANDHLYP functional. Our calculations confirm that the concerted mechanism should be discarded and that only the stepwise mechanism occurs. The critical points of the singlet and triplet state potential energy surfaces (S- and T-PES) of the thermolysis reaction of both compounds have been determined. The calculated activation energies of the different steps vary linearly with the number of methyl groups of the methyl-tetroxanes series. The mechanism for the S-PES leads to a diradical O···O open structure, which leads to a C···O dissociation in the second step and the production of the first acetaldehyde/acetone molecule. This last one yields a second C···O dissociation, producing O2 and another acetone/acetaldehyde molecule. The O2 molecule is in the singlet state. A quasi-parallel mechanism for the T-PES from the open diradical to products is also found. Most of the critical points of both PES are linear with the number of methyl groups. Reaction in the triplet state is much more exothermic than the singlet state mechanism. Transitions from the singlet ground state, S0 and low-lying singlet states S1-3, to the low-lying triplet excited states, T1-4, (chemical excitation) in the family of methyl tetroxanes are also studied at the CASSCF/CASPT2 level. Two possible mechanisms are possible here: (i) from S0 to T3 by strong spin orbit coupling (SOC) and subsequent fast internal conversion to the excited T1 state and (ii) from S0 to S2 from internal conversion and subsequent S2 to T1 by SOC. From these experimental and theoretical results, the additivity effect of the methyl groups in the thermolysis reaction of the methyl tetroxane derivatives is clearly highlighted. This information will have a great impact for controlling these processes in the laboratory and chemical industries.

Novel Insights into the Catalytic Mechanism of Collagenolysis by Zn(II)-Dependent Matrix Metalloproteinase-1.

Collagen hydrolysis, catalyzed by Zn(II)-dependent matrix metalloproteinases (MMPs), is a critical physiological process. Despite previous computational investigations into the catalytic mechanisms of MMP-mediated collagenolysis, a significant knowledge gap in understanding remains regarding the influence of conformational sampling and entropic contributions at physiological temperature on enzymatic collagenolysis. In our comprehensive multilevel computational study, employing quantum mechanics/molecular mechanics (QM/MM) metadynamics (MetD) simulations, we aimed to bridge this gap and provide valuable insights into the catalytic mechanism of MMP-1. Specifically, we compared the full enzyme-substrate complex in solution, clusters in solution, and gas-phase to elucidate insights into MMP-1-catalyzed collagenolysis. Our findings reveal significant differences in the catalytic mechanism when considering thermal effects and the dynamic evolution of the system, contrasting with conventional static potential energy surface QM/MM reaction path studies. Notably, we observed a significant stabilization of the critical tetrahedral intermediate, attributed to contributions from conformational flexibility and entropy. Moreover, we found that protonation of the scissile bond nitrogen occurs via proton transfer from a Zn(II)-coordinated hydroxide rather than from a solvent water molecule. Following C-N bond cleavage, the C-terminus remains coordinated to the catalytic Zn(II), while the N-terminus forms a hydrogen bond with a solvent water molecule. Subsequently, the release of the C-terminus is facilitated by the coordination of a water molecule. Our study underscores the pivotal role of protein conformational dynamics at physiological temperature in stabilizing the transition state of the rate-limiting step and key intermediates, compared to the corresponding reaction in solution. These fundamental insights into the mechanism of collagen degradation provide valuable guidance for the development of MMP-1-specific inhibitors.

Uncovering the sensing mechanism of a zinc ion sensor: Fluorescence enhancement induced by the elimination of the TICT state

Precise detection of zinc ion is of fundamental importance in the fields of environment protection and food safety. A comprehensive understanding of the sensing mechanism will help to the design of such sensors. The detailed photophysical process of a zinc ion sensor as well as the sensing mechanism are uncovered with the aid of density functional theory (DFT) and time-dependent density functional theory (TDDFT). Both the ground state and first excited state potential energy surfaces (PES) of the sensor are carefully explored to reveal the photo-physical process of the sensor. Excited state intramolecular proton transfer (ESIPT) is observed on the S1 state PES. Then, the twist motion of C=N double bond is triggered after the ESIPT process, which leads to a twisted intramolecular charge transfer (TICT) state. This TICT state is found to make the sensor non-emissive. With the addition of Zn2+, the TICT state is eliminated which greatly enhances the fluorescence of the sensor and achieves zinc ion detection. The interaction of the sensor with Cd2+ and Hg2+ are also explored, which well explains the good selectivity of the sensor.

Ultrafast terahertz Stark spectroscopy reveals the excited-state dipole moments of retinal in bacteriorhodopsin

The photoinduced all-trans to 13-cis isomerization of the retinal Schiff base represents the ultrafast first step in the reaction cycle of bacteriorhodopsin (BR). Extensive experimental and theoretical work has addressed excited-state dynamics and isomerization via a conical intersection with the ground state. In conflicting molecular pictures, the excited state potential energy surface has been modeled as a pure S[Formula: see text] state that intersects with the ground state, or in a 3-state picture involving the S[Formula: see text] and S[Formula: see text] states. Here, the photoexcited system passes two crossing regions to return to the ground state. The electric dipole moment of the Schiff base in the S[Formula: see text] and S[Formula: see text] state differs strongly and, thus, its measurement allows for assessing the character of the excited-state potential. We apply the method of ultrafast terahertz (THz) Stark spectroscopy to measure electric dipole changes of wild-type BR and a BR D85T mutant upon electronic excitation. A fully reversible transient broadening and spectral shift of electronic absorption is induced by a picosecond THz field of several megavolts/cm and mapped by a 120-fs optical probe pulse. For both BR variants, we derive a moderate electric dipole change of 5 [Formula: see text] 1 Debye, which is markedly smaller than predicted for a neat S[Formula: see text]-character of the excited state. In contrast, S[Formula: see text]-admixture and temporal averaging of excited-state dynamics over the probe pulse duration gives a dipole change in line with experiment. Our results support a picture of electronic and nuclear dynamics governed by the interaction of S[Formula: see text] and S[Formula: see text] states in a 3-state model.

The reaction of methylidyne with methane: role of the quartet electronic state

The reaction of methylidyne (CH) with methane, which plays an important role in combustion and planetary atmospheres, is characterised with high accuracy coupled-cluster calculations. The reaction on the quartet state potential energy surface (PES) is predicted to lead to direct H-abstraction with a small barrier of 2.2 kcal mol−1, yielding triplet CH2 and CH3. The reaction on the doublet state PES is barrierless, involving an insertion/elimination mechanism that leads to the vibrationally excited (hot) intermediate C2H5, which promptly dissociates to H + C2H4. The results obtained from a statistical chemical kinetics analysis show that the quartet channel makes a negligible contribution to the overall kinetics below 500 K, but becomes important at the higher temperatures characteristic of hydrocarbon flames. Reaction rate coefficients calculated from first principles are in good agreement (within 20%) with experimental values where they are available (200–800 K). For the higher temperatures of combustion environments (where experimental data are lacking), the present study provides high-level results for kinetics modelling purposes.

Impact of random nanoscale roughness on gas-scattering dynamics.

The impact of nanoscale wall roughness on rarefied gas transport is widely acknowledged, yet the associated scattering dynamics largely remain elusive. In this paper, we develop a scattering kernel for surfaces having nanoscale roughness that distinctly characterizes the two major types of interactions between gas molecules and rough surfaces. Namely these are (a) the weak perturbations arising from the thermal motion of wall atoms, essentially gas-phonon collisions, which are captured by the well-established Cercignani-Lampis model, and (b) the hard collisions owing to the irregularities of the rough, static potential energy surface, which are generally described by the fully diffuse model. Drawing an analogy between wave-surface and gas-surface scattering, a pseudo Debye-Waller factor is incorporated into the modeling as a weighting coefficient to allow the transition between smooth and rough surface conditions. The proposed scattering kernel is validated through high-fidelity molecular dynamics simulations that are performed for systems with varying roughness, temperature, and gas-surface combinations. The results indicate that the model well captures the scattering dynamics of gas molecular beams impinging on surfaces at different velocities, specifically for the accommodation coefficients and reflection patterns. Additionally, in flow and heat transport cases, it accurately predicts macroscopic quantities such as velocity slip and temperature jumps across the range of tested conditions.

Theoretical study of the photoisomerization of 1,2-bispyrazinyl-ethylene and the halogen ion salts of 1-Pyrazinyl-2-(4’-methylpyrazinyl)ethylene

ContextIt has been reported that photoexcitation of azastilbene compounds like E-1,2-bispyrazinyl-ethylene (bpe) can undergo E → Z photoisomerization of its quaternary salts via the excited triplet state. However, experimentally it is possible to get low fluorescence and photoisomerisation quantum yields in a state with higher internal conversion than intersystem crossing. We modelled bpe and its methylated derivative (bpeMe), as well as its quaternary halogen salts (bpeMeX with X = F−, Cl−, Br− and I−) to study levels of fluorescence, phosphorescence and excited state potential energy surfaces (PES). Results support experimental observations of molecules where the anion of a salt is an efficient electron donor, that molecules with weak electron-donating anions like Cl− to give increased fluorescence and photoisomerization, as compared to molecules with stronger electron-donating anions like I−, which are dominated by competing electron transfer. The fluorescence of bpeMeF and bpeMeCl was found to be stronger than bpeMeBr and bpeMeI. A deep well in the triplet excited state of bpeMeI is considered responsible for the decreased photoisomerization, compared to what was experimentally observed for bpeMeCl. Uniquely, the bpeMeI molecule is characterised by near-zero splitting of the s1 and t1 excited states that can enhance charge transfer. The quaternary salt of bpeMe with stronger electron-donating Br− anion was observed to undergo fluorescence and phosphorescence at much lower energy compared to those with weak electron-donating F− and Cl− anions. This research shows how to control the excited state fluorescence, phosphorescence and isomerization of quaternary halogen salts of methyl derivatives of 1,2-bispyrazinyl-ethylene, which aids experimental design where excited state isomerization is considered.MethodsGeometry optimization, molecular electrostatic potential (MESP), and time dependent density functional theory (TDDFT) calculations were conducted utilizing Gaussian 16 with the B3LYP functional and the 6–31 + G(d,p) basis set. The minimum energy path (MEP) for the E to Z isomerization of the molecules was established employing the Nudged-Elastic-Band (NEB) method, implemented in Orca 4.2. Precise energies of the E → Z isomerization reaction path were determined employing CASSCF and a more accurate multireference method, NEVPT2.

Multiple retinal isomerizations during the early phase of the bestrhodopsin photoreaction

Bestrhodopsins constitute a class of light-regulated pentameric ion channels that consist of one or two rhodopsins in tandem fused with bestrophin ion channel domains. Here, we report on the isomerization dynamics in the rhodopsin tandem domains of Phaeocystis antarctica bestrhodopsin, which binds all-trans retinal Schiff-base (RSB) absorbing at 661 nm and, upon illumination, converts to the meta-stable P540 state with an unusual 11-cis RSB. The primary photoproduct P682 corresponds to a mixture of highly distorted 11-cis and 13-cis RSB directly formed from the excited state in 1.4 ps. P673 evolves from P682 in 500 ps and contains highly distorted 13-cis RSB, indicating that the 11-cis fraction in P682 converts to 13-cis. Next, P673 establishes an equilibrium with P595 in 1.2 µs, during which RSB converts to 11-cis and then further proceeds to P560 in 48 µs and P540 in 1.0 ms while remaining 11-cis. Hence, extensive isomeric switching occurs on the early ground state potential energy surface (PES) on the hundreds of ps to µs timescale before finally settling on a metastable 11-cis photoproduct. We propose that P682 and P673 are trapped high up on the ground-state PES after passing through either of two closely located conical intersections that result in 11-cis and 13-cis RSB. Co-rotation of C11=C12 and C13=C14 bonds results in a constricted conformational landscape that allows thermal switching between 11-cis and 13-cis species of highly strained RSB chromophores. Protein relaxation may release RSB strain, allowing it to evolve to a stable 11-cis isomeric configuration in microseconds.

Beyond MD17: the reactive xxMD dataset

System specific neural force fields (NFFs) have gained popularity in computational chemistry. One of the most popular datasets as a bencharmk to develop NFF models is the MD17 dataset and its subsequent extension. These datasets comprise geometries from the equilibrium region of the ground electronic state potential energy surface, sampled from direct adiabatic dynamics. However, many chemical reactions involve significant molecular geometrical deformations, for example, bond breaking. Therefore, MD17 is inadequate to represent a chemical reaction. To address this limitation in MD17, we introduce a new dataset, called Extended Excited-state Molecular Dynamics (xxMD) dataset. The xxMD dataset involves geometries sampled from direct nonadiabatic dynamics, and the energies are computed at both multireference wavefunction theory and density functional theory. We show that the xxMD dataset involves diverse geometries which represent chemical reactions. Assessment of NFF models on xxMD dataset reveals significantly higher predictive errors than those reported for MD17 and its variants. This work underscores the challenges faced in crafting a generalizable NFF model with extrapolation capability.

Quantum dynamics study of C—H stretching vibrational excitation in the F+CHD<sub>3</sub> → HF+CD<sub>3</sub> reaction

In recent decades, significant progress has been made in the precise theoretical investigation of gas-phase chemical reactions. Presently, a major challenge in the field of quantum dynamics is to develop the precise methodologies for studying chemical reactions involving more than four atoms. As a typical multi-atomic reaction system, the F+CH<sub>4</sub> reaction and its isotopic substitution reactions have attracted widespread attention from both experimental and theoretical perspectives in recent years. Experimental studies on the reaction of F+CHD<sub>3</sub> have revealed that the stretching vibration excitation of the C—H bond inhibits the bond dissociation, favoring the formation of DF+CHD<sub>2</sub> product channels. In this study, we use a seven-dimensional quantum time-dependent wave packet method to investigate the dynamics of the F+CHD<sub>3</sub> reaction in both the reactant vibrational ground state and the first stretching excited state of the C—H bond. In this work, the reaction probabilities under different vibrational conditions are analyzed, showing that when the collision energy is below 0.06 eV, the reaction probability curves exhibit numerous fast-oscillating peaks, supporting the experimentally suggested phenomena of dynamic resonance. In a collision energy range from 0.06 eV to 0.3 eV, the reaction probability for the HF product channel in the vibrational excited state is lower than that in the ground state, which is consistent with experimental observation. Through the analysis of the time-independent wave functions of product channels under low-energy collision conditions, it is found that for reactions involving vibrational ground states, the HF products in the product asymptotic region and the reaction transition state region are in the <i>v'</i> = 2 excited state and <i>v'</i> = 3 excited state of stretching vibration, respectively, which are consistent with previous experimental observations and six-dimensional quantum wave packet simulations. For reactions involving the first excited state of C—H stretching vibration, the HF products in the product asymptotic region and the reaction transition state region are both in the <i>v'</i> = 3 excited state of stretching vibration, which are consistent with the results obtained based on energy analysis. Simulation results indicate that in the case of low-energy collisions, the time-independent wave function for the C—H stretching vibrational excited state tends to be closer to the D atom side in the transition state region. This phenomenon is attributed to the more significant energy advantage of the vibrational excited state potential energy surface in the large collision angle region, explaining the inhibitory effect of stretching vibration excitation on the HF product channel. This study offers important theoretical support for explaining experimental results and contributes to a more in-depth understanding of the influence of vibrational mode excitations on the dynamical processes in poly-atomic reactions.

Tracking the conical intersection dynamics for the photoinduced Jahn-Teller switch of a Mn(iii) complex.

Octahedral Mn(iii) ions predominantly exhibit an axially elongated Jahn-Teller (JT) distortion, which is responsible for their large uniaxial magnetic anisotropy. As a result, they are often used in the synthesis of single-molecule magnets (SMMs). Modulation of the JT distortion using femtosecond laser pulses could offer a route to controlling the magnetisation direction, and therefore is promising for the development of data storage devices that work on ultrafast timescales. Photoinduced switching of the distortion from an axially elongated to an axially compressed structure has been demonstrated for various Mn(iii) complexes. However, the dynamics around the region of the conical intersection for the photoinduced JT switch remains unclear. Here, ultrafast transient absorption spectra were recorded for solutions of tris(2,2,6,6-tetramethyl-3,5-heptanedionato)manganese(iii) (Mn(dpm)3) in ethanol to further explore the dynamics of the photoinduced JT switch. We observe the generation of a vibrational wavepacket on the excited state surface, which has a frequency of approximately 155 cm-1 and encompasses a JT-active vibrational mode. This coherent motion is maintained after passage through the conical intersection back to the ground state, which launches wavepackets along the ground state potential energy surface (PES) with frequencies of approximately 180 and 110 cm-1 that we assign to the elongated and compressed state, respectively. Inspection of the relative phases of the frequencies reveals phase shifts that are consistent with a one-mode reaction coordinate, and passes through the conical intersection at 1/4 and 3/4 of the excited state vibrational period. Our results provide direct insights into the non-adiabatic dynamics of Mn(iii) complexes, which can be used to guide the synthesis of optically controlled SMMs.

Combined Quantum Mechanical and Quasi-Classical State-to-State Dynamical Study on the Isotopic Effect in H/D + LiH+/LiD+ → H2/HD/D2 + Li+ Reactions.

Coriolis-coupled quantum mechanical (QM-CC) and quasi-classical trajectory (QCT) calculations are carried out to investigate the dynamics of the H(D) + LiH+(v = 0, j = 0) → H2(HD) (v', j') + Li+ reactions on the ground electronic state potential energy surface reported by Martinazzo et al. (Martinazzo et al., J. Chem. Phys. 2003, 119, 11241). The QM-CC and QCT results at the initial state-selected and state-to-state levels are used to investigate the validity and accuracy of the QCT method for these exoergic barrierless reactions. Furthermore, the QCT method is used to understand the isotopic effects on reaction observables like total and state-to-state integral cross section, differential cross section, product energy disposal, and rate constants of H(D) + LiH+(v = 0, j = 0) → H2(HD) (v', j') + Li+ and H(D) + LiD+(v = 0, j = 0) → HD(D2) (v', j') + Li+ reactions. Attempts are also made to understand the impact of the isotopic substitution on the reaction mechanism. It is observed that QM-CC and QCT results closely follow each other at the initial state-selected and state-to-state levels. Noticeable kinematic effects of reagents on the reactivity and mechanism of the reactions are also observed.

A global 2A″ state potential energy surface for the Al (2P) + O2 (Σg−3) → AlO (2Σ+) + O (3P) reaction based on the doubly hybrid functional XYG3

In this paper, a global and full-dimensional potential energy surface at the 2A″ ground state for the Al + O2 → AlO + O reaction was constructed, for the first time, based on extensive electronic structure calculations using the doubly hybrid functional XYG3 and potential energy surface fittings by neural networks. Details of the reaction paths have been analyzed. It was found that both two intermediates, the cyclic-AlO2 and the linear-OAlO, were able to dissociate to the AlO + O products, and the isomerization process between these two intermediates was controlled by conical intersections between two 2A″ states. Ro-vibrational state resolved integral cross sections have also been calculated at collision energies from 1.0 to 10.0 kcal/mol. The results support the harpooning mechanism in this metal-oxidant-involved reaction.

Ab initio molecular dynamics study of intersystem crossing dynamics for MH2 (M = Si, Ge, Sn, Pb) on spin-pure and spin-mixed potential energy surfaces.

Recently, surface-hopping ab initio molecular dynamics (SH-AIMD) simulations have come to be used to discuss the mechanisms and dynamics of excited-state chemical reactions, including internal conversion and intersystem crossing. In dynamics simulations involving intersystem crossing, there are two potential energy surfaces (PESs) governing the motion of nuclei: PES in a spin-pure state and PES in a spin-mixed state. The former gives wrong results for molecular systems with large spin-orbit coupling (SOC), while the latter requires a potential gradient that includes a change in SOC at each point, making the computational cost very high. In this study, we systematically investigate the extent to which the magnitude of SOC affects the results of the spin-pure state-based dynamics simulations for the hydride MH2 (M = Si, Ge, Sn, Pb) by performing SH-AIMD simulations based on spin-pure and spin-mixed states. It is clearly shown that spin-mixed state PESs are indispensable for the dynamics simulation of intersystem crossing in systems containing elements Sn and Pb from the fifth period onward. Furthermore, in addition to the widely used Tully's fewest switches (TFS) algorithm, the Zhu-Nakamura (ZN) global switching algorithm, which is computationally less expensive, is applied to SH for comparison. The results from TFS- and ZN-SH-AIMD methods are in qualitative agreement, suggesting that the less expensive ZN-SH-AIMD can be successfully utilized to investigate the dynamics of photochemical reactions based on quantum chemical calculations.

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods.

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.

Fantastical excited state optimized structures and where to find them.

The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximation, where the nuclear and electronic degrees of freedom separate. This work demonstrates through abinitio calculations and simple nonadiabatic dynamics that some excited state minimum structures are fantastical: they appear to exist as stable configurations only as a consequence of the PES construct, rather than being physically observable. Each fantastical structure exhibits an unphysically high predicted harmonic frequency and associated force constant. This fact can serve as a valuable diagnostic of when an optimized excited state structure is non-observable. The origin of this phenomenon can be attributed to the coupling between different electronic states. As PESs approach one another, the upper surface can form a minimum that is very close to a near-touching point. The force constant, evaluated at this minimum, relates to the strength of the electronic coupling rather than to any characteristic excited state vibration. Nonadiabatic dynamics results using a Landau-Zener model illustrate that fantastical excited state structures have extremely short lifetimes on the order of a few femtoseconds. Their appearance in a calculation signals the presence of a nearby conical intersection through which the system will rapidly cross to a lower surface.

Strongly Bound Polynuclear Anions Comprising Scandium Fluoride Building Blocks.

The stability of polynuclear anions composed of ScF3 building blocks was studied by using ab initio and density functional theory electronic structure methods and flexible basis sets. Thorough exploration of ground state potential energy surfaces of (Sc2F7)-, (Sc3F10)-, and (Sc4F13)- anions which may be viewed as comprising ScF3 fragments and the additional fluorine atom led to determining the isomeric structures thereof. It was found that the most stable isomers which are predicted to dominate at room temperature correspond to the compact structures enabling the formation of a large number of Sc-F-Sc bridging linkages rather than to the chain-like structures. The vertical electron detachment energies of the (ScnF3n+1)- anions were found to be very large (spanning the 10.85-12.29 eV range) and increasing with the increasing number of scandium atoms (n) and thus the ScF3 building blocks involved in the structure. Thermodynamic stability of (ScnF3n+1)- anions (i.e., their susceptibility to fragmentation) was also verified and discussed.

Computational chemistry‐assisted design of hydrazine‐based fluorescent molecular rotor for viscosity sensors

AbstractDeep understanding of the fluorescence quenching mechanisms of probes plays a crucial role in developing their practical applications. The fluorescence quenching mechanism of hydrazine‐based fluorescence probes needs to be clarified up to the present. Herein, we designed and synthesized a new hydrazine‐based fluorescence probe (HA‐Na) based on the naphthalimide skeleton. We clarified the molecular origin of the non‐fluorescence of this probe with the aid of computational chemistry and spectroscopic analysis. We showed that the significant rotation of the hydrazine group in the excited state potential energy surface, which caused the complete charge separation, was responsible for the fluorescence quenching of the probe in an organic solvent. Once the rotation was prevented in an aggregative state or high‐viscosity solution, the fluorescence of the probe recovered. In other words, the fluorescence quenching mechanism of hydrazine‐based fluorescence probes was attributed to the formation of a twisted intramolecular charge transfer (TICT) state. More importantly, we demonstrated that this fluorescence molecular rotor could be used to monitor the autophagy process in living cells by detecting lysosomal viscosity changes during starvation. Altogether, this work provides an essential theoretical basis for the developing potential hydrazine‐based fluorescence molecular rotors.