Search

The ab initio multiple spawning method was used to calculate the absorption and resonance Raman excitation profiles of ethylene after varPi-varPi excitation.

Full multiple spawning method is used in investigating the photoinduced cis-trans isomerization of all-trans retinal protonated Schiff base (RPSB). The study seeks to find out if equilibrium is a condition for directed excitation transport given an electronic funnel. RPSB is known for fast quenching of less than 500 fs which guarantees that equilibrium among the excited electronic states cannot be achieved before electronic quenching. Results show quenching and exciton transfer are two distinct events although the latter is a highly directed event.

We describe a new multistate empirical valence bond (MS-EVB) model of OH−in aqueous solutions. This model is based on the recently proposed “charged ring” parameterization for the intermolecular interaction of hydroxyl ion with water [Ufimtsev, et al., Chem. Phys. Lett., 2007, 442, 128] and is suitable for classical molecular simulations of OH−solvation and transport. The model reproduces the hydration structure of OH−(aq) in good agreement with experimental data and the results of ab initiomolecular dynamics simulations. It also accurately captures the major structural, energetic, and dynamic aspects of the proton transfer processes involving OH−(aq). The model predicts an approximately two-fold increase of the OH−mobility due to proton exchange reactions. [ABSTRACT FROM AUTHOR]

It is generally accepted that the anomalous diffusion of the aqueous hydroxide ion results from its ability to accept a proton from a neighboring water molecule; yet, many questions exist concerning the mechanism for this process. What is the solvation structure of the hydroxide ion? In what way do water hydrogen bond dynamics influence the transfer of a proton to the ion? We present the results of femtosecond pump-probe and 2D infrared experiments that probe the O-H stretching vibration of a solution of dilute HOD dissolved in NaOD/D2O. Upon the addition of NaOD, measured pump-probe transients and 2D IR spectra show a new feature that decays with a 110-fs time scale. The calculation of 2D IR spectra from an empirical valence bond molecular dynamics simulation of a single NaOH molecule in a bath of H2O indicates that this fast feature is due to an overtone transition of Zundel-like H3O+ states, wherein a proton is significantly shared between a water molecule and the hydroxide ion. Given the frequency of vibration of shared protons, the observations indicate the shared proton state persists for 2-3 vibrational periods before the proton localizes on a hydroxide. Calculations based on the EVB-MD model argue that the collective electric field in the proton transfer direction is the appropriate coordinate to describe the creation and relaxation of these Zundel-like transition states. [ABSTRACT FROM AUTHOR]

Excited state intramolecular proton transfer (ESIPT) and electronic relaxation in o-hydroxybenzaldehyde (OHBA) are simulated using ab initio multiple spawning (AIMS) dynamics, and the results are compared with those obtained for the closely related molecules malonaldehyde and methyl salicylate. The role of electronic state ordering (more specifically, the placement of spectroscopically dark states) is examined in the context of tautomerisation yield. When a state of nπ* character lies between the optically accessible state and S0, the extent of ESIPT is hindered. [ABSTRACT FROM AUTHOR]

Abstract: The effect of vibrational Zero Point Energy (ZPE) on the torsional barriers of trans-stilbene is studied in the adiabatic approximation. The two torsional modes corresponding to phenyl rotation are explicitly separated, and the remaining modes are treated as normal coordinates. ZPE reduces the adiabatic barrier along the in-phase torsion from 198 to 13cm−1. A one-dimensional adiabatic potential for the anti-phase torsion, including the ZPE of the in-phase torsion, reduces the adiabatic barrier from 260 to 58cm−1. Comparison with recent electronic structure benchmark calculations suggests that vibrational corrections play a significant role in trans-stilbene’s experimentally observed planar structure. [Copyright &y& Elsevier]

Examines the role of intersection topography in bond selectivity of cis-trans photoisomerization. Use of ab initio method to characterize the ground; Impact of the absence of specific interactions in explaining the observed bond selectivity; Attribution of the origin of bond selectivity.

We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions. [ABSTRACT FROM AUTHOR]

