Search

Environmental effects in excitation energy transfer have mostly been modeled by baths of harmonic oscillators, but to what extent such modeling provides a reliable description of actual interactions between molecular systems and environments remains an open issue. We address this issue by investigating fluctuations in the excitation energies of the light harvesting 2 complex using a realistic all-atomistic simulation of the potential energy surface. Our analyses reveal that molecular motions exhibit significant anharmonic features, even for underdamped intramolecular vibrations. In particular, we find that the anharmonicity contributes to the broadening of spectral densities and substantial overlaps between neighboring peaks, which complicates the meaning of mode frequencies constituting a bath model. Thus, we develop a strategy to construct a minimally underdamped harmonic bath that has a clear connection to all-atomistic dynamics by utilizing actual normal modes of molecules but optimizing their frequencies such that the resulting bath model can best reproduce the all-atomistic simulation results. By subtracting the underdamped contribution from the entire fluctuations, we also show that identifying a residual spectral density representing all other contributions with overdamped behavior is possible. We find that this can be fitted well with a well-established analytic form of a spectral density function or, alternatively, modeled as explicit time dependent fluctuations with muti-exponential or power law type correlation functions. We provide an assessment and the implications of these possibilities. The approach presented here can also serve as a general strategy to construct a simplified bath model that can effectively represent the underlying all-atomistic bath dynamics. [ABSTRACT FROM AUTHOR]

This document is an erratum for a previous paper titled "Modified Fermi's golden rule rate expressions." The erratum corrects an error in equation (41) by adding a missing factor of 1/2. The correction does not change any major conclusions or implications of the paper, but it does require the correction of figures 1-4, which show results of model calculations. The corrected figures and their captions are provided in the erratum. [Extracted from the article]

Peptides are promising therapeutic agents for various biological targets due to their high efficacy and low toxicity, and the design of peptide ligands with high binding affinity to the target of interest is of utmost importance in peptide‐based drug design. Introducing a conformational constraint to a flexible peptide ligand using a side‐chain lactam‐bridge is a convenient and efficient method to improve its binding affinity to the target. However, in general, such a small structural modification to a flexible ligand made with the intent of lowering the configurational entropic penalty for binding may have unintended consequences in different components of the binding enthalpy and entropy, including the configurational entropy component, which are still not clearly understood. Toward probing this, we examine different components of the binding enthalpy and entropy as well as the underlying structure and dynamics, for a side‐chain lactam‐bridged peptide inhibitor and its flexible analog forming complexes with vascular endothelial growth factor (VEGF), using all‐atom molecular dynamics simulations. It is found that introducing a side‐chain lactam‐bridge constraint into the flexible peptide analog led to a gain in configurational entropy change but losses in solvation entropy, solute internal energy, and solvation energy changes upon binding, pinpointing the opportunities and challenges in drug design. The present study features an interplay between configurational and solvation entropy changes, as well as the one between binding enthalpy and entropy, in ligand‐target binding upon imposing a conformational constraint into a flexible ligand. [ABSTRACT FROM AUTHOR]

A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.

Vascular endothelial growth factor 165 (VEGF165) is a prominent isoform of the VEGF-A protein that plays a crucial role in various angiogenesis-related diseases. It is homodimeric, and each of its monomers is composed of two domains connected by a flexible linker. DNA aptamers, which have emerged as potent therapeutic molecules for many proteins with high specificity and affinity, can also work for VEGF165. A DNA aptamer heterodimer composed of monomers of V7t1 and del5-1 connected by a flexible linker (V7t1:del5-1) exhibits a greater binding affinity with VEGF165 compared to either of the two monomers alone. Although the structure of the complex formed between the aptamer heterodimer and VEGF165 is unknown due to the highly flexible linkers, gaining structural information will still be valuable for future developments. Toward this end of accessing structural information, we adopt an ensemble docking approach here. We first obtain an ensemble of structures for both VEGF165 and the aptamer heterodimer by considering both small- and large-scale motions. We then proceed through an extraction process based on ensemble docking, molecular dynamics simulations, and binding free energy calculations to predict the structures of the VEGF165/V7t1:del5-1 complex. Through the same procedures, we reach a new aptamer heterodimer that bears a locked nucleic acid-modified counterpart of V7t1, namely RNV66:del5-1, which also binds well with VEGF165. We apply the same protocol to the monomeric units V7t1, RNV66, and del5-1 to target VEGF165. We observe that V7t1:del5-1 and RNV66:del5-1 show higher binding affinities with VEGF165 than any of the monomers, consistent with experiments that support the notion that aptamer heterodimers are more effective anti-VEGF165 aptamers than monomeric aptamers. Among the five different aptamers studied here, the newly designed RNV66:del5-1 shows the highest binding affinity with VEGF165. We expect that our ensemble docking approach can help in de novo designs of homo/heterodimeric anti-angiogenic drugs to target the homodimeric VEGF165. [ABSTRACT FROM AUTHOR]

