Search

Modeling chemical reactions and complicated molecular systems has been proposed as the "killer application" of a future quantum computer. Accurate calculations of derivatives of molecular eigenenergies are essential toward this end, allowing for geometry optimization, transition state searches, predictions of the response to an applied electric or magnetic field, and molecular dynamics simulations. In this work, we survey methods to calculate energy derivatives, and present two new methods: one based on quantum phase estimation, the other on a low-order response approximation. We calculate asymptotic error bounds and approximate computational scalings for the methods presented. Implementing these methods, we perform geometry optimization on an experimental quantum processor, estimating the equilibrium bond length of the dihydrogen molecule to within 0.014 Å of the full configuration interaction value. Within the same experiment, we estimate the polarizability of the H 2 molecule, finding agreement at the equilibrium bond length to within 0.06 a.u. ( 2 % relative error). [ABSTRACT FROM AUTHOR]

A density functional theory is used to determine the solvent effect on nuclear magnetic shielding parameters by using example of sensitive nitrogen shift of acetonitrile immersed in a selected set of solvents, like water, chloroform, and cyclohexane. The structural errors introduced can be largely avoided by employing the more expensive and accurate Car-Parrinello molecular dynamics simulations method to generate the ensemble of configurations.

The studies of the wavelength-resolved multiple photon photodissociation of discrete, gas-phase uranyl complexes containing a single anionic ligand with or without ligated solvent molecule is presented. The comparison between the uranyl frequencies and [[U[O.sub.2]AS].sup.+] species with different solvent neutrals showed the values to decrease with increase in neutral nucleophilicity.

The Resolution of Identity (RI) technique has been employed to speed up the use of hybrid exchange-correlation (xc) functionals at the TDDFT level using the Hybrid Diagonal Approximation. The RI has been implemented within the polTDDFT algorithm (a complex damped polarization method) in the AMS/ADF suite of programs. A speedup factor of 30 has been obtained with respect to a previous numerical implementation, albeit with the same level of accuracy. Comparison of TDDFT simulations with the experimental photoabsorption spectra of the cluster series Au 8 n +4 (SR) 4 n +8 ( n = 3-6; R = C 6 H 5 ) showed the excellent accuracy and efficiency of the method. Results were compared with those obtained via the more simplified and computationally cheaper TDDFT+TB and sTDDFT methods. The present method represents an accurate as well as computationally affordable approach to predict photoabsorption spectra of complex species, realizing an optimal compromise between accuracy and computational efficiency, and is suitable for applications to large metal clusters with sizes up to several hundreds of atoms.

Low-order scaling GW implementations for molecules are usually restricted to approximations with diagonal self-energy. Here, we present an all-electron implementation of quasiparticle self-consistent GW for molecular systems. We use an efficient algorithm for the evaluation of the self-energy in imaginary time, from which a static non-local exchange-correlation potential is calculated via analytical continuation. By using a direct inversion of iterative subspace method, fast and stable convergence is achieved for almost all molecules in the GW100 database. Exceptions are systems which are associated with a breakdown of the single quasiparticle picture in the valence region. The implementation is proven to be starting point independent and good agreement of QP energies with other codes is observed. We demonstrate the computational efficiency of the new implementation by calculating the quasiparticle spectrum of a DNA oligomer with 1,220 electrons using a basis of 6,300 atomic orbitals in less than 4 days on a single compute node with 16 cores. We use then our implementation to study the dependence of quasiparticle energies of DNA oligomers consisting of adenine-thymine pairs on the oligomer size. The first ionization potential in vacuum decreases by nearly 1 electron volt and the electron affinity increases by 0.4 eV going from the smallest to the largest considered oligomer. This shows that the DNA environment stabilizes the hole/electron resulting from photoexcitation/photoattachment. Upon inclusion of the aqueous environment via a polarizable continuum model, the differences between the ionization potentials reduce to 130 meV, demonstrating that the solvent effectively compensates for the stabilizing effect of the DNA environment. The electron affinities of the different oligomers are almost identical in the aqueous environment., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Förster and Visscher.)

Light-harvesting complex II (LHCII) is a pigment-protein complex present in higher plants and green algae. LHCII represents the main site of light absorption, and its role is to transfer the excitation energy toward the photosynthetic reaction centers, where primary energy conversion reactions take place. The optical properties of LHCII are known to depend on protein conformation. However, the relation between the structural and spectroscopic properties of the pigments is not fully understood yet. In this respect, previous classical molecular dynamics simulations of LHCII in a model membrane [ Sci. Rep. 2015 , 5 , 1-10] have shown that the configuration and excitonic coupling of a chlorophyll (Chl) dimer functioning as the main terminal emitter of the complex are particularly sensitive to conformational changes. Here, we use quantum chemistry calculations to investigate in greater detail the effect of pigment-pigment interactions on the excited-state landscape. While most previous studies have used a local picture in which electrons are localized on single pigments, here we achieve a more accurate description of the Chl dimer by adopting a supramolecular picture where time-dependent density functional theory is applied to the whole system at once. Our results show that specific dimer configurations characterized by shorter inter-pigment distances can result in a sizable intensity decrease (up to 36%) of the Chl absorption bands in the visible spectral region. Such a decrease can be predicted only when accounting for Chl-Chl charge-transfer excitations, which is possible using the above-mentioned supramolecular approach. The charge-transfer character of the excitations is quantified by two types of analyses: one focusing on the composition of the excitations and the other directly on the observable total absorption intensities.

Argon is often considered as an innocent probe that can be attached and detached to study the structure of a particular species without perturbing the species too much. We have investigated whether this assumption also holds for small copper cationic clusters and demonstrated that small but significant charge transfer from argon to metal changes the remaining binding positions, leading in general, to weaker binding of other argon atoms. The exception is binding to just one copper ion, where the binding of the first argon facilitates the binding of the second.