Jiang and co-workers use high resolution STM to investigate the reaction and self-assembly of (3,6-dibromo-9,10-phenanthrenequinone, or DBPQ) molecules on Ag (100) and Ag (110) surfaces in order to understand the mechanism of bottom-up assembly on surfaces.[31] They show that, through the inclusion of multiple functional groups within a precursor molecule, it becomes possible to fabricate new low-dimensional materials with unique chemical, physical, and electronic properties. Herbst and Fransson consider the core-valence separation approximation that is often used in the calculation of core-level spectra.[5] They show how to quantify the errors in this approximation, thereby opening the door to error-quantified predictions relevant to x-ray spectroscopy. 153(16), 164108 (2020).10.1063/5.0019557 5 M. F. Herbst and T. Fransson, "Quantifying the error of the core-valence separation approximation", J. Chem. Phys. Zhu and co-workers tackle this problem for a model system containing a 2D semiconductor heterojunction and show convincingly the efficient hot electron transfer from photoexcited MoTe SB 2 sb to WS SB 2 sb .[30] This finding provides important insight into the competition between hot electron cooling and transfer at 2D semiconductor interfaces and suggests an intriguing possibility for the exploration of hot electron devices. [Extracted from the article]

The article discusses a study which examines the isomerization of molecule retinal in rhodopsin passing through a canonical intersection seam using optical spectroscopy. It involves observation of the molecule's stimulated light emission and absorption, and evaluation of its energy difference between excited and ground electronic states. It also performs several computational simulations of excited retinal. Results of the study show the rapid de-excitation of retinal caused by the seam.

We demonstrate the existence of a strictly non-adiabatic control pathway in deprotonation of the acetylene dication. This pathway is identified experimentally by measuring a kinetic energy shift in an ion coincidence experiment. We use a time dependent Schrödinger equation simulation to identify which properties most strongly affect our control. We find that resonant control around conical intersections is limited by the speed of non-adiabatic dynamics. [ABSTRACT FROM AUTHOR]

Abstract: We propose a new ‘charged ring’ (CR) model to better represent the delocalized electron density of the oxygen atom of hydrated OH− in classical molecular simulations. This approach corrects a major defect of traditional point-charge models, which inaccurately predict that OH− accepts predominantly 5–6 hydrogen bonds in aqueous solution. Classical molecular dynamics simulations with the CR model of OH− result in ∼4 accepted H-bonds, in agreement with recent neutron diffraction data and ab initio molecular dynamics results. The observed competition between 3- and 4-coordinate OH−(aq) structures is also in good agreement with ab initio results. [Copyright &y& Elsevier]

Dynamic effects are an important determinant of chemical reactivity and selectivity, but the deliberate manipulation of atomic motions during a chemical transformation is not straightforward. Here, we demonstrate that extrinsic force exerted upon cyclobutanes by stretching pendant polymer chains influences product selectivity through force-imparted nonstatistical dynamic effects on the stepwise ring-opening reaction. The high product stereoselectivity is quantified by carbon-13 labeling and shown to depend on external force, reactant stereochemistry, and intermediate stability. Computational modeling and simulations show that, besides altering energy barriers, the mechanical force activates reactive intramolecular motions nonstatistically, setting up “flyby trajectories” that advance directly to product without isomerization excursions. A mechanistic model incorporating nonstatistical dynamic effects accounts for isomer-dependent mechanochemical stereoselectivity. [ABSTRACT FROM AUTHOR]

The article focuses on a study of the primary steps of singlet fission in a pentacene film by using a combination of TG and 2D electronic spectroscopy. It mentions that coherent vibrational dynamics induces the ultrafast transition from the singlet excited electronic state to the triplet-pair state by a degeneracy of potential energy surfaces. It also mentions that significant vibronic coupling of the electronic wave packet to a few key intermolecular rocking modes in the low-frequency region.

Transition state structures are central to the rates and outcomes of chemical reactions, but their fleeting existence often leaves their properties to be inferred rather than observed. By treating polybutadiene with a difluorocarbene source, we embedded gem-difluorocyclopropanes (gDFCs) along the polymer backbone. We report that mechanochemical activation of the polymer under tension opens the gDFCs and traps a 1,3-diradical that is formally a transition state in their stress-free electrocyclic isomerization. The trapped diradical lives long enough that we can observe its noncanonical participation in bimolecular addition reactions. Furthermore, the application of a transient tensile force induces a net isomerization of the trans-gDK into its less-stable cis isomer, leading to the counterintuitive result that the gDK contracts in response to a transient force of extension. [ABSTRACT FROM AUTHOR]

