PubMed Health⌕ Search

Biomedical subjects

Vladimir G Malkin

Publications and source records attributed to Vladimir G Malkin.

9 recordsLinked to original sources
PubMed2007

Density functional calculations of NMR shielding tensors for paramagnetic systems with arbitrary spin multiplicity: validation on 3d metallocenes.

The calculation of nuclear shieldings for paramagnetic molecules has been implemented in the ReSpect program, which allows the use of modern density functional methods with accurate treatments of spin-orbit effects for all relevant terms up to order Omicron(alpha4) in the fine structure constant. Compared to previous implementations, the methodology has been extended to compounds of arbitrary spin multiplicity. Effects of zero-field splittings in high-spin systems are approximately accounted for. Validation of the new implementation is carried out for the 13C and 1H NMR signal shifts of the 3d metallocenes 4VCp2, 3CrCp2, 2MnCp2, 6MnCp2, 2CoCp2, and 3NiCp2. Zero-field splitting effects on isotropic shifts tend to be small or negligible. Agreement with experimental isotropic shifts is already good with the BP86 gradient-corrected functional and is further improved by admixture of Hartree-Fock exchange in hybrid functionals. Decomposition of the shieldings confirms the dominant importance of the Fermi-contact shifts, but contributions from spin-orbit dependent terms are frequently also non-negligible. Agreement with 13C NMR shift tensors from solid-state experiments is of similar quality as for isotropic shifts.

Journal Article↗
PubMed2006

Calculation of zero-field splitting parameters: comparison of a two-component noncolinear spin-density-functional method and a one-component perturbational approach.

Two different sets of approaches for the density-functional calculation of the spin-orbit contributions to zero-field splitting (ZFS) parameters of high-spin systems have been implemented within the same quantum chemistry code ReSpect and have been validated and compared for a series of model systems. The first approach includes spin-orbit coupling variationally in a two-component calculation, using either an all-electron Douglas-Kroll-Hess ansatz or two-component relativistic pseudopotentials. The ZFS parameters are computed directly from energy differences between different relativistic states. Additionally, an approximate second-order perturbation theory approach has been implemented, based on nonrelativistic or scalar relativistic wave functions. For a series of group 16 triplet diatomics and for the octet GdH3 molecules, two-component density functional calculations underestimate the zero-field splitting D systematically by a factor of 2. This may be rationalized readily by the incomplete description of states with absolute value MJ < J by a single-determinantal wave function built from two-component spinors. In the case of two 3d transition metal complexes and for GdH3, the results depend furthermore sensitively on exchange-correlation functional. Results of the alternative one-component approach agree strikingly with the two-component data for systems with small spin-orbit effects and start to deviate from them only for heavier systems with large spin-orbit effects. These results have fundamental implications for the achievable accuracy of one-component density-functional approaches used widely to compute ZFS parameters in the field of molecular magnetism. Possible refinements of both one-and two-component approaches are discussed.

Journal Article↗
PubMed2006

Scalar relativistic calculations of hyperfine coupling tensors using the Douglas-Kroll-Hess method with a finite-size nucleus model.

A scalar relativistic method to calculate hyperfine coupling tensors at the Douglas-Kroll-Hess level has been extended to incorporate a finite-size nucleus model using a Gaussian charge and magnetic moment distribution. Density functional calculations at gradient-corrected and hybrid functional levels have been carried out for the group 11 atoms and for a set of small group 12 molecules, comparing nonrelativistic as well as scalar relativistic results at second-order Douglas-Kroll-Hess level with and without finite-size nucleus. While nonrelativistic calculations underestimate isotropic hyperfine couplings increasingly with increasing nuclear charge, scalar relativistic calculations with point nucleus provide somewhat overestimated values. Inclusion of the finite-size nuclear model in the calculation of the wavefunction, and in the transformed hyperfine operators both decrease the magnitude of the hyperfine couplings. The effects, which are cumulative, improve agreement with experiment.

Journal Article↗
PubMed2006

Resolution of identity Dirac-Kohn-Sham method using the large component only: Calculations of g-tensor and hyperfine tensor.

A new relativistic two-component density functional approach, based on the Dirac-Kohn-Sham method and an extensive use of the technique of resolution of identity (RI), has been developed and is termed the DKS2-RI method. It has been applied to relativistic calculations of g and hyperfine tensors of coinage-metal atoms and some mercury complexes. The DKS2-RI method solves the Dirac-Kohn-Sham equations in a two-component framework using explicitly a basis for the large component only, but it retains all contributions coming from the small component. The DKS2-RI results converge to those of the four-component Dirac-Kohn-Sham with an increasing basis set since the error associated with the use of RI will approach zero. The RI approximation provides a basis for a very efficient implementation by avoiding problems associated with complicated integrals otherwise arising from the elimination of the small component. The approach has been implemented in an unrestricted noncollinear two-component density functional framework. DKS2-RI is related to Dyall's [J. Chem. Phys. 106, 9618 (1997)] unnormalized elimination of the small component method (which was formulated at the Hartree-Fock level and applied to one-electron systems only), but it takes advantage of the local Kohn-Sham exchange-correlation operators (as, e.g., arising from local or gradient-corrected functionals). The DKS2-RI method provides an attractive alternative to existing approximate two-component methods with transformed Hamiltonians (such as Douglas-Kroll-Hess [Ann. Phys. 82, 89 (1974); Phys. Rev. A 33, 3742 (1986)] method, zero-order regular approximation, or related approaches) for relativistic calculations of the structure and properties of heavy-atom systems. In particular, no picture-change effects arise in the property calculations.

