PubMed Health⌕ Search

Biomedical subjects

Kenneth Ruud

Publications and source records attributed to Kenneth Ruud.

At least 19 recordsLinked to original sources
PubMed2006

Combined density functional/polarizable continuum model study of magnetochiral birefringence: can theory and experiment be brought to agreement?

The magnetic-field-induced axial birefringence (magnetochirality) of five closed-shell chiral molecules (three substituted oxiranes, carvone, and limonene) is studied at the density functional theory level using Becke's 3-parameter Lee-Yang-Parr functional and frequency-dependent quadratic response theory. The influence of the environment and the conformational distribution on the property is also studied. The environment effects are described by the polarizable continuum model in its integral-equation formulation. The effect of the conformational distribution is investigated by performing calculations on several conformers-for carvone and limonene-followed by Boltzmann averaging. The calculated values for the magnetochiral birefringence are compared to previous ab initio results and experimental data where available. The refined model presented here brings the ab initio values closer to experiment. Still, disagreements remain in some cases and it appears difficult to resolve these discrepancies.

Journal Article↗
PubMed2006

Parallelization of the integral equation formulation of the polarizable continuum model for higher-order response functions.

We present a parallel implementation of the integral equation formalism of the polarizable continuum model for Hartree-Fock and density functional theory calculations of energies and linear, quadratic, and cubic response functions. The contributions to the free energy of the solute due to the polarizable continuum have been implemented using a master-slave approach with load balancing to ensure good scalability also on parallel machines with a slow interconnect. We demonstrate the good scaling behavior of the code through calculations of Hartree-Fock energies and linear, quadratic, and cubic response function for a modest-sized sample molecule. We also explore the behavior of the parallelization of the integral equation formulation of the polarizable continuum model code when used in conjunction with a recent scheme for the storage of two-electron integrals in the memory of the different slaves in order to achieve superlinear scaling in the parallel calculations.

Journal Article↗
PubMed2006

Can Raman optical activity separate axial from local chirality? A theoretical study of helical deca-alanine.

Motivated by experimental work on the distinction of protein secondary structure motifs by Raman optical activity (ROA) spectroscopy, we demonstrate using density functional theory that axial chirality in structures with different local chirality can be filtered out by ROA spectroscopy. To this purpose, two diastereomers of right-handed helical deca-alanine, the (all-S) and the (R,S,R,S,R,S,R,S,R,S) form, are compared. Furthermore, we suggest to interpret calculated ROA spectra of large molecules in terms of vibrational bands rather than individual peaks. This is due to the non-homogeneous effect of the harmonic approximation as well as of the chosen electronic structure method onto the vibrational frequencies, which in a dense region of many vibrations will strongly determine the shape of the spectrum. In addition, the calculated ROA spectrum of (all-S)-deca-alanine is compared to the experimental spectrum of poly-(L)-alanine in solution.

Alanine↗
PubMed2006

Solvent effects on Raman optical activity spectra calculated using the polarizable continuum model.

The integral equation formulation of the polarizable continuum model (IEFPCM) has been extended to the calculation of solvent effects on vibrational Raman optical activity spectra. Gauge-origin independence of the differential scattering intensities of right and left circularly polarized light is ensured through the use of London atomic orbitals. Density functional theory (DFT) calculations have been carried out for bromochlorofluoromethane, methyloxirane, and epichlorhydrin. The results indicate that solvent effects on the ROA differential scattering intensities can be substantial, and vary in sign and magnitude for different vibrational modes. It is demonstrated that both direct and indirect effects are important in determining the total solvent effects on the ROA differential scattering intensities. Local field effects are shown to be in general small, whereas electronic nonequilibrium solvation has a profound effect on the calculated solvent effects compared to an equilibrium solvation model. For molecules with several conformations, the changes in the relative stability of the different conformers also lead to noticeable changes in the ROA spectra.

Journal Article↗
PubMed2006

Superlinear scaling in master-slave quantum chemical calculations using in-core storage of two-electron integrals.

We describe the implementation of a parallel, in-core, integral-direct Hartree-Fock and density functional theory code for the efficient calculation of Hartree-Fock wave functions and density functional theory. The algorithm is based on a parallel master-slave algorithm, and the two-electron integrals calculated by a slave are stored in available local memory. To ensure the greatest computational savings, the master node keeps track of all integral batches stored on the different slaves. The code can reuse undifferentiated two-electron integrals both in the wave function optimization and in the evaluation of second-, third-, and fourth-order molecular properties. Superlinear scaling is achieved in a series of test examples, with speedups of up to 55 achieved for calculations run on medium-sized molecules on 16 processors with respect to the time used on a single processor.