Ion imaging reveals distinct photodissociation dynamics for propanal cations initially prepared in either the cis or gauche conformation, even though these isomers differ only slightly in energy and face a small interconversion barrier. The product kinetic energy distributions for the hydrogen atom elimination channels are bimodal, and the two peaks are readily assigned to propanoyl cation or hydroxyallyl cation coproducts. Ab initio multiple spawning dynamical calculations suggest that distinct ultrafast dynamics in the excited state deposit each conformer in isolated regions of the ground-state potential energy surface, and, from these distinct regions, conformer interconversion does not effectively compete with dissociation. [ABSTRACT FROM AUTHOR]

Electronic structure theory, i.e., quantum chemistry, is the fundamental building block for many problems in computational chemistry. We present a new distributed computing framework (BigChem), which allows for an efficient solution of many quantum chemistry problems in parallel. BigChem is designed to be easily composable and leverages industry-standard middleware (e.g., Celery, RabbitMQ, and Redis) for distributed approaches to large scale problems. BigChem can harness any collection of worker nodes, including ones on cloud providers (such as AWS or Azure), local clusters, or supercomputer centers (and any mixture of these). BigChem builds upon MolSSI packages, such as QCEngine to standardize the operation of numerous computational chemistry programs, demonstrated here with Psi4, xtb, geomeTRIC, and TeraChem. BigChem delivers full utilization of compute resources at scale, offers a programable canvas for designing sophisticated quantum chemistry workflows, and is fault tolerant to node failures and network disruptions. We demonstrate linear scalability of BigChem running computational chemistry workloads on up to 125 GPUs. Finally, we present ChemCloud, a web API to BigChem and successor to TeraChem Cloud. ChemCloud delivers scalable and secure access to BigChem over the Internet., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Even though cluster perturbation theory has been shown to be a robust noniterative alternative to coupled cluster theory, it is still plagued by high order polynomial computational scaling and the storage of higher order tensors. We present a proof-of-concept strategy for implementing a cluster perturbation theory ground-state energy series for the coupled cluster singles and doubles energy with N 4 computational scaling using tensor hypercontraction (THC). The reduction in computational scaling by two orders is achieved by decomposing two electron repulsion integrals, doubles amplitudes and multipliers, as well as selected double intermediates to the THC format. Using the outlined strategy, we showcase that the THC pilot implementations retain numerical accuracy to within 1 kcal/mol relative to corresponding conventional and density fitting implementations, and we empirically verify the N 4 scaling.

Directly imaging structural dynamics involving hydrogen atoms by ultrafast diffraction methods is complicated by their low scattering cross sections. Here we demonstrate that megaelectronvolt ultrafast electron diffraction is sufficiently sensitive to follow hydrogen dynamics in isolated molecules. In a study of the photodissociation of gas phase ammonia, we simultaneously observe signatures of the nuclear and corresponding electronic structure changes resulting from the dissociation dynamics in the time-dependent diffraction. Both assignments are confirmed by ab initio simulations of the photochemical dynamics and the resulting diffraction observable. While the temporal resolution of the experiment is insufficient to resolve the dissociation in time, our results represent an important step towards the observation of proton dynamics in real space and time.

Photomechanical crystals composed of three-dimensionally ordered and densely packed photochromes hold promise for high-performance photochemical actuators. However, bulk crystals with high structural ordering are severely limited in their flexibility, resulting in poor processibility and a tendency to fragment upon light exposure, while previous nano- or microcrystalline composites have lacked global alignment. Here we demonstrate a photon-fuelled macroscopic actuator consisting of diarylethene microcrystals in a polyethylene terephthalate host matrix. These microcrystals survive large deformations and show a high degree of three-dimensional ordering dictated by the anisotropic polyethylene terephthalate, which critically also has a similar stiffness. Overall, these ordered and compliant composites exhibit rapid response times, sustain a performance of over at least hundreds of cycles and generate work densities exceeding those of single crystals. Our composites represent the state-of-the-art for photochemical actuators and enable properties unattainable by single crystals, such as controllable, reversible and abrupt jumping (photosalient behaviour)., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

