scholarly journals Analytical calculation of pressure for confined atomic and molecular systems using the eXtreme‐Pressure Polarizable Continuum Model

2018 ◽  
Vol 39 (26) ◽  
pp. 2243-2250 ◽  
Author(s):  
Roberto Cammi ◽  
Bo Chen ◽  
Martin Rahm
2021 ◽  
Author(s):  
Ariel Gale ◽  
Eugen Hruska ◽  
Fang Liu

Pressure plays essential roles in chemistry by altering structures and controlling chemical reactions. The extreme-pressure polarizable continuum model (XP-PCM) is an emerging method with an efficient quantum mechanical description of small and medium-size molecules at high pressure (on the order of GPa). However, its application to large molecular systems was previously hampered by CPU computation bottleneck: the Pauli repulsion potential unique to XP-PCM requires the evaluation of a large number of electric field integrals, resulting in significant computational overhead compared to the gas-phase or standard-pressure polarizable continuum model calculations. Here, we exploit advances in Graphical Processing Units (GPUs) to accelerate the XP-PCM integral evaluations. This enables high-pressure quantum chemistry simulation of proteins that used to be computationally intractable. We benchmarked the performance using 18 small proteins in aqueous solutions. Using a single GPU, our method evaluates the XP-PCM free energy of a protein with over 500 atoms and 4000 basis functions within half an hour. The time taken by the XP-PCM-integral evaluation is typically 1\% of the time taken for a gas-phase density functional theory (DFT) on the same system. The overall XP-PCM calculations require less computational effort than that for their gas-phase counterpart due to the improved convergence of self-consistent field iterations. Therefore, the description of the high-pressure effects with our GPU accelerated XP-PCM is feasible for any molecule tractable for gas-phase DFT calculation. We have also validated the accuracy of our method on small molecules whose properties under high pressure are known from experiments or previous theoretical studies.


2021 ◽  
Author(s):  
Bo Chen ◽  
K. N. Houk ◽  
Roberto Cammi

Quantum chemical calculations are reported for the thermal dimerizations of 1,3-cyclohexadiene at 1 atm and high pressures up to 6 GPa. Previous experiments [Klärner et al. Angew. Chem. Int. Ed. 1986, 25, 108], based on measured activation energies and activation volumes, suggested concerted mechanisms for the formation of the endo [4+2] cycloadduct and a [6+4]-ene adduct, and stepwise mechanisms for the formation of the exo [4+2] cycloadduct and two [2+2] cycloadducts. Computed activation enthalpies (ωB97XD, CCSD(T) and SC-NEVPT2) of plausible dimerization pathways at 1 atm agree well with the experiment activation energies and the values from previous calculations [Ess et al. J. Org. Chem. 2008, 73, 7586]. High-pressure reaction profiles, computed by the recently-developed extreme pressure-polarizable continuum model (XP-PCM), show that the reduction of reaction barrier is more profound in concerted reactions than in stepwise reactions, which is rationalized on the basis of the volume profiles of different mechanisms. A clear shift of the transition state towards the reactant by high pressure is revealed for the [6+4]-ene reaction by the calculations. The computed activation volumes by XP-PCM agree excellently with the experimental values, confirming the existence of competing mechanisms in the thermal dimerizations of 1,3-cyclohexadiene.


Symmetry ◽  
2021 ◽  
Vol 13 (7) ◽  
pp. 1298
Author(s):  
Ilya G. Shenderovich ◽  
Gleb S. Denisov

The isotopically enriched cyanide anion, (13C≡15N)−, has a great potential as the NMR probe of non-covalent interactions. However, hydrogen cyanide is highly toxic and can decompose explosively. It is therefore desirable to be able to theoretically estimate any valuable results of certain experiments in advance in order to carry out experimental studies only for the most suitable molecular systems. We report the effect of hydrogen bonding on NMR properties of 15N≡13CH···X and 13C≡15NH···X hydrogen bonding complexes in solution, where X = 19F, 15N, and O=31P, calculated at the ωB97XD/def2tzvp and the polarizable continuum model (PCM) approximations. In many cases, the isotropic 13C and 15N chemical shieldings of the cyanide anion are not the most informative NMR properties of such complexes. Instead, the anisotropy of these chemical shieldings and the values of scalar coupling constants, including those across hydrogen bonds, can be used to characterize the geometry of such complexes in solids and solutions. 1J(15N13C) strongly correlates with the length of the N≡C bond.


2021 ◽  
Vol 20 (1) ◽  
pp. 59-68
Author(s):  
Zohreh Khanjari ◽  
Bita Mohtat ◽  
Reza Ghiasi ◽  
Hoorieh Djahaniani ◽  
Farahnaz Kargar Behbahani

This research examined the effects of solvent polarity and temperature on the tautomerization of a carbonitrile molecule at CAM-B3LYP/6-311G (d,p) level of theory. The selected solvents were n-hexane, diethyl ether, pyridine, ethanol, methanol, and water. The solvent effects were examined by the self-consistent reaction field theory (SCRF) based on conductor-like polarizable continuum model (CPCM). The solvent effects were explored on the energy barrier, frontier orbitals energies, and HOMO-LUMO gap. Dependencies of thermodynamic parameters (ΔG and ΔH) on the dielectric constants of solvents were also tested. Specifically, the temperature dependencies of the thermodynamics parameters were studied within 100–1000 K range. The rate constant of the tautomerism reaction was computed from 300 to 1200 K, in the gas phase.


Sign in / Sign up

Export Citation Format

Share Document