Vibrational circular dichroism (VCD) is one of the major spectroscopic tools to study peptides. Nevertheless, a full understanding of what determines the signs and intensities of VCD bands of these compounds in the amide I and amide II spectral regions is still far from complete. In the present work, we study the origin of these VCD signals using the general coupled oscillator (GCO) analysis, a novel approach that has recently been developed. We apply this approach to the ForValNHMe model peptide in both α-helix and β-sheet configurations. We show that the intense VCD signals observed in the amide I and amide II spectral regions essentially have the same underlying mechanism, namely, the through-space coupling of electric dipoles. The crucial role played by intramolecular hydrogen bonds in determining VCD intensities is also illustrated. Moreover, we find that the contributions to the rotational strengths, considered to be insignificant in standard VCD models, may have sizable magnitudes and can thus not always be neglected. In addition, the VCD robustness of the amide I and II modes has been investigated by monitoring the variation of the rotational strength and its contributing terms during linear transit scans and by performing calculations with different computational parameters. From these studies-and in particular, the decomposition of the rotational strength made possible by the GCO analysis-it becomes clear that one should be cautious when employing measures of robustness as proposed previously.

Cationic cobalt clusters complexed with water Con(+)-H2O (n = 6-20) are produced through laser ablation and investigated via infrared multiple photon dissociation (IR-MPD) spectroscopy in the 200-1700 cm(-1) spectral range. All spectra exhibit a resonance close to the 1595 cm(-1) frequency of the free water bending vibration, indicating that the water molecule remains intact upon adsorption. For n = 6, the frequency of this band is blue-shifted, but it gradually converges to the free water value with increasing cluster size. In the lower-frequency range (200-650 cm(-1)) the spectra contain several bands which show a very regular frequency evolution, suggesting that the exact cluster geometry has little effect on the water-surface interaction. Density functional theory (DFT) calculations are carried out at the OPBE/TZP level for three representative sizes (n = 6, 9, 13) and indicate that the vibrations responsible for the resonances correspond to bending and torsional modes between the cluster and water moieties. The potential energy surfaces describing these interactions are very shallow, making the calculated harmonic frequencies and IR intensities very sensitive to small geometrical perturbations. We conclude that harmonic frequency calculations on (local) minima structures provide insufficient information for these types of cluster complexes and need to be complemented with calculations that provide a more extensive sampling of the potential energy surface.

Four-component relativistic time-dependent density functional theory (TD-DFT) is used to study charge-transfer (CT) excitation energies of the uranyl molecule as well as the uranyl tetrachloride complex. Adiabatic excitation energies and vibrational frequencies of the excited states are calculated for the lower energy range of the spectrum. The results for TD-DFT with the CAM-B3LYP exchange-correlation functional for the [UO(2)Cl(4)](2-) system are in good agreement with the experimentally observed spectrum of this species and agree also rather well with other theoretical data. Use of the global hybrid B3LYP gives qualitatively correct results, while use of the BLYP functional yields results that are qualitatively wrong due to the too low CT states calculated with this functional. The applicability of the overlap diagnostic of Peach et al. (J. Chem. Phys.2008, 128, 044118) to identify such CT excitations is investigated for a wide range of vertical transitions using results obtained with three different approximate exchange-correlation functionals: BLYP, B3LYP, and CAM-B3LYP.

Isoflavones are known to show fluorescence with intensities that depend strongly on the solvent properties and exhibit Stokes' shifts as large as 1.4 eV. While some of this behavior can be explained by (excited state) deprotonation, this mechanism does not apply for all isoflavones. The aim of this study is to computationally and experimentally investigate the reasons for this anomalous behavior of neutral isoflavones, taking the daidzein molecule as a model compound. We find that the absence in fluorescence in aprotic solvents and the weak fluorescence in protic solvents can be explained by a change of order of the lowest singlet states in which a fluorescent charge-transfer state lies below the nonfluorescent locally excited state in water but not in acetonitrile. The large Stokes' shift is partly due to a significant rotation among the chromone-phenyl bond in the excited state.

The Free-Electron Laser for Infrared Experiments (FELIX) was used to study the wavelength-resolved multiple photon photodissociation of discrete, gas-phase uranyl (UO22+) complexes containing a single anionic ligand (A), with or without ligated solvent molecules (S). The uranyl antisymmetric and symmetric stretching frequencies were measured for complexes with general formula [UO2A(S)n]+, where A was hydroxide, methoxide, or acetate; S was water, ammonia, acetone, or acetonitrile; and n = 0-3. The values for the antisymmetric stretching frequency for uranyl ligated with only an anion ([UO2A]+) were as low or lower than measurements for [UO2]2+ ligated with as many as five strong neutral donor ligands and are comparable to solution-phase values. This result was surprising because initial DFT calculations predicted values that were 30-40 cm(-1) higher, consistent with intuition but not with the data. Modification of the basis sets and use of alternative functionals improved computational accuracy for the methoxide and acetate complexes, but calculated values for the hydroxide were greater than the measurement regardless of the computational method used. Attachment of a neutral donor ligand S to [UO2A]+ produced [UO2AS]+, which produced only very modest changes to the uranyl antisymmetric stretch frequency, and did not universally shift the frequency to lower values. DFT calculations for [UO2AS]+ were in accord with trends in the data and showed that attachment of the solvent was accommodated by weakening of the U-anion bond as well as the uranyl. When uranyl frequencies were compared for [UO2AS]+ species having different solvent neutrals, values decreased with increasing neutral nucleophilicity.