Calculation of optical rotation with time-periodic magnetic-field-dependent basis functions in approximate time-dependent density-functional theory

We report the implementation of a method for the calculation of optical rotation. This method is based on the time-dependent density-functional theory and utilizes time-periodic magnetic-field-dependent basis functions. The calculations are based on a density fit. It is demonstrated that additional terms in the analytical expression appearing from derivatives of the approximated Coulomb potential are necessary to provide the gauge-origin independence of the results within a given numerical accuracy. Contributions from these terms also restore the symmetry between the electric and magnetic perturbations in the optical rotation tensor.     

Calculation of static and dynamic linear magnetic response in approximate time-dependent density functional theory

We report implementations and results of time-dependent density functional calculations (i) of the frequency-dependent magnetic dipole-magnetic dipole polarizability, (ii) of the (observable) translationally invariant linear magnetic response, and (iii) of a linear intensity differential (LID) which includes the dynamic dipole magnetizability. The density functional calculations utilized density fitting. For achieving gauge-origin independence we have employed time-periodic magnetic-field-dependent basis functions as well as the dipole velocity gauge, and have included explicit density-fit related derivatives of the Coulomb potential. We present the results of calculations of static and dynamic magnetic dipole-magnetic dipole polarizabilities for a set of small molecules, the LID for the SF6 molecule, and dispersion curves for M-hexahelicene of the origin invariant linear magnetic response as well as of three dynamic polarizabilities: magnetic dipole-magnetic dipole, electric dipole-electric dipole, and electric dipole-magnetic dipole. We have also performed comparison of the linear magnetic response and magnetic dipole-magnetic dipole polarizability over a wide range of frequencies for H2O and SF6.     

An analytical derivative procedure for the calculation of vibrational Raman optical activity spectra

We present an analytical time-dependent Hartree-Fock algorithm for the calculation of the derivatives of the electric dipole-magnetic dipole polarizability with respect to atomic Cartesian coordinates. Combined with analogous procedures to determine the derivatives of the electric dipole-electric dipole and electric dipole-electric quadrupole polarizabilities, it enables a fully analytical evaluation of the three frequency-dependent vibrational Raman optical activity (VROA) invariants within the harmonic approximation. The procedure employs traditional non-London atomic orbitals, and the gauge-origin dependence of the VROA intensities has, therefore, been assessed for the commonly used aug-cc-pVDZ and rDPS:3-21G basis sets.     

Gauge-origin independent calculation of magnetizabilities and rotational g tensors at the coupled-cluster level

An implementation of the gauge-origin independent calculation of magnetizabilities and rotational g tensors at the coupled-cluster (CC) level is presented. The properties of interest are obtained as second derivatives of the energy with respect to the external magnetic field (in the case of the magnetizability) or with respect to magnetic field and rotational angular momentum (in the case of the rotational g tensor), while gauge-origin independence and fast basis-set convergence are ensured by using gauge-including atomic orbitals (London atomic orbitals) as well as their extension to treat rotational perturbations (rotational London atomic orbitals). The implementation within our existing CC analytic second-derivative code is described, focusing on the required modifications concerning integral evaluation and treatment of the unperturbed and perturbed two-particle density matrices. An extensive set of test calculations for LiH and BH (up to the full configuration-interaction limit), for a series of simple hydrides (HF, H(2)O, NH(3), and CH(4)) as well as the more challenging molecules CO, N(2), and O(3) [employing the CC singles and doubles (CCSD) and the CCSD approximation augmented by a perturbative treatment of triple excitations] demonstrates the importance of electron correlation for high-accuracy predictions of magnetizabilities and rotational g tensors.     

Dynamic polarizabilities and raman intensities of CO,N2, HCl and Cl2

We have performed first-order (identical to coupled Hartree-Fock) and second-order polarization propagator calculations of the dynamic dipole polarizability tensor for the CO, N₂, HCl and Cl₂ molecules. The derivatives of the polarizability tensor with respect to the internuclear distance at the equilibrium internuclear separation are compared with related data obtained from non-resonance Raman spectra. In most cases the correlation contributions beyond the coupled Hartree-Fock approximation have a more pronounced effect on the derivatives of the polarizability tensor than on the polarizability tensor itself. At both levels of approximation we found that derivatives of the dynamic polarizability tensor with respect to the internuclear separation increase with 10–15% for variation of the frequency from 0 (static polarizability) to about 28500 cm^?1. The depolarization ratio calculated from the polarizability derivatives shows no variation with frequency in the same frequency range.     