Light-driven molecular motors have generated considerable interest due to their potential applications in material and biological systems. Recently, Greb and Lehn reported a new class of molecular motors, chiral N -alkyl imines, which undergo unidirectional rotation induced by light and heat. The mechanism of unidirectional motion in molecular motors containing a C═N group has been assumed to consist of photoinduced torsion about the double bond. In this work, we present a computational study of the photoisomerization dynamics of a chiral N -alkyl imine motor. We find that the location and energetics of minimal energy conical intersections (MECIs) alone are insufficient to understand the mechanism of the motor. Furthermore, a key part of the mechanism consists of out-of-plane distortions of the N atom (followed by isomerization about the double bond). Dynamic effects and out-of-plane distortions are critical to understand the observed (rather low) quantum yield for photoisomerization. Our results provide hints as to how the photoisomerization quantum yield might be increased, improving the efficiency of this class of molecular motors.

Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromo-uracil in full dimensionality, i.e. , taking all the vibrational modes into consideration. The propagation starts on the resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.

Water is one of the most important, yet least understood, liquids in nature. Many anomalous properties of liquid water originate from its well-connected hydrogen bond network 1 , including unusually efficient vibrational energy redistribution and relaxation 2 . An accurate description of the ultrafast vibrational motion of water molecules is essential for understanding the nature of hydrogen bonds and many solution-phase chemical reactions. Most existing knowledge of vibrational relaxation in water is built upon ultrafast spectroscopy experiments 2-7 . However, these experiments cannot directly resolve the motion of the atomic positions and require difficult translation of spectral dynamics into hydrogen bond dynamics. Here, we measure the ultrafast structural response to the excitation of the OH stretching vibration in liquid water with femtosecond temporal and atomic spatial resolution using liquid ultrafast electron scattering. We observed a transient hydrogen bond contraction of roughly 0.04 Å on a timescale of 80 femtoseconds, followed by a thermalization on a timescale of approximately 1 picosecond. Molecular dynamics simulations reveal the need to treat the distribution of the shared proton in the hydrogen bond quantum mechanically to capture the structural dynamics on femtosecond timescales. Our experiment and simulations unveil the intermolecular character of the water vibration preceding the relaxation of the OH stretch., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

We present the first vacuum ultraviolet time-resolved photoelectron spectroscopy (VUV-TRPES) study of photoisomerization dynamics in the paradigmatic molecule cis -stilbene. A key reaction intermediate in its dynamics, known as the phantom state, has often been invoked but never directly detected in the gas phase. We report direct spectral signatures of the phantom state in isolated cis -stilbene, observed and characterized through a combination of VUV-TRPES and ab initio multiple spawning (AIMS) nonadiabatic dynamics simulations of the channel-resolved observable. The high VUV probe photon energy tracks the complete excited-state dynamics via multiple photoionization channels, from initial excitation to its return to the "hot" ground state. The TRPES was compared with AIMS simulations of the dynamics from initial excitation, to the phantom-state intermediate (an S 1 minimum), through to the ultimate electronic decay to the ground state. This combination revealed the unique spectral signatures and time-dependent dynamics of the phantom-state intermediate, permitting us to report here its direct observation.

We present the first investigation of excited state dynamics by resonant Auger-Meitner spectroscopy (also known as resonant Auger spectroscopy) using the nucleobase thymine as an example. Thymine is photoexcited in the UV and probed with X-ray photon energies at and below the oxygen K-edge. After initial photoexcitation to a ππ* excited state, thymine is known to undergo internal conversion to an nπ* excited state with a strong resonance at the oxygen K-edge, red-shifted from the ground state π* resonances of thymine (see our previous study Wolf, et al., Nat. Commun., 2017, 8, 29). We resolve and compare the Auger-Meitner electron spectra associated both with the excited state and ground state resonances, and distinguish participator and spectator decay contributions. Furthermore, we observe simultaneously with the decay of the nπ* state signatures the appearance of additional resonant Auger-Meitner contributions at photon energies between the nπ* state and the ground state resonances. We assign these contributions to population transfer from the nπ* state to a ππ* triplet state via intersystem crossing on the picosecond timescale based on simulations of the X-ray absorption spectra in the vibrationally hot triplet state. Moreover, we identify signatures from the initially excited ππ* singlet state which we have not observed in our previous study.

