Kestutis Aidas

The Dalton quantum chemistry program system

Dalton is a powerful general‐purpose program system for the study of molecular electronic structure at the Hartree–Fock, Kohn–Sham, multiconfigurational self‐consistent‐field, Møller–Plesset, configuration‐interaction, and coupled‐cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic‐structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge‐origin‐invariant manner. Frequency‐dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one‐, two‐, and three‐photon processes. Environmental effects may be included using various dielectric‐medium and quantum‐mechanics/molecular‐mechanics models. Large molecules may be studied using linear‐scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

On the performance of quantum chemical methods to predict solvatochromic effects: The case of acrolein in aqueous solution

The performance of the Hartree-Fock method and the three density functionals B3LYP, PBE0, and CAM-B3LYP is compared to results based on the coupled cluster singles and doubles model in predictions of the solvatochromic effects on the vertical n-->pi* and pi-->pi* electronic excitation energies of acrolein. All electronic structure methods employed the same solvent model, which is based on the combined quantum mechanics/molecular mechanics approach together with a dynamical averaging scheme. In addition to the predicted solvatochromic effects, we have also performed spectroscopic UV measurements of acrolein in vapor phase and aqueous solution. The gas-to-aqueous solution shift of the n-->pi* excitation energy is well reproduced by using all density functional methods considered. However, the B3LYP and PBE0 functionals completely fail to describe the pi-->pi* electronic transition in solution, whereas the recent CAM-B3LYP functional performs well also in this case. The pi-->pi* excitation energy of acrolein in water solution is found to be very dependent on intermolecular induction and nonelectrostatic interactions. The computed excitation energies of acrolein in vacuum and solution compare well to experimental data.

On the Accuracy of Density Functional Theory to Predict Shifts in Nuclear Magnetic Resonance Shielding Constants due to Hydrogen Bonding.

We present the first systematic investigation of shifts in the nuclear magnetic resonance (NMR) shielding constant due to hydrogen bonding using either the series of wave function based methods, Hartree-Fock (HF), second-order Møller-Plesset perturbation theory (MP2), Coupled Cluster Singles and Doubles (CCSD) and CCSD extended with an approximate description of triples (CCSD(T)), or Density Functional Theory (DFT) employing either the B3LYP, PBE0, or KT3 exchange correlation (xc) functionals. The molecular systems considered are (i) the water dimer and (ii) formaldehyde in complex with two water molecules. Specially for the (17)O in formaldehyde we observe significant differences between the DFT and CCSD(T) predictions. However, the extent of these deviations depends crucially on the applied xc functional. Compared to CCSD(T) we find the KT3 functional to provide accurate results, whereas both B3LYP and PBE0 are in significant error. Potential consequences of this observation are discussed in the context of general predictions of NMR shielding constants in condensed phase.

Modeling the structure and absorption spectra of stilbazolium merocyanine in polar and nonpolar solvents using hybrid QM/MM techniques.

We have performed Car-Parrinello mixed quantum mechanics/molecular mechanics (CP-QM/MM) calculations for stilbazolium merocyanine (SM) in polar and nonpolar solvents in order to explore the role of solute molecular geometry, solvation shell structure, and different interaction mechanisms on the absorption spectra and its dependence on solvent polarity. On the basis of the average bond length values and group charge distributions, we find that the SM molecule remains in a neutral quinonoid form in chloroform (a nonpolar solvent) while it transforms to a charge-separated benzenoid form in water (a polar solvent). Based on a quantum mechanical/molecular mechanical response technique, with different MM descriptions for the water environment, absorption spectra were obtained as averages over configurations derived from the CP-QM/MM simulations. We show that for SM in water the solute polarization plays a major role in predictions of the λ(max) and solvatochromic shift and that once this effect is included the contributions from solvent polarization and intermolecular charge transfer become less important. For SM in chloroform and water solvents, we have also performed absorption spectra calculations using a polarizable continuum model in order to address its relative performance compared to the QM/MM response technique. In the case of SM in water, our study supports the notion that, in order to predict accurate absorption spectra and solvatochromic shifts, it is important to use a discrete description of the solvent when it, as in water, is involved in site-specific interaction with the solute molecule. The technique is thus shown to outperform the more conventional polarizable continuum model in predicting the solvatochromic shift.

