Theoretical analysis of the hydrogen bond interaction between acetone and water

Probing microhydration-induced effects on carbonyl compounds

Characterizing the microhydration of organic molecules is a crucial step in understanding many phenomena relevant to atmospheric, biological, and industrial applications. However, its precise experimental and theoretical description remains a challenge. For four organic solutes containing a CO bond, and included in the recent HyDRA challenge [T. L. Fischer, M. Bödecker, A. Zehnacker-Rentien, R. A. Mata and M. A. Suhm, Phys. Chem. Chem. Phys., 2022, 24, 11442-11454.], we performed a detailed study of different monohydrate isomers and their properties; these were cyclooctanone (CON), 1,3-dimethyl-2-imidazolidinon (DMI), methyl lactate (MLA), and 2,2,2-trifluoroacetophenone (TPH) molecules. As reported in the literature, the O-H elongation shift of the water molecule appears to be a good candidate for characterizing complexation-induced effects. We also show that CO elongation shift and UV-vis spectroscopy can be successfully used for these purposes. Besides, we present a comparative analysis of the strengths of non-covalent interactions within these monohydrated complexes based on interpretative tools of quantum chemistry, including topological analysis of electron density (ρ), topological analysis of electron pairing function, and analysis of the core-valence bifurcation index (CVBI), which exhibits a close linear dependency on ρ. Accordingly, a classification of intermolecular water-solute interactions is proposed.

Quantum Chemical Approaches to the Calculation of NMR Parameters: From Fundamentals to Recent Advances

Quantum chemical methods for the calculation of indirect NMR spin–spin coupling constants and chemical shifts are always in progress. They never stay the same due to permanently developing computational facilities, which open new perspectives and create new challenges every now and then. This review starts from the fundamentals of the nonrelativistic and relativistic theory of nuclear magnetic resonance parameters, and gradually moves towards the discussion of the most popular common and newly developed methodologies for quantum chemical modeling of NMR spectra.

The cooperative effect of the sorption of volatile molecules into metal-oxide frameworks as a function of the dielectric constant

The separation of volatile molecules into metal-oxide frameworks (MOxF) is possible with the combination of two approaches: spatial confinement and sorption-desorption equilibria. We have shown that nanoporous MOxF, namely giant POMs with Keplerate or toroidal structures, have a specific sorption behaviour towards different volatile substances that can be generalized through the dependence of the limit sorption on the dielectric constant and molecular polarizability. Assisted by the Hildebrand and Hansen solubility parameters and the Snyder polarity index, an analysis of sorption behaviour shows that the contributions of hydrogen bonding, hydrophobic forces and polar and dispersion interactions can be expressed through a cooperative phenomenon (the formation of a quasi-liquid layer) correlated with dielectric constant as a major parameter. The obtained data can be used to predict the separation of volatile substances using MOxF-embedded hybrid materials, such as selective membranes.

Aldehydes as powerful initiators for photochemical transformations

Photochemistry, the use of light to promote organic transformations, has been known for more than a century but only recently has revolutionized the way modern chemists are thinking. Except from transition metal-based complexes, small organic molecules have been introduced as catalysts or initiators. In this review, we summarize the potential that (aromatic or aliphatic) aldehydes have as photoinitiators. The photophysical properties and photoreactivity of benzaldehyde are initially provided, followed by applications of aldehydes as initiators for polymerization reactions. Finally, the applications to date regarding aldehydes as photoinitiators in organic synthesis are presented.

Prolongation of secondary drying step of phospholipid lyophilization greatly improves acidolysis reactions catalyzed by immobilized lecitase ultra