Woodward and Hoffman once jested that a very powerful Maxwell demon could seize a molecule of cyclobutene at its methylene groups and tear it open in a disrotatory fashion to obtain butadiene (Woodward, R. B.; Hoffmann, R. The Conservation of Orbital Symmetry. Angew. Chem., Int. Ed . 1969 , 8 , 781-853). Nearly 40 years later, that demon was discovered, and the field of covalent polymer mechanochemistry was born. In the decade since our demon was befriended, many fundamental investigations have been undertaken to build up our understanding of force-modified pathways for electrocyclic ring-opening reactions. Here, we seek to extend that fundamental understanding by exploring substituent effects in allowed and forbidden ring-opening reactions of cyclobutene (CBE) and benzocyclobutene (BCB) using a combination of single-molecule force spectroscopy (SMFS) and computation. We show that, while the forbidden ring-opening of cis -BCB occurs at a lower force than the allowed ring-opening of trans -BCB on the time scale of the SMFS experiment, the opposite is true for cis - and trans -CBE. Such a reactivity flip is explained through computational analysis and discussion of the so-called allowed/forbidden gap.

Direct irradiation of crystalline cis,cis -1,4-diphenyl-1,3-butadiene (cc-DPB) forms trans,trans -1,4-diphenyl-1,3,-butadiene via a concerted two-bond isomerization called the bicycle pedal (BP) mechanism. However, little is known about photoisomerization pathways in the solid state and there has been much debate surrounding the interpretation of volume-conserving isomerization mechanisms. The bicycle pedal photoisomerization is investigated using the quantum mechanics/molecular mechanics complete active space self-consistent field/Amber force-field method. Important details about how the steric environment influences isomerization mechanisms are revealed including how the one-bond flip and hula-twist mechanisms are suppressed by the crystal cavity, the nature of the seam space in steric environments, and the features of the bicycle pedal mechanism. Specifically, in the bicycle pedal, the phenyl rings of cc-DPB are locked in place and the intermolecular packing allows a passageway for rotation of the central diene in a volume-conserving manner. In contrast, the bicycle pedal rotation in the gas phase is not a stable pathway, so single-bond rotation mechanisms become operative instead. Furthermore, the crystal BP mechanism is an activated process that occurs completely on the excited state; the photoproduct can decay to the ground state through radiative and non-radiative pathways. The present models, however, do not capture the quantitative activation barriers, and more work is needed to better model reactions in crystals. Last, the reaction barriers of the different crystalline conformations within the unit cell of cc-DPB are compared to investigate the possibility for conformation-dependent isomerization. Although some difference in reaction barriers is observed, the difference is most likely not responsible for the experimentally observed periods of fast and slow conversion.

The dynamics of proton transfer to the aprotic solvent 1-methylimidazole (MeIm, proton acceptor) from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid trisodium salt (HPTS) was investigated using fast fluorescence measurements. The closely related molecule, 8-methoxypyrene-1,3,6-trisulfonic acid trisodium salt (MPTS), which is not a photoacid, was also studied for comparison. Following optical excitation, the wavelength-dependent population dynamics of HPTS in MeIm resulting from the deprotonation process were collected over the entire fluorescence emission window. Analysis of the time-dependent fluorescence spectra revealed four distinct fluorescence bands that appear and decay on different time scales. We label these four states as protonated (P), associated I (A I ), associated II (A II ), and deprotonated (D). We find that the simple kinetic scheme of P → A I → A II → D is not consistent with the data. Instead, the kinetic scheme that describes the data has P decaying into A I , which mainly goes on to deprotonation (D), but A I can also feed into A II . A II can return to A I or decay to the ground state, but does not deprotonate within experimental error. Quantum chemistry and excited state QM/MM Born-Oppenheimer molecular dynamics simulations indicate that A I and A II are two H-bonding conformations of MeIm to the HPTS hydroxyl, axial, and equatorial, respectively.

A multilayer multi-configurational theory framework is adapted to simulate circuit-based quantum computers. Quantum addition of superpositions of an exponential number of summands is performed in polynomial time with high accuracy. We demonstrate numerically accurate calculations including up to one million qubits for entangling benchmarks. Simulation cost can be assessed by entropy-based entanglement measures. For the considered systems, we show that the entanglement only grows weakly with the system size. The present simulations demonstrate how quantum algorithms in low-entropy regimes can be used efficiently on classically simulated quantum computers.