We present an analysis of the performances of a parameter free density functional model (PBE0) obtained combining the so called PBE generalized gradient functional with a predefined amount of exact exchange. The results obtained for structural, thermodynamic, kinetic and spectroscopic (magnetic, infrared and electronic) properties are satisfactory and not far from those delivered by the most reliable functionals including heavy parameterization. The way in which the functional is derived and the lack of empirical parameters fitted to specific properties make the PBE0 model a widely applicable method for both quantum chemistry and condensed matter physics.

Toward reliable density functional methods without adjustable parameters: The PBE0 model

6.2.2. Definition of Effective Properties 3064 6.3. Response Properties to Magnetic Fields 3066 6.3.1. Nuclear Shielding 3066 6.3.2. Indirect Spin−Spin Coupling 3067 6.3.3. EPR Parameters 3068 6.4. Properties of Chiral Systems 3069 6.4.1. Electronic Circular Dichroism (ECD) 3069 6.4.2. Optical Rotation (OR) 3069 6.4.3. VCD and VROA 3070 7. Continuum and Discrete Models 3071 7.1. Continuum Methods within MD and MC Simulations 3072

/pdf/quantum-mechanical-continuum-solvation-models-3edv963g14.pdf

Quantum mechanical continuum solvation models.

The halogen bond occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a halogen atom in a molecular entity and a nucleophilic region in another, or the same, molecular entity. In this fairly extensive review, after a brief history of the interaction, we will provide the reader with a snapshot of where the research on the halogen bond is now, and, perhaps, where it is going. The specific advantages brought up by a design based on the use of the halogen bond will be demonstrated in quite different fields spanning from material sciences to biomolecular recognition and drug design.

The Halogen Bond

Continuum solvation models are appealing because of the simplified yet accurate description they provide of the solvent effect on a solute, described either by quantum mechanical or classical methods. The polarizable continuum model (PCM) family of solvation models is among the most widely used, although their application has been hampered by discontinuities and singularities arising from the discretization of the integral equations at the solute-solvent interface. In this contribution we introduce a continuous surface charge (CSC) approach that leads to a smooth and robust formalism for the PCM models. We start from the scheme proposed over ten years ago by York and Karplus and we generalize it in various ways, including the extension to analytic second derivatives with respect to atomic positions. We propose an optimal discrete representation of the integral operators required for the determination of the apparent surface charge. We achieve a clear separation between “model” and “cavity” which, together with simple generalizations of modern integral codes, is all that is required for an extensible and efficient implementation of the PCM models. Following this approach we are now able to introduce solvent effects on energies, structures, and vibrational frequencies (analytical first and second derivatives with respect to atomic coordinates), magnetic properties (derivatives with respect of magnetic field using GIAOs), and in the calculation more complex properties like frequency-dependent Raman activities, vibrational circular dichroism, and Raman optical activity.

Continuous surface charge polarizable continuum models of solvation. I. General formalism

Halogen bonding refers to the non-covalent interactions of halogen atoms X in some molecules, RX, with negative sites on others. It can be explained by the presence of a region of positive electrostatic potential, the sigma-hole, on the outermost portion of the halogen's surface, centered on the R-X axis. We have carried out a natural bond order B3LYP analysis of the molecules CF(3)X, with X = F, Cl, Br and I. It shows that the Cl, Br and I atoms in these molecules closely approximate the [Formula: see text] configuration, where the z-axis is along the R-X bond. The three unshared pairs of electrons produce a belt of negative electrostatic potential around the central part of X, leaving the outermost region positive, the sigma-hole. This is not found in the case of fluorine, for which the combination of its high electronegativity plus significant sp-hybridization causes an influx of electronic charge that neutralizes the sigma-hole. These factors become progressively less important in proceeding to Cl, Br and I, and their effects are also counteracted by the presence of electron-withdrawing substituents in the remainder of the molecule. Thus a sigma-hole is observed for the Cl in CF(3)Cl, but not in CH(3)Cl.

Halogen bonding: the σ-hole

The direct (recomputation of two‐electron integrals) implementation of the gauge‐including atomic orbital (GIAO) and the CSGT (continuous set of gauge transformations) methods for calculating nuclear magnetic shielding tensors at both the Hartree‐Fock and density functional levels of theory are presented. Isotropic 13C, 15N, and 17O magnetic shielding constants for several molecules, including taxol (C47H51NO14 using 1032 basis functions) are reported. Shielding tensor components determined using the GIAO and CSGT methods are found to converge to the same value at sufficiently large basis sets; however, GIAO shielding tensor components for atoms other than carbon are found to converge faster with respect to basis set size than those determined using the CSGT method for both Hartree‐Fock and DFT. For molecules where electron correlation effects are significant, shielding constants determined using (gradient‐corrected) pure DFT or hybrid methods (including a mixture of Hartree‐Fock exchange and DFT exchange...

/pdf/a-comparison-of-models-for-calculating-nuclear-magnetic-20e4t1dzjd.pdf

A comparison of models for calculating nuclear magnetic resonance shielding tensors

Properties of atoms in molecules: atomic volumes

A new theory of solvent effects on the optical rotations of chiral molecules is presented. The frequency-dependent electric dipole−magnetic dipole polarizability, βαβ(ν), is calculated using density functional theory (DFT). Solvent effects are included using the polarizable continuum model (PCM). DFT/PCM calculations of sodium D line specific rotations, [α]D, have been carried out for seven conformationally rigid chiral organic molecules (fenchone, camphor, α-pinene, β-pinene, camphorquinone, verbenone, and methyloxirane) for a diverse set of seven solvents (cyclohexane, carbon tetrachloride, benzene, chloroform, acetone, methanol, and acetonitrile). The predicted variation in [α]D for the solvents cyclohexane, acetone, methanol, and acetonitrile are in excellent agreement with experiment for all seven molecules. For the solvents carbon tetrachloride, benzene, and chloroform, agreement is much poorer. Since only electrostatic solute−solvent interactions are included in the PCM, our results lead to the con...

Polarizable Continuum Model (PCM) Calculations of Solvent Effects on Optical Rotations of Chiral Molecules

Ab initio calculation of atomic axial tensors and vibrational rotational strengths using density functional theory

Ab initio Hartree−Fock (HF) and density functional theory (DFT) calculations of the electric dipole−magnetic dipole polarizability βαβ are reported for the chiral molecules methyloxirane (1) and tr...

James R. Cheeseman

Papers

A comparison of models for calculating nuclear magnetic resonance shielding tensors

Properties of atoms in molecules: atomic volumes

Polarizable Continuum Model (PCM) Calculations of Solvent Effects on Optical Rotations of Chiral Molecules

Ab initio calculation of atomic axial tensors and vibrational rotational strengths using density functional theory

Hartree−Fock and Density Functional Theory ab Initio Calculation of Optical Rotation Using GIAOs: Basis Set Dependence