The realization of quantum advantage with noisy-intermediate-scale quantum (NISQ) machines has become one of the major challenges in computational sciences. Maintaining coherence of a physical system with more than ten qubits is a critical challenge that motivates research on compact system representations to reduce algorithm complexity. Toward this end, the variational quantum eigensolver (VQE) used to perform quantum simulations is considered to be one of the most promising algorithms for quantum chemistry in the NISQ era. We investigate reduced mapping of one spatial orbital to a single qubit to analyze the ground state energy in a way that the Pauli operators of qubits are mapped to the creation/annihilation of singlet pairs of electrons. To include the effect of non-bosonic (or non-paired) excitations, we introduce a simple correction scheme in the electron correlation model approximated by the geometrical mean of the bosonic (or paired) terms. Employing it in a VQE algorithm, we assess ground state energies of H2O, N2, and Li2O in good agreement with full configuration interaction (FCI) models respectively, using only 6, 8, and 12 qubits with quantum gate depths proportional to the squares of the qubit counts. With the adopted seniority-zero approximation that uses only one half of the qubit counts of a conventional VQE algorithm, we find that our non-bosonic correction method reaches reliable quantum chemistry simulations at least for the tested systems. [ABSTRACT FROM AUTHOR]

Understanding nonradiative transition mechanisms is important in various situations. However, compared with radiative processes where temporal profiles of photon emission can be monitored in a straightforward manner, experimentally accessing the rate information may not be an easy task with nonradiative transitions. Hence, applying theoretical tools toward predicting the rates can be a useful tactic. Such predictions become very useful in designing optoelectronic materials as in the molecules adopted for constructing organic light‐emitting diodes (OLEDs). The correlation function formalism is a method that can fulfill the purpose of designing OLED materials. The formalism requires information regarding the vibrational normal modes of the two electronic states before and after the transition. Because the method is also based on harmonic oscillator approximation, it can actually fail to provide high reliability when there is a large geometric distortion between the initial and the final states. In fact, the harmonic normal mode picture is more prone to lose reliability in the Cartesian coordinates than in the internal ones even at a small distortion, and hence adopting internal coordinates may be more preferable for practical calculations. This is because normal mode mixing becomes less severe when molecular coordinates are described with internal degrees of freedom such as bond stretching, bending, and torsion. In this regard, how much more reliable the nonradiative rate predictions in OLED materials become with the use of internal coordinate system deserves a close inspection. In this account, we review on the derivation of the correlation function formalism and provide how it can be adapted toward the use of the internal coordinates. As a demonstration, we evaluate the intersystem crossing and the internal conversion rates of a series of thermally activated delayed fluorescence (TADF) molecules with both Cartesian and internal coordinate systems. Overall, handling transitions involving substantial structural changes is improved indeed with the internal coordinates. However, limitations are still apparent for the TADF systems with a flexible donor–acceptor type construct especially when the inevitable inter‐domain twisting takes place with the electronic transition. Future prospect for handling the issue is commented as a concluding remark. [ABSTRACT FROM AUTHOR]

Expressions for analytical molecular gradients of core-excited states have been derived and implemented for the hierarchy of algebraic diagrammatic construction (ADC) methods up to extended second-order within the core–valence separation (CVS) approximation. We illustrate the use of CVS-ADC gradients by determining relaxed core-excited state potential energy surfaces and optimized geometries for water, formic acid, and benzene. For water, our results show that in the dissociative lowest core-excited state, a linear configuration is preferred. For formic acid, we find that the O K-edge lowest core-excited state is non-planar, a fact that is not captured by the equivalent core approximation where the core-excited atom with its hole is replaced by the "Z + 1" neighboring atom in the periodic table. For benzene, the core-excited state gradients are presented along the Jahn–Teller distorted geometry of the 1s → π* excited state. Our development may pave a new path to studying the dynamics of molecules in their core-excited states. [ABSTRACT FROM AUTHOR]