Journal Article↗
PubMed2005

Relativistic two-component calculations of electronic g-tensors that include spin polarization.

The first two-component relativistic density-functional approach for the calculation of electronic g-tensors is reported that includes spin polarization using noncollinear spin-density functionals. The method is based on the relativistic Douglas-Kroll-Hess Hamiltonian and has been implemented into the ReSpect program package. Using three self-consistent-field calculations with orthogonal orientations of total magnetization J, the full g-matrix may be obtained. In contrast to previous spin-restricted two-component treatments, results with the new approach agree excellently with spin-polarized one-component calculations for light-atom radicals. Additionally, unlike one-component approaches, the method also reproduces successfully the negative deltag(parallel)-values of heavy-atom 2sigma radicals and the negative deltag(perpendicular) components in cysteinyl. The new method removes effectively the dilemma existing up to now regarding the simultaneous inclusion of spin polarization and higher-order spin-orbit effects in g-tensor calculations. It is straightforwardly applicable to higher than doublet spin multiplicities and has been implemented with hybrid functionals.

Journal Article↗
PubMed2005

Quantum chemical investigation of hyperfine coupling constants on first coordination sphere water molecule of gadolinium(III) aqua complexes.

Hyperfine interactions (HFI) on the nuclei of the first coordination sphere water molecules in a model [Gd(H(2)O)(8)](3+) aqua complex and in the magnetic resonance imaging contrast agent [Gd(DOTA)(H(2)O)](-) were studied theoretically. Density functional theory (DFT) calculations combined with classical molecular dynamics (MD) simulations have been used in order to take into account dynamic effects in aqueous solution. DFT relativistic calculations show a strong spin-polarization of the first coordination sphere water molecules. This spin-polarization leads to a positive (17)O isotropic hyperfine coupling constant (A(iso)((17)O) = 0.58 +/- 0.11 MHz) and to a significant increase of the effective distance (r(eff)(Gd-O) = 2.72 +/- 0.06 A) of dipolar interaction compared to the mean internuclear distance (r(Gd-O) = 2.56 +/- 0.06 A) obtained from the MD trajectory of [Gd(DOTA)(H(2)O)](-) in aqueous solution. The point-dipole model for anisotropic hyperfine interaction overestimates therefore the longitudinal relaxation rate of the (17)O nucleus by approximately 45%. The (1)H isotropic hyperfine coupling constant of the bound water molecule is predicted to be very small (A(iso)((1)H) = 0.03 +/- 0.02 MHz), and the point-dipole approximation for first coordination sphere water protons holds. The calculated hyperfine parameters are in good agreement with available experimental data.

Journal Article↗
PubMed2002

Density functional calculations of electronic g-tensors for semiquinone radical anions. The role of hydrogen bonding and substituent effects.

A recently developed density functional approach has been used to carry out a systematic computational study of electronic g-tensors for a series of 1,4-semiquinone radical anions. Good agreement with high-field EPR data in frozen 2-propanol is achieved only after taking into account the significant reduction of g-tensor anisotropy caused by hydrogen bonding to solvent molecules. The comparison of various model systems for the first solvation shell suggests two hydrogen bonds from 2-propanol molecules to each of the carbonyl groups of the radical anions, and one additional hydrogen bond to each of the methoxy groups in ubiquinone systems. 2-Propanol makes stronger hydrogen bonds than water and thus influences g-tensor anisotropy more strongly. Substituent effects at the semiquinone are reproduced quantitatively by the calculations. The g-tensor anisotropy is influenced significantly by the conformations of methyl and methoxy substituents, with opposite contributions. Analyses and interpretations of the interrelations between structure, bonding, and spectroscopic data are provided. The relevance of the computational results for the EPR spectroscopy of semiquinone radical anions in photosynthetic reaction centers is discussed.

1-Propanol↗
PubMed2002

Calculation of electronic g-tensors for transition metal complexes using hybrid density functionals and atomic meanfield spin-orbit operators.

We report the first implementation of the calculation of electronic g-tensors by density functional methods with hybrid functionals. Spin-orbit coupling is treated by the atomic meanfield approximation. g-Tensors for a set of small main group radicals and for a series of ten 3d and two 4d transition metal complexes have been compared using the local density approximation (VWN functional), the generalized gradient approximation (BP86 functional), as well as B3-type (B3PW91) and BH-type (BHPW91) hybrid functionals. For main group radicals, the effect of exact-exchange mixing is small. In contrast, significant differences between the various functionals arise for transition metal complexes. As has been shown previously, local and in particular gradient-corrected functionals tend to underestimate the "paramagnetic" contributions to the g-tensors in these cases and thereby recover only about 40-50% of the range of experimental g-tensor components. This is improved to ca. 60% by the B3PW91 functional, which also gives slightly reduced standard deviations. The range increases to almost 100% using the half-and-half functional BHPW91. However, the quality of the correlation with experimental data worsens due to a significant overestimate of some intermediate g-tensor values. The worse performance of the BHPW91 functional in these cases is accompanied by spin contamination. Although none of the functionals tested thus appears to be ideal for the treatment of electronic g-tensors in transition metal complexes, the B3PW91 hybrid functional exhibited the overall most satisfactory performance. Apart from the validation of hybrid functionals, some aspects in the treatment of spin-orbit contributions to the g-tensor are discussed.

Journal Article↗