NMR and Quantum Chemistry Study of Mesoscopic Effects in Ionic Liquids

(1)H, (13)C, and (81)Br NMR spectra of the neat room-temperature ionic liquid (RTIL), namely, 1-decyl-3-methyl-imidazolium bromide ([C(10)mim][Br]) as well as its solutions in acetonitrile, dichloromethane, methanol, and water have been investigated. The most important observation of the present work is the significant broadening of (81)Br NMR signal in the solutions of [C(10)mim][Br] in organic solvents, which molecules tend to associate into hydrogen bond networks and the appearance of the complex contour of (81)Br NMR signal in the neat RTIL as well as in the liquid crystalline (LC) ionogel formed in RTIL/water solution. The complex structure of (81)Br signal changes upon heating and dilution in water. It disappears at ca. 353 K and in the aqueous solution below ca. 0.1 mol fraction of RTIL. Several new (1)H NMR signals appear at the [C(10)mim][Br]/water compositions just before the solidification of the sample (approximately 0.3 mol fraction of [C(10)mim][Br]). These additional peaks can be attributed to the H(2)O protons placed in inhomogeneous regions of the sample or due to the appearance of nonequivalent water sites in LC ionogel, the exchange between which is highly restricted or even frozen. The complex shape of (81)Br NMR signal can originate from the presence of supra-molecular structures (mesoscopic domains) that live over the period of the NMR time-scale due to a very high viscosity of [C(10)mim][Br]. These domains exhibit some features of partially disordered solids (liquid- or plastic crystals). To evaluate the static and dynamic contributions into the relaxation rate of (81)Br nuclei, the quantum chemistry calculations of the electronic structure, magnetic shielding, and electric field gradient (EFG) tensors of [C(10)mim][Br] and related model systems (Br(-).6H(2)O cluster, with addition of the dipoles (hydrogen fluoride) and charged particles - cations: H(+) or C(1)mim(+)) were performed.

Prediction of spin-spin coupling constants in solution based on combined density functional theory/molecular mechanics

We present theory and implementation of calculation of spin-spin coupling constants within combined quantum mechanics/molecular mechanics methods. Special attention is given to the role of explicit solvent polarization as well as the molecular consequences due to hydrogen bonding. The model is generally applicable but is here implemented for the case of density functional theory. First applications to liquid water and acetylene in aqueous solution are presented. Good agreement between theory and experiment is obtained in both cases, thereby showing the strength of our approach. Finally, spin-spin coupling constants across hydrogen bonds are discussed considering for the first time the role of an explicit solvent on this class of spin-spin couplings.

Amyloid Fibril-Induced Structural and Spectral Modifications in the Thioflavin-T Optical Probe

Motivated by future possibilities to design target molecules for fibrils with diagnostic or therapeutic capability related to amyloidosis diseases, we investigate in this work the dielectric nature of amyloid fibril microenvironments in different binding sites using an optical probe, thioflavin-T (THT), which has been used extensively to stain such fibrils. We study the fibril-environment-induced structural and spectral changes of THT at each binding site and compare the results to the fibril-free situation in aqueous solution. All binding sites are found to show a similar effect with respect to the conformational changes of THT; in the presence of the fibril, its molecular geometry tends to become planarized. However, depending on the dielectric nature of the specific binding site, a red shift, blue shift, or no shift in the absorption spectra of THT is predicted. Interestingly, the experimentally measured red shift in the spectra is seen only when THT binds to one of the core or surface-binding sites. It is found that the dielectric nature of the microenvironment in the fibril is strongly nonhomogeneous.

Computational protocols for prediction of solute NMR relative chemical shifts. A case study of L‐tryptophan in aqueous solution