Conjugated linolenic acid (CLA) has important health effects, and their phospholipids forms (PL) are advantageous vehicles of this bioactive agent. Acidolysis of soybean phosphatidylcholine (PC) with conjugated linolenic acid (CLA) catalyzed by Lecitase® Ultra immobilized on Duolite A658 was studied. This reaction is typically limited by hydrolysis, producing 60-90 % of lyso- and glycero-PC and yielding low the process efficiency. Drying the amphiphilic PC material was found the key factor for maximal diacylglycerol phosphatidylcholines (PC) production and different drying approaches were studied to maximize the formation of PC rich in CLA in a solvent free process. PC lyophilization for 24 h getting dry solid appearance (PC; 783 ± 11 mg water/Kg PC) or other standard protocols to reduce water content/activity of reaction medium, did not improve the reaction performance. However, adding 4 extra days to the second step of lyophilization at high vacuum (1 Pa) and moderate temperature (20 °C), followed by further PC dehydration with molecular sieves, drastically reduced the hydrolysis process by achieving a extensive PC dehydration (279 ± 4 mg water/Kg PC), obtaining for the first time >99% molar yield of diacyl-PC product. After 24 h of reaction, a diacyl-PC product with 72.3% CLA content was obtained. PC molecules containing two CLA were the major species formed.

Multiple active zones in hybrid QM/MM molecular dynamics simulations for large biomolecular systems

A new QM/MM molecular dynamics approach that can deal with the dynamics of large real systems involving several simultaneous active zones is presented.

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.

A quantum mechanical/molecular dynamics/mean field study of acrolein in aqueous solution: analysis of H bonding and bulk effects on spectroscopic properties

A novel molecular dynamics methodology recently proposed by our group [Rega et al., Chem. Phys. Lett. 422, 367 (2006)], which is based on an integrated hybrid potential rooted in high level quantum mechanical methods using localized basis functions and nonperiodic boundary conditions, has been applied to study acrolein in aqueous solution. The solute structural rearrangement and its hydrogen-bonding pattern due to the interactions with water have been analyzed in some detail. Moreover, the solvent effects on the UV n-->pi* vertical transition and on the NMR 13C and 17O shielding constants of acrolein have been investigated theoretically by performing a posteriori quantum mechanical calculations on a statistically significant number of snapshots extracted from both gas-phase and aqueous solution simulations. Results show that such effective computational strategy can be successfully used to improve our understanding, at atomic level, of important spectroscopic observables.

Coupled cluster calculation of the n --> pi* electronic transition of acetone in aqueous solution

The combined linear response coupled cluster/molecular mechanics (CC/MM) scheme including mutual polarization effects in the coupling Hamiltonian is applied together with supermolecular CC methods to the study of the gas-to-aqueous solution blue shift of the n --> pi* excitation energy in acetone. The aug-cc-pVDZ basis set is found to be adequate for the calculation of this excitation energy. In the condensed phase, the shift in the excitation energy is obtained by statistical averaging over 800 solute-solvent configurations extracted from a molecular dynamics simulation. We find the shift to be around 1100-1200 cm(-1) depending on the specific model used to describe solvent polarization. The importance of including explicit polarization in both the molecular dynamics simulation as well as the CC/MM calculations is emphasized. Furthermore, the significant dependence of the excitation energy on the CO bond length of acetone is discussed.

Probing supercritical water with the n-π* transition of acetone: A Monte Carlo/quantum mechanics study

The n-pi(*) electronic transition of acetone is a convenient and important probe to study supercritical water. The solvatochromic shift of this transition in supercritical water (adopting the experimental condition of P=340.2 atm and T=673 K) has been studied theoretically using Metropolis NPT Monte Carlo (MC) simulation and quantum mechanics (QM) calculations based on INDO/CIS and TDDFT-B3LYP6-31+G(d) methods. MC simulations are used to analyze hydration shells, solute-solvent interaction, and for generating statistically relevant configurations for subsequent QM calculations of the n-pi(*) transition of acetone. The results show that the average number of hydrogen bonds between acetone and water is essentially 13 of that in normal water condition of temperature and pressure. But these hydrogen bonds have an important contribution in the solute stabilization and in the solute-solvent interaction. In addition, they respond for nearly half of the solvatochromic shift. The INDO/CIS calculations explicitly considering all valence electrons of the water molecules, using different solvation shells, up to the third shell (170 water molecules), give a solvatochromic shift of 670+/-36 cm(-1) in very good agreement with the experimentally inferred result of 500-700 cm(-1). It is found that the solvatochromic effect on n-pi(*) transition of acetone in the supercritical condition is essentially given by the first solvation shell. The time-dependent density-functional theory (TDDFT) calculations are also performed including all solvent molecules up to the third shell, now represented by point charges. This TDDFT-B3LYP6-31+G(d) also gives a good but slightly overestimated result of 825+/-65 cm(-1). For comparison the same study is also made for acetone in water at normal condition. Finally, all average results reported here are statistically converged.

