Molecular Magnetizabilities: Zero-Point Vibrational Effects and the Breakdown of Pascal's Rule

Calculation of Diamagnetic Susceptibility Tensors of Organic Crystals: From Coronene to Pharmaceutical Polymorphs

Understanding why crystallization in strong magnetic fields can lead to new polymorphs requires methods to calculate the diamagnetic response of organic molecular crystals. We develop the calculation of the macroscopic diamagnetic susceptibility tensor, χ^cryst, for organic molecular crystals using periodic density functional methods. The crystal magnetic susceptibility tensor, χ^cryst, for all experimentally known polymorphs, and its molecular counterpart, χ^mol, are calculated for flexible pharmaceuticals such as carbamazepine, flufenamic acid, and chalcones, and rigid molecules, such as benzene, pyridine, acridine, anthracene, and coronene, whose molecular magnetic properties have been traditionally studied. A tensor addition method is developed to approximate the crystal diamagnetic susceptibility tensor, χ^cryst, from the molecular one, χ^mol, giving good agreement with those calculated directly using the more costly periodic density functional method for χ^cryst. The response of pharmaceutical molecules and crystals to magnetic fields, as embodied by χ^cryst, is largely determined by the packing in the crystal, as well as the molecular conformation. The anisotropy of χ^cryst can vary considerably between polymorphs though the isotropic terms are fairly constant. The implications for developing a computational method for predicting whether crystallization in a magnetic field could produce a novel or different polymorph are discussed.

NMR shielding constants in group 15 trifluorides

By combining large basis and complete basis set (CBS) extrapolations of the coupled-cluster equilibrium geometry results with rovibrational and relativistic corrections, we demonstrate that it is possible to achieve near-quantitative accuracy for the NMR shielding constants in three group 15 trifluorides – NF₃, PF₃ and AsF₃.

On the calculation of second-order magnetic properties using subsystem approaches in a relativistic framework

The implementation of second-order magnetic properties in a frozen density embedding scheme in a four component relativistic framework is outlined and applied to model H₂X–H₂O systems (X = Se, Te, Po).

Magnitude of zero-point vibrational corrections to optical rotation in rigid organic molecules: a time-dependent density functional study

The zero-point vibrational corrections (ZPVCs) to the optical rotation of 22 rigid organic molecules have been calculated using time-dependent density functional theory with the B3LYP hybrid functional. We outline an implementation for calculating ZPVCs that can be applied with a variety of quantum chemistry programs and methods. It is shown that the ZPVCs to optical rotation have a wide range of values and can be quite significant depending on the molecule. On average, it has been determined that vibrational corrections can account for about 20% of the optical rotation for the equilibrium value. It is also concluded that vibrational effects alone cannot be the only factor in improving the calculated values of optical rotation with respect to experimental data measured in solution.

Vibrational Corrections to Magneto-Optical Rotation: A Computational Study

Vibrational corrections to the Verdet constants of nine molecules (H2, N2, CO, H2O, CH4, benzene, toluene, p-xylene, and o-xylene) were calculated with pure density functional theory (DFT), hybrid DFT, and an approximate coupled-cluster theory. Comparisons are made for the accuracy of the vibrational averages among different methods and with respect to experimental data where available. It is found that vibrational corrections to magneto-optical rotation can be as large as 10% of the equilibrium value. Hybrid DFT with the B3LYP hybrid functional offers reasonable accuracy at a relatively inexpensive computational cost for accurate calculations of vibrationally averaged Verdet constants.

Gauge-origin-independent magnetizabilities of solvated molecules using the polarizable continuum model

We present an implementation of the polarizable continuum model in its integral equation formulation for the calculation of the magnetizabilities of solvated molecules. The gauge-origin independence of the calculated magnetizabilities and the fast basis set convergence are ensured through the use of London atomic orbitals. Our implementation can use Hartree-Fock and multiconfigurational self-consistent-field (MCSCF) wave functions as well as density-functional theory including hybrid functionals such as B3LYP. We present the results of dielectric continuum effects on water and pyridine using MCSCF wave functions, as well as dielectric medium effects on the magnetizability of the aromatic amino acids as a model for how a surrounding protein environment affects the magnetizability of these molecules. It is demonstrated that the dielectric medium effects on the magnetizability anisotropies of the aromatic amino acids may be substantial, being as large as 25% in the case of tyrosine.

A computational study of some electric and magnetic properties of gaseous BF3 and BCl3

We present the results of an extended computational study of the electric and magnetic properties connected to Cotton-Mouton birefringences, on the trifluoro- and trichloroborides in the gas phase. The electric dipole polarizabilities, magnetizabilities, quadrupole moments, and higher-order hypersusceptibilities--expressed as quadratic and cubic frequency-dependent response functions--are computed within Hartree-Fock, density-functional, and coupled-cluster response theories employing singly and doubly augmented correlation-consistent basis sets and London orbitals in the magnetic property calculations. The results, which illustrate the capability of time-dependent density-functional theory for electron-rich systems, are compared with available experimental data. Revised values of both experimentally derived quadrupole moment of BF3, 2.72 +/- 0.15 a.u., and magnetizability anisotropy of BCl3, -0.45 +/- 0.09 a.u., both obtained in birefringence experiments that neglect the effects of higher-order hypersusceptibilities, are presented. In the theoretical limit the traceless quadrupole moments of BF3 and BCl3 are determined to be 3.00 +/- 0.01 and 0.71 +/- 0.01 a.u., respectively.

Density-functional theory study of electric and magnetic properties of hexafluorobenzene in the vapor phase

A series of electric and magnetic properties of hexafluorobenzene have been calculated, including the electric dipole polarizability, magnetizability, electric quadrupole moment, and nonlinear mixed electric dipole-magnetic dipole-electric quadrupole hyperpolarizabilities needed to obtain estimates of the Kerr, Cotton-Mouton, Buckingham, Jones, and magnetoelectric birefringences in the vapor phase. Time-dependent density-functional theory was employed for the calculation of linear-, quadratic, and cubic response functions. A number of density functionals have been considered, along with Sadlej's triple-zeta basis set and the augmented correlation-consistent polarized valence double zeta and augmented correlation-consistent polarized valence triple zeta basis sets. Comparisons have been made with experiment where possible. The analysis of results allows for an assessment of the capability of time-dependent density-functional theory for high-order electromagnetic properties of an electron-rich system such as hexafluorobenzene.