In this study, we have applied two different spanning protocols for obtaining the molecular conformations of L‐tryptophan in aqueous solution, namely a molecular dynamics simulation and a molecular mechanics conformational search with subsequent geometry re‐optimization of the stable conformers using a quantum mechanically based method. These spanning protocols represent standard ways of obtaining a set of conformations on which NMR calculations may be performed. The results stemming from the solute–solvent configurations extracted from the MD simulation at 300 K are found to be inferior to the results stemming from the conformations extracted from the MM conformational search in terms of replicating an experimental reference as well as in achieving the correct sequence of the NMR relative chemical shifts of L‐tryptophan in aqueous solution. We find this to be due to missing conformations visited during the molecular dynamics run as well as inaccuracies in geometrical parameters generated from the classical molecular dynamics simulations.

Proton transfer in hydrogen-bonded pyridine/acid systems: the role of higher aggregation

In this work the role of higher molecular aggregation in the proton transfer processes within hydrogen bond (H-bond) is investigated. The H-bonded complex consisting of 4-cyanopyridine (CyPy) with trichloroacetic acid (TCA) has been studied in the solutions of acetonitrile, carbon tetrachloride, chloroform and dichloroethane as solvent by FTIR spectroscopy and quantum chemical DFT calculations. In order to illustrate the effect of increasing H-bond strength FTIR investigations have also been performed on solutions of CyPy with H(2)O, acetic-, trifluoroacetic- and methanesulfonic acids. Proton states in the H-bond have been monitored using vibrational CyPy ring modes in FTIR spectra. The stabilization of the CyPy/TCA complex in its protonated form upon increasing polarity of the solvent has been evidenced. It was shown that formation of the CyPy/(TCA)(2) aggregates in the solutions favors the proton transfer process. An X-ray diffraction study has been performed on a single 1 : 2 co-crystal of pyridine/3,5-dinitrobenzoic acid. The H-bond motif found in this system exhibits the same connectivity by strong hydrogen bonds N-H(+)[dot dot dot]O(-) and O-H[dot dot dot]O as that in the CyPy/(TCA)(2) complex predicted by DFT calculation. Certain discrepancies are observed in C-H[dot dot dot]O connectivity only. The networks of H-bonds in both assemblies differ from those usually pictured for 1 : 2 base/carboxylic acid complexes in the literature.

NMR and Raman Spectroscopy Monitoring of Proton/Deuteron Exchange in Aqueous Solutions of Ionic Liquids Forming Hydrogen Bond: A Role of Anions, Self-Aggregation, and Mesophase Formation

The H/D exchange process in the imidazolium-based room temperature ionic liquids (RTILs) 1-decyl-3-methyl-imidazolium bromide- and chloride ([C10mim][Br] and [C10mim][Cl]) in D2O solutions of various concentrations was studied applying (1)H, (13)C NMR, and Raman spectroscopy. The time dependencies of integral intensities in NMR spectra indicate that the H/D exchange in [C10mim][Br] at very high dilution (10(-4) mole fraction of RTIL) runs only slightly faster than in [C10mim][Cl]. The kinetics of this process drastically changes above critical aggregation concentration (CAC). The time required to reach the apparent reaction saturation regime in the solutions of 0.01 mole fraction of RTIL was less 10 h for [C10mim][Br], whereas no such features were seen for [C10mim][Cl] even tens of days after the sample was prepared. The H/D exchange was not observed in the liquid crystalline gel mesophase. The role of anions, self-aggregation (micellization), and mesophase formation has been discussed. Crucial influence of Br(-) and Cl(-) anions on the H/D exchange rates above CAC could be related to the short-range ordering and molecular microdynamics, in particular that of water molecules. The concept of the conformational changes coupled with the H/D exchange in imidazolium-based ionic liquids with longer hydrocarbon chains can be rejected in the light of (13)C NMR experiment. The revealed changes in (13)C NMR spectra are caused by the secondary ((13)C) isotope effects not being the signal shifts due to the conformational trans-gauche transition.