Electro-optical parameters of bond polarizability model for aluminosilicates

Electro-optical parameters (EOPs) of bond polarizability model (BPM) for aluminosilicate structures were derived from quantum-chemical DFT calculations of molecular models. The tensor of molecular polarizability and the derivatives of the tensor with respect to the bond length are well reproduced with the BPM, and the EOPs obtained are in a fair agreement with available experimental data. The parameters derived were found to be transferable to larger molecules. This finding suggests that the procedure used can be applied to systems with partially ionic chemical bonds. The transferability of the parameters to periodic systems was tested in molecular dynamics simulation of the polarized Raman spectra of alpha-quartz. It appeared that the molecular Si-O bond EOPs failed to reproduce the intensity of peaks in the spectra. This limitation is due to large values of the longitudinal components of the bond polarizability and its derivative found in the molecular calculations as compared to those obtained from periodic DFT calculations of crystalline silica polymorphs by Umari et al. (Phys. Rev. B 2001, 63, 094305). It is supposed that the electric field of the solid is responsible for the difference of the parameters. Nevertheless, the EOPs obtained can be used as an initial set of parameters for calculations of polarizability related characteristics of relevant systems in the framework of BPM.     

Polarizable continuum model study of solvent effects on electronic circular dichroism parameters

We present an implementation of the polarizable continuum model for the calculation of solvent effects on electronic circular dichroism spectra. The computational model used is density functional theory in the length-gauge formulation, and gauge-origin independence is ensured through the use of London atomic orbitals. Results of calculations carried out for methyloxirane and bicyclic ketones, camphor, norcamphor, norbornenone, and fenchone are presented, and the theoretically obtained solvent effects are compared with experimental observations.     

Time-dependent density functional calculations of optical rotatory dispersion including resonance wavelengths as a potentially useful tool for determining absolute configurations of chiral molecules