Journal Article↗
PubMed2005

Absolute configuration of C76 from optical rotatory dispersion.

The absolute configuration of C76 has been determined as (+)589-(fC)-C76 , for the first time, by comparing the experimental and predicted optical rotatory dispersion (ORD) patterns. The experimental ORD pattern was derived from the experimental electronic circular dichroism (ECD) spectrum using the Kramers-Kronig (KK) transform. The theoretical ORD spectra were calculated in the resonant region using linear response theory, and also using the KK transform of the theoretical ECD spectrum, at different theoretical levels, namely BHLYP/6-31G*, B3LYP/6-31G*, BLYP/6-31G*, and HF/6-31G*. Good agreement noted between experimental and predicted spectra allows for an unambiguous determination of the absolute configuration.

Journal Article↗
PubMed2005

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.

Journal Article↗
PubMed2005

Second-harmonic generation of solvated molecules using multiconfigurational self-consistent-field quadratic response theory and the polarizable continuum model.

We present the first implementation of the quadratic response function for multiconfigurational self-consistent-field wave functions of solvated molecules described by a polarizable continuum model employing a molecule-shaped cavity. We apply the methodology to the first hyperpolarizability beta and, in particular, the second-harmonic generation process for a series of conjugated push-pull oligomers, as well as for para-nitroaniline. The effect of solvation on the dispersion of the hyperpolarizability and the change in the hyperpolarizability for increasing chain length of the oligomers in vacuum and in solution is considered. The effect of a correlated description is analyzed by comparing the Hartree-Fock hyperpolarizabilities to the multiconfigurational self-consistent-field hyperpolarizabilities. The effect of geometry relaxation in the solvent on the properties of the solvated molecules are also investigated.

Journal Article↗
PubMed2005

Basis set and density functional dependence of vibrational Raman optical activity calculations.

We present an extensive investigation of the dependence of the scattering intensity difference of right and left circularly polarized light observed in vibrational Raman optical activity (VROA) on the choice of basis set and exchange-correlation functional. These dependencies are investigated for five molecules for which accurate experimental data are available: (S)-methyloxirane, (R)-epichlorhydrin, (S)-glycidol, (M)-spiro[2,2]pentane-1,4-diene, and (M)-sigma-[4]-helicene. Calculations are presented using the SVWN exchange-correlation functional (LDA), the BLYP exchange-correlation functional, and the B3LYP hybrid functional, using six different basis sets: the cc-pVDZ, cc-pVTZ, aug-cc-pVDZ, aug-cc-pVTZ, Sadlej's polarized basis set, and a minimal VROA basis set recently proposed by Zuber and Hug. It is demonstrated that results from pure gradient-corrected and hybrid functionals are comparable and that the aug-cc-pVDZ and aug-cc-pVTZ basis sets yield similar results. Furthermore, the combination of the small basis set by Zuber and Hug with an accurate force field represents the best compromise between computational accuracy and computational efficiency.

Journal Article↗
PubMed2005

Density functional theory study of indirect nuclear spin-spin coupling constants with spin-orbit corrections.

This work outlines the calculation of indirect nuclear spin-spin coupling constants with spin-orbit corrections using density functional response theory. The nonrelativistic indirect nuclear spin-spin couplings are evaluated using the linear response method, whereas the relativistic spin-orbit corrections are computed using quadratic response theory. The formalism is applied to the homologous systems H2X (X=O,S,Se,Te) and XH4 (X=C,Si,Ge,Sn,Pb) to calculate the indirect nuclear spin-spin coupling constants between the protons. The results confirm that spin-orbit corrections are important for compounds of the H2X series, for which the electronic structure allows for an efficient coupling between the nuclei mediated by the spin-orbit interaction, whereas in the case of the XH4 series the opposite situation is encountered and the spin-orbit corrections are negligible for all compounds of this series. In addition we analyze the performance of the density functional theory in the calculations of nonrelativistic indirect nuclear spin-spin coupling constants.

Journal Article↗
PubMed2005

Leading-order relativistic effects on nuclear magnetic resonance shielding tensors.