A sequential Monte Carlo quantum mechanics study of the hydrogen-bond interaction and the solvatochromic shift of the n–π* transition of acrolein in water

The sequential Monte Carlo (MC) quantum mechanics (QM) methodology, using time-dependent density-functional theory (TD-DFT), is used to study the solvatochromic shift of the n-pi* transition of trans-acrolein in water. Using structures obtained from the isothermal-isobaric Metropolis MC simulation TD-DFT calculations, within the B3LYP functional, are performed for the absorption spectrum of acrolein in water. In the average acrolein makes one hydrogen bond with water and the hydrogen-bond shell is responsible for 30% of the total solvatochromic shift, considerably less than the shift obtained for the minimum-energy configurations. MC configurations are sampled after analysis of the statistical correlation and 100 configurations are extracted for subsequent QM calculations. All-electron TD-DFT B3LYP calculations of the absorption transition including acrolein and all explicit solvent molecules within the first hydration shell, 26 water molecules, give a solvatochromic shift of 0.18 +/- 0.11 eV. Using simple point charges to represent the solvent the shifts are calculated for the first, second, and third solvation shells. The results converge for the calculated shift of 0.20 +/- 0.10 eV in very good agreement with the experimentally inferred result of 0.20 +/- 0.05 eV. All average results presented are statistically converged.

A Sequential Molecular Mechanics/Quantum Mechanics Study of the Electronic Spectra of Amides

We report gas-phase electronic spectra of formamide, N-methyformamide, acetamide, and N-methylacetamide at 300 K calculated using a combination of classical molecular dynamics and time-dependent density functional theory (TDDFT). In comparison to excitation energies computed using the global minima structures, the valence npi* and pi(nb)pi* states show a significant red-shift of 0.1-0.35 eV, while smaller shifts are found for the n3s and pi(nb)3s Rydberg states. In this work, we have identified the physical origin of these shifts arising from variations of the molecular structure. We present simple relationships between key geometrical parameters and spectral shifts. Consequently, electronic spectra can be generated directly from ground-state structures, without additional quantum chemical calculations. The electronic spectrum of formamide in aqueous solution is computed using TDDFT using an explicit solvent model. This provides a quantitative determination of the condensed-phase spectrum. In general, this study shows that temperature effects can change the predicted excitation energies significantly and demonstrates how electronic spectra at elevated temperatures can be computed in a computationally efficient way.

A computational study of microsolvation effect on ethylene glycol by density functional method

This study focuses on the conformational analysis of ethylene glycol-(water)n (n=1-3) complex by using density functional theory method and the basis set 6-311++G*. Different conformers are reported and the basis set superposition error corrected total energy is -306.767 5171, -383.221 3135, and -459.694 1528 for lowest energy conformer with 1, 2, and 3 water molecules, respectively, with corresponding binding energy -7.75, -15.43, and -36.28 kcal/mol. On applying many-body analysis it has been found that relaxation energy, two-body, three-body energy have significant contribution to the binding energy for ethylene glycol-(water)3 complex whereas four-body energies are negligible. The most stable conformers of ethylene glycol-(water)n complex are the cyclic structures in which water molecules bridge between the two hydroxyl group of ethylene glycol.