The optical rotations for six organic molecules (verbenone, fenchone, camphor, nopinone, Tr?ger's base, dimethyl-cyclopropane) and the transition metal complex [Co(en)(3)](3+) were calculated as a function of wavelength using time-dependent density functional theory (TDDFT). In the calculations, a realistic behavior of the optical rotation in the vicinity of an electronic transition was obtained by using a phenomenological damping parameter of the order of 0.2 eV (0.007 au). In comparison with experiment, for the molecules studied here the sign and order of magnitude of the optical rotation as well as the excitation energies were reasonably well reproduced in most computations. These findings apply to the investigated wavelength ranges typically between about 200 and 650 nm even when using comparatively small basis sets. Such calculations might therefore routinely be applied to help assigning the absolute configurations of chiral molecules. Supplementary calculations of the circular dichroism (CD) and comparison with experimental CD were used for further assessment of the optical rotation calculations. In particular, a combined study of optical rotation and CD turned out to be useful in cases where the optical rotatory dispersion in a specific energy range exhibits a considerable blue or red shift or where it is difficult to reproduce because of an interplay of several competing Cotton effects. The influence of basis set, density functional, and the damping parameter was also investigated.     

Resonance Raman spectra of uracil based on Kramers-Kronig relations using time-dependent density functional calculations and multireference perturbation theory

The use of time-dependent density functional calculations for the optimization of excited-state structures and the subsequent calculation of resonance Raman intensities within the transform-theory framework is compared to calculations of Hartree-Fock/configuration interaction singles-type (CIS). The transform theory of resonance Raman scattering is based on Kramers-Kronig relations between polarizability tensor components and the optical absorption. Stationary points for the two lowest excited singlet states of uracil are optimized and characterized by means of numerical differentiation of analytical excited-state gradients. It is shown that the effect of electron correlation leads to substantial modifications of the relative intensities. Calculations of vibrational frequencies for ground and excited states are carried out, which show that the neglect of Duschinsky mixing and the assumption of equal wave numbers for ground and excited state are not in all cases good approximations. We also compare the transform-theory resonance Raman intensities with those obtained within a simple approximation from excited-state gradients at the ground-state equilibrium position, and find that they are in qualitative agreement in the case of CIS, but show some important differences in calculations based on density functional theory. Since the results from CIS calculations are in better agreement with experiment, we also present approximate resonance Raman spectra obtained using excited-state gradients from multireference perturbation theory calculations, which confirm the CIS gradients.     

A physical model for the longitudinal polarizabilities of polymer chains

The aim of this work is to provide a physical model to relate the polarizability per unit cell of oligomers to that of their corresponding infinite polymer chains. For this we propose an extrapolation method for the polarizability per unit cell of oligomers by fitting them to a physical model describing the dielectric properties of polymer chains. This physical model is based on the concept of a dielectric needle in which we assume a polymer chain to be well described by a cylindrically shaped nonconducting rod with a radius much smaller than its length. With this model we study in which way the polarizability per unit cell approaches the limit of the infinite chain. We show that within this model the macroscopic contribution of the induced electric field to the macroscopic electric field vanishes in the limit of an infinite polymer chain, i.e., there is no macroscopic screening. The macroscopic electric field becomes equal to the external electric field in this limit. We show that this identification leads to a relation between the polarizability per unit cell and the electric susceptibility of the infinite polymer chain. We test our dielectric needle model on the polarizability per unit cell of oligomers of the hydrogen chain and polyacetylene obtained earlier using time-dependent current-density-functional theory in the adiabatic local-density approximation and with the Vignale-Kohn functional. We also perform calculations using the same theory on truly infinite polymer chains by employing periodic boundary conditions. We show that by extrapolating the oligomer results according to our dielectric needle model we get good agreement with our results from calculations on the corresponding infinite polymer chains.     

Nuclear magnetic resonance J coupling constant polarizabilities of hydrogen peroxide: A basis set and correlation study

In this article, we present the so far most extended investigation of the calculation of the coupling constant polarizability of a molecule. The components of the coupling constant polarizability are derivatives of the nuclear magnetic resonance (NMR) indirect nuclear spin–spin coupling constant with respect to an external electric field and play an important role for both chiral discrimination and solvation effects on NMR coupling constants. In this study, we illustrate the effects of one‐electron basis sets and electron correlation both at the level of density functional theory as well as second‐order polarization propagator approximation for the small molecule hydrogen peroxide, which allowed us to perform calculations with the largest available basis sets optimized for the calculation of NMR coupling constants. We find a systematic but rather slow convergence with the one‐electron basis set and that augmentation functions are required. We observe also large and nonsystematic correlation effects with significant differences between the density functional and wave function theory methods. © 2012 Wiley Periodicals, Inc.     

Time-dependent density functional theory for calculating origin-independent optical rotation and rotatory strength tensors

An approach to calculate origin-independent electronic chiroptical property tensors using time-dependent density functional theory (TDDFT) and gauge-including atomic orbital (GIAO) basis sets is evaluated. Computations of origin-dependent optical rotation tensors and of rotatory strengths needed to simulate circular dichroism spectra are presented. The optical rotation tensor computations employ solutions of coupled perturbed Kohn-Sham equations for a dynamic electric field and a static magnetic field. Because the magnetic field is time independent, the GIAO treatment is somewhat simplified compared to a previously reported method, at some added computational cost if hybrid functionals are employed. GIAO rotatory strengths are also calculated, using transition density matrices from a standard TDDFT excitation energy module. A new implementation in the NWChem quantum chemistry package is employed for representative computations of origin-invariant chiroptical response tensors for methyloxirane, norbornenone, and the ketosteroid androstadienone. For the steroid molecule the vibrational structure of the CD spectrum is modeled explicitly by using calculated Franck-Condon factors. The agreement with experiment is favorable.     

Lagrangian approach to molecular vibrational Raman intensities using time-dependent hybrid density functional theory

The authors propose a new route to vibrational Raman intensities based on analytical derivatives of a fully variational polarizability Lagrangian. The Lagrangian is constructed to recover the negative frequency-dependent polarizability of time-dependent Hartree-Fock or adiabatic (hybrid) density functional theory at its stationary point. By virtue of the variational principle, first-order polarizability derivatives can be computed without using derivative molecular orbital coefficients. As a result, the intensities of all Raman-active modes within the double harmonic approximation are obtained at approximately the same cost as the frequency-dependent polarizability itself. This corresponds to a reduction of the scaling of computational expense by one power of the system size compared to a force constant calculation and to previous implementations. Since the Raman intensity calculation is independent of the harmonic force constant calculation more, computationally demanding density functionals or basis sets may be used to compute the polarizability gradient without much affecting the total time required to compute a Raman spectrum. As illustrated for fullerene C60, the present approach considerably extends the domain of molecular vibrational Raman calculations at the (hybrid) density functional level. The accuracy of absolute and relative Raman intensities of benzene obtained using the PBE0 hybrid functional is assessed by comparison with experiment.     

Static electric dipole polarizability of small sodium aggregates

We present pseudopotential local-spin-density calculations of the static electric polarizability of sodium dimers and trimers and their respective cations. The electronic polarizabilities are obtained from self-consistent calculations in the presence of an external electric field, which is kept sufficiently small to avoid non-linear effects. The calculated polarizability tensor has a strong anisotropy directly related to the geometric and electronic structures of the molecules, the anisotropy being larger for the neutral clusters. The polarizabilities are averaged over the vibrational motion and rotations of the aggregates in order to be compared with the experimental measurements. The obtained values show an improvement in the agreement with experiment with respect to the values calculated in the spherical approximation.     

Density-functional theory calculations of optical rotatory dispersion in the nonresonant and resonant frequency regions

The complex linear response function, which can be employed for calculations of second-order molecular properties in regions of strong absorption, is here extended to encompass the mixed electric-dipole-magnetic-dipole polarizability. The mixed electric-dipole-magnetic-dipole polarizability determines the optical rotation and, when absorption is taken into account, the full anomalous optical rotatory dispersion (ORD) spectra of chiral molecules can be calculated using first-principle quantum-chemical methods. Gauge-origin independence of the results is ensured through the use of London atomic orbitals. To illustrate the importance of taking the absorption process properly into account, we here apply this methodology to the study of the anomalous ORD of hydrogen peroxide, 3R-methylcyclohexanone, 4R-1,1-dimethyl-[3]-(1,2)-ferrocenophan-2-on, and the D(2) isomer of the C(84) fullerene.     

Dynamic polarizability of many-electron systems within a time-dependent density-functional theory

Using the time-dependent extension, proposed by us recently, of the Hohenberg—Kohn—Sham density-functional theory, in the presence of an oscillating electric field, we suggest a Karplus—Kolker-type variation—perturbation method for the calculation of dynamic 2^L-pole polarizability of many-electron systems. As an illustration of the present density-functional formalism, the frequency-dependent dipole polarizability of He atom has been calculated in the frequency range 0 ? ω ? 0.65 au, using local density forms of the exchange and correlation potentials. For ω = 0, the results are numerically better than recent density-functional calculations of the static dipole polarizability of He. The corresponding hydrodynamical formulation, which employs the single-particle density as a basic variable, is also presented.     

Geometrical derivatives of dipole moments and polarizabilities

Molecular dipole moments and polarizabilities, as well as their geometrical derivatives, are given analytical expressions for multiconfiguration self-consistent-field and configuration interaction wavefunctions. By considering the response of the electronic wavefunction induced by electric field and geometrical displacement terms in the Hamiltonian, the response of the total electronic energy to these terms is analyzed. The dipole moment and polarizability are then identified through the factors in the energy which are linear and quadratic in the electric field, respectively. Derivatives with respect to molecular deformation are obtained by identifying factors in these moments which are linear, quadratic, etc., in the distortion parameter. The analytical derivative expressions obtained here are compared to those which arise through finite-difference calculations, and it is shown how previous configuration-interaction-based finite difference dipole moment and polarizability derivatives are wrong. The proper means of treating such derivatives are detailed.     

A Density Functional Theory Study of Optical Rotation in Some Aziridine and Oxirane Derivatives

We present time-dependent density functional theory (TDDFT) calculations of the electronic optical rotation (ORP) for seven oxirane and two aziridine derivatives in the gas phase and in solution and compare the results with the available experimental values. For seven of the studied molecules it is the first time that their optical rotation was studied theoretically and we have therefore investigated the influence of several settings in the TDDFT calculations on the results. This includes the choice of the one-electron basis set, the exchange-correlation functional or the particular polarizable continuum model (PCM). We can confirm that polarized quadruple zeta basis sets augmented with diffuse functions are necessary for converged results and find that the aug-pc-3 basis set is a viable alternative to the frequently employed aug-cc-pVQZ basis set. Based on our study, we cannot recommend the generalized gradient functional KT3 for calculations of the ORP in these compounds, whereas the hybrid functional PBE0 gives results quite similar to the long-range correct CAM−B3LYP functional. Finally, we observe large differences in the solvent effects predicted by the integral equation formalism of PCM and the SMD variant of PCM. For the majority of solute/solvent combinations in this study, we find that the SMD model in combination with the PBE0 functional and the aug-pc-3 basis set gives the best agreement with the experimental values.