We present perturbational ab initio calculations of the nuclear-spin-dependent relativistic corrections to the nuclear magnetic resonance shielding tensors that constitute, together with the other relativistic terms reported by us earlier, the full leading-order perturbational set of results for the one-electron relativistic contributions to this observable, based on the (Breit-)Pauli Hamiltonian. These contributions are considered for the H(2)X (X = O,S,Se,Te,Po) and HX (X = F,Cl,Br,I,At) molecules, as well as the noble gas (Ne, Ar, Kr, Xe, Rn) atoms. The corrections are evaluated using the relativistic and magnetic operators as perturbations on an equal footing, calculated using analytical linear and quadratic response theory applied on top of a nonrelativistic reference state provided by self-consistent field calculations. The (1)H and heavy-atom nuclear magnetic shielding tensors are compared with four component, nearly basis-set-limit Dirac-Hartree-Fock calculations that include positronic excitations, as well as available literature data. Besides the easy interpretability of the different contributions in terms of familiar nonrelativistic concepts, the accuracy of the present perturbational scheme is striking for the isotropic part of the shielding tensor, for systems including elements up to Xe.

Journal Article↗
PubMed2005

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.

Journal Article↗
PubMed2004

Electronic g-tensors of solvated molecules using the polarizable continuum model.

We present the implementation of density functional response theory combined with the polarizable continuum model (PCM), enabling first principles calculations of molecular g-tensors of solvated molecules. The calculated g-tensor shifts are compared with experimental g-tensor shifts obtained from electron paramagnetic resonance spectra for a few solvated species. The results indicate qualitative agreement between the calculations and the experimental data for aprotic solvents, whereas PCM fails to reproduce the electronic g-tensor behavior for protic solvents. This failure of PCM for protic solvents can be resolved by including into the model those solvent molecules which are involved in hydrogen bonding with the solute. The results for the protic solvents show that the explicit inclusion of the solvent molecules of the first solvation sphere is not sufficient in order to reproduce the behavior of the electronic g-tensor in protic solvents, and that better agreement with experimental data can be obtained by including the long-range electrostatic effects accounted for by the PCM approach on top of the explicit hydrogen-bonded complexes.

Journal Article↗
PubMed2004

Electric field effects on the shielding constants of noble gases: a four-component relativistic Hartree-Fock study.

Second derivatives of nuclear shielding constants with respect to an electric field, i.e., shielding polarizabilities, have been calculated for the noble gas atoms from helium to xenon. The calculations have been carried out using the four-component relativistic Hartree-Fock method. In order to assess the importance of the individual relativistic corrections, the shielding polarizabilities have also been calculated at the nonrelativistic Hartree-Fock level, with spin-orbit and scalar (Darwin and mass-velocity) effects having been established by perturbative methods. Electron correlation effects have been estimated using the second-order polarization propagator approach. The relativistic effects on the tensor components of the shielding polarizabilities are found to be larger and changing less regularly with the atomic number than for the shielding constant itself. However, there is a partial cancellation of the contributions to the parallel and perpendicular components of the shielding polarizability and as a consequence the mean shielding polarizability is far less affected than the individual components.

Journal Article↗
PubMed2004

Perturbational relativistic theory of electron spin resonance g-tensor.

We carry out a complete treatment of the leading-order relativistic one-electron contributions, arising from the Breit-Pauli Hamiltonian, to the g-tensor of electron spin resonance spectroscopy. We classify the different terms and discuss their interpretation as well as give numerical ab initio estimates for the F2(-), Cl2(-), Br2(-), and I2(-) series, using analytical response theory calculations with a multiconfigurational self-consistent field reference state. The results are compared to available experimental data. (c) 2004 American Institute of Physics

Journal Article↗
PubMed2004

Atomic dipole moments calculated using analytical molecular second-moment gradients.

We have implemented analytical second-moment gradients for Hartree-Fock and multiconfigurational self-consistent-field wave functions. The code is used to calculate atomic dipole moments based on the generalized atomic polar tensor (GAPT) formalism [Phys. Rev. Lett. 62, 1469 (1989)], and the proposal of Dinur and Hagler (DH) for the calculation of atomic multipoles [J. Chem. Phys. 91, 2949 (1989)]. Both approaches display smooth basis-set convergence toward a well-defined basis-set limit and give reasonable electron correlation effects on the calculated atomic properties. However, the atomic charges and atomic dipole moments obtained from the GAPT partitioning scheme are unable to provide even qualitatively meaningful molecular quadrupole moments for some molecules, and thus the atomic multipole moments calculated in this scheme cannot be considered well suited for analyzing the electron density in molecules and for calculating intermolecular interaction energies. In contrast, the DH approach gives atomic charges and dipole moments that by definition exactly reproduce the molecular quadrupole moments. The approach of DH is, however, restricted to planar molecules and thus suffers from not being applicable to molecules of arbitrary shape. Both the GAPT and DH approaches give rather poor results for octupole and hexadecapole moments, indicating that at least atomic quadrupole moments are required for an accurate representation of the molecular charge distribution in terms of atomic electric moments.

Journal Article↗
PubMed2004

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.

Journal Article↗