Mapping basis solutions provide efficient ways for simulating mixed quantum–classical (MQC) dynamics in complex systems by matching multiple quantum states of interest to some fictitious physical states. Recently, various MQC methods were devised such that two harmonic oscillators are employed to represent each electronic state, showing improvements over one-oscillator-based methods. Here, we introduce and analyze newly modified mapping approximations of the quantum–classical Liouville equation (QCLE) using two oscillators for each electronic state. We design two separate mapping relations that we can adopt toward simulating dynamics and computing expectation values. Through the process, two MQC methods can be constructed, one of which actually reproduces the population dynamics of the forward and backward trajectory solution of QCLE. By applying the methods to spin–boson systems with a range of parameters, we find out that the choice of mapping relations greatly affects the simulation results. We also show that further improvement is possible through using modified identity operator formulations. Our findings may be helpful in constructing improved MQC methods in the future. [ABSTRACT FROM AUTHOR]

In open quantum system dynamics, rich information about the major energy relaxation channels and corresponding relaxation rates can be elucidated by monitoring the vibrational energy flow among individual bath modes. However, such calculations often become tremendously difficult as the complexity of the subsystem–bath coupling increases. In this paper, we attempt to make this task feasible by using a mixed quantum–classical method, the Poisson-bracket mapping equation with non-Hamiltonian modification (PBME-nH) [H. W. Kim and Y. M. Rhee, J. Chem. Phys. 140, 184106 (2014)]. For a quantum subsystem bilinearly coupled to harmonic bath modes, we derive an expression for the mode energy in terms of the classical positions and momenta of the nuclei, while keeping consistency with the energy of the quantum subsystem. The accuracy of the resulting expression is then benchmarked against a numerically exact method by using relatively simple models. Although our expression predicts a qualitatively correct dissipation rate for a range of situations, cases involving a strong vibronic resonance are quite challenging. This is attributed to the inherent lack of quantum back reaction in PBME-nH, which becomes significant when the subsystem strongly interacts with a small number of bath modes. A rigorous treatment of such an effect will be crucial for developing quantitative simulation methods that can handle generic subsystem–bath coupling. [ABSTRACT FROM AUTHOR]

Cohosts based on hole transporting and electron transporting materials often act as exciplexes in the form of intermolecular charge transfer complexes. Indeed, exciplex-forming cohosts have been widely developed as the host materials for efficient phosphorescent organic light-emitting diodes (OLEDs). In host–guest systems of OLEDs, the guest can be excited by two competing mechanisms, namely, excitation energy transfer (EET) and charge transfer (CT). Experimentally, it has been reported that the EET mechanism is dominant and the excitons are primarily formed in the host first and then transferred to the guest in phosphorescent OLEDs based on exciplex-forming cohosts. With this, exciplex-forming cohosts are widely employed for avoiding the formation of trapped charge carriers in the phosphorescent guest. However, theoretical studies are still lacking toward elucidating the relative importance between EET and CT processes in exciting the guest molecules in such systems. Here, we obtain the kinetics of guest excitation processes in a few trimer model systems consisting of an exciplex-forming cohost pair and a phosphorescent guest. We adopt the Förster resonance energy transfer (FRET) rate constants for the electronic transitions between excited states toward solving kinetic master equations. The input parameters for calculating the FRET rate constants are obtained from density functional theory (DFT) and time-dependent DFT. The results show that while the EET mechanism is important, the CT mechanism may still play a significant role in guest excitations. In fact, the relative importance of CT over EET depends strongly on the location of the guest molecule relative to the cohost pair. This is understandable as both the coupling for EET and the interaction energy for CT are strongly influenced by the geometric constraints. Understanding the energy transfer pathways from the exciplex state of cohost to the emissive state of guest may provide insights for improving exciplex-forming materials adopted in OLEDs. [ABSTRACT FROM AUTHOR]