PubMed Health⌕ Search

Biomedical subjects

Juha Vaara

Publications and source records attributed to Juha Vaara.

16 recordsLinked to original sources
PubMed2007

Nuclear magnetic resonance chemical shifts and quadrupole couplings for different hydrogen-bonding cases occurring in liquid water: a computational study.

Nuclear magnetic resonance (NMR) parameters are determined theoretically for the oxygen and hydrogen/deuterium nuclei of differently hydrogen-bonded water molecules in liquid water at 300 K. The parameters are the chemical shift, the shielding anisotropy, the asymmetry parameter of shielding, the nuclear quadrupole coupling constant, and the asymmetry parameter of the nuclear quadrupole coupling. We sample instantaneous configurations from a Car-Parrinello molecular dynamics simulation and feed nuclear coordinates into a quantum chemical program for the calculation of NMR parameters using density-functional theory with the three-parameter hybrid exchange-correlation (B3LYP) functional. In the subsequent analysis, molecules are divided into groups according to the number of hydrogen bonds they possess, and the full average NMR tensors are calculated separately for each group. The classification of the hydrogen-bonding cases is performed using a simple distance-based criterion. The analysis reveals in detail how the NMR tensors evolve as the environment changes gradually from gas to liquid upon increasing the number of hydrogen bonds to the molecule of interest. Liquid-state distributions of the instantaneous values of the NMR properties show a wide range of values for each hydrogen-bonding species with significant overlap between the different cases. Our study shows how local changes in the environment, along with classical thermal averaging, affect the NMR parameters in liquid water. For example, a broken or alternatively extra hydrogen bond induces major changes in the NMR tensors, and the effect is more pronounced for hydrogen or deuterium than for oxygen. The data sheds light on the usefulness of NMR experiments in investigating the local coordination of liquid water.

Journal Article↗
PubMed2006

Experimental and quantum-chemical determination of the (2)H quadrupole coupling tensor in deuterated benzenes.

Deuterium Quadrupole Coupling Constant (DQCC) in benzene was determined both experimentally by Nuclear Magnetic Resonance spectroscopy in Liquid Crystalline solutions (LC NMR) and theoretically by ab initio electronic structure calculations. DQCCs were measured for benzene-d(1) and 1,3,5-benzene-d(3) using several different liquid crystalline solvents and taking vibrational and deformational corrections into account in the analysis of experimental dipolar couplings, used to determine the orientational order parameter of the dissolved benzene. The experimental DQCC results for the isotopomers benzene-d(1) and 1,3,5-benzene-d(3) are found to be 187.7 kHz and 187.3 kHz, respectively, which are essentially equal within the experimental accuracy (+/-0.4 kHz). Theoretical results were obtained at different C-D bond lengths, and by applying corrections for electron correlation and rovibrational motion on top of large-basis-set Hartree-Fock results. The computations give a consistent DQCC of ca. 189 kHz for three different isotopomers; benzene-d(1), 1,3,5-benzene-d(3), and benzene-d(6), revealing that isotope effects are not detectable within the present experimental accuracy. Calculations carried out using a continuum solvation model to account for intermolecular interaction effects result in very small changes as compared to the data obtained in vacuo. The comparison of theoretical and experimental results points out the selection of the underlying molecular geometry as the most likely source of the remaining discrepancy of less than 2 kHz. Such an agreement between the calculated and the experimental DQCC results can only be achieved if rovibrational effects are considered on one hand in the experimental direct dipolar coupling data, and on the other hand in the theoretical property calculation, as is done presently.

Benzene↗
PubMed2006

Relativistic heavy-atom effects on heavy-atom nuclear shieldings.

The principal relativistic heavy-atom effects on the nuclear magnetic resonance (NMR) shielding tensor of the heavy atom itself (HAHA effects) are calculated using ab initio methods at the level of the Breit-Pauli Hamiltonian. This is the first systematic study of the main HAHA effects on nuclear shielding and chemical shift by perturbational relativistic approach. The dependence of the HAHA effects on the chemical environment of the heavy atom is investigated for the closed-shell X(2+), X(4+), XH(2), and XH(3) (-) (X=Si-Pb) as well as X(3+), XH(3), and XF(3) (X=P-Bi) systems. Fully relativistic Dirac-Hartree-Fock calculations are carried out for comparison. It is necessary in the Breit-Pauli approach to include the second-order magnetic-field-dependent spin-orbit (SO) shielding contribution as it is the larger SO term in XH(3) (-), XH(3), and XF(3), and is equally large in XH(2) as the conventional, third-order field-independent spin-orbit contribution. Considering the chemical shift, the third-order SO mechanism contributes two-thirds of the difference of approximately 1500 ppm between BiH(3) and BiF(3). The second-order SO mechanism and the numerically largest relativistic effect, which arises from the cross-term contribution of the Fermi contact hyperfine interaction and the relativistically modified spin-Zeeman interaction (FC/SZ-KE), are isotropic and practically independent of electron correlation effects as well as the chemical environment of the heavy atom. The third-order SO terms depend on these factors and contribute both to heavy-atom shielding anisotropy and NMR chemical shifts. While a qualitative picture of heavy-atom chemical shifts is already obtained at the nonrelativistic level of theory, reliable shifts may be expected after including the third-order SO contributions only, especially when calculations are carried out at correlated level. The FC/SZ-KE contribution to shielding is almost completely produced in the s orbitals of the heavy atom, with values diminishing with the principal quantum number. The relative contributions converge to universal fractions for the core and subvalence ns shells. The valence shell contribution is negligible, which explains the HAHA characteristics of the FC/SZ-KE term. Although the nonrelativistic theory gives correct chemical shift trends in present systems, the third-order SO-I terms are necessary for more reliable predictions. All of the presently considered relativistic corrections provide significant HAHA contributions to absolute shielding in heavy atoms.

Journal Article↗
PubMed2006

Density functional calculations of 3He chemical shift in endohedral helium fullerenes: Neutral, anionic, and di-helium species.

We report density functional calculations of 3He nuclear magnetic resonance chemical shifts in a series of experimentally known endohedral helium fullerenes, He(n)@Cm(q) (n = 1, 2; m = 60, 70, 76, 78; q = 0, 6-), including for the first time anionic and di-helium species. Despite the lack of dispersion in the density functional model, the results are in promising agreement with experiment. Density functional theory performs better than Hartree-Fock for the anionic systems. In the di-helium species confined in the small C60 cage, besides the atomic displacements from the center position, the direct He-He interactions contribute to the 3He shift.

Journal Article↗
PubMed2006

Calculations of nuclear quadrupole coupling in noble gas-noble metal fluorides: interplay of relativistic and electron correlation effects.

The nuclear quadrupole coupling constants (NQCCs) of noble gas and noble metal nuclei in the recently found noble gas-noble metal fluorides (NgMF, where Ng=Ar,Kr,Xe and M=Cu,Ag,Au) are obtained theoretically by high-level ab initio calculations, where both relativistic and electron correlation effects are included, and compared to experimental results. Fully relativistic four-component Dirac-Hartree-Fock (DHF) calculations are carried out at the basis set limit for electric field gradient that couples with the electric quadrupole moment of the nucleus, and uncorrelated relativistic effects are extracted by comparing DHF results to nonrelativistic (NR) HF calculations. Electron correlation effects are investigated both at fully relativistic second-order Moller-Plesset (DMP2) and at NR MP2 levels of theory, as well as at the NR coupled-cluster singles and doubles with perturbational triples [CCSD(T)] level. The validity of the approximation where relativistic effects, on the one hand, and nonrelativistically obtained correlation effects, on the other hand, are evaluated separately and assumed to be additive, is investigated by comparison with the DMP2 results. Inclusion of relativistic effects is shown to be necessary for obtaining the correct NQCC trends as the nucleus of interest and/or its neighbors become heavier. Electron correlation treatment is needed for approaching quantitative agreement with the experimental NQCCs. The assumption of additive electron correlation and relativistic effects, corresponding to the NR correlation treatment added on top of relativistic DHF data, gives qualitatively correct noble gas NQCCs. For noble metal NQCCs, correlation treatment at the relativistic level of theory is mandatory for reaching agreement with experimental results. Current work also confirms the experimental trends of NQCCs, which have been taken as an indication of nearly covalent interaction between noble gas and noble metal in the heaviest present systems, especially in XeAuF.

Journal Article↗
PubMed2006

Systematic Gaussian basis-set limit using completeness-optimized primitive sets. A case for magnetic properties.

We discuss the connection between the completeness of a basis set, measured by the completeness profiles introduced by Chong (Can J Chem 1995, 73, 79) at a certain exponent interval, and the possibility of reproducing molecular properties that arise either in the region close to the atomic nuclei or in the valence region. We present a scheme for generating completeness-optimized Gaussian basis sets, in which a preselected range of exponents is covered to an arbitrary accuracy. This is done by requiring Gaussian functions, the exponents of which are selected without reference to the atomic structure, to span the range with completeness profile as close to unity as wanted with as few functions as possible. The initial exponent range can be chosen suitable for calculations of molecular energetics or other valence-like properties. By extending the exponent range, properties requiring augmentation of the basis at a given angular momentum value and/or in a given distance range from the nucleus may be straightforwardly and systematically treated. In this scheme a universal, element-independent exponent set is generated in an automated way. The relation of basis-set completeness and performance in the calculation of magnetizability, nuclear magnetic shielding, and spin-spin coupling is tested with the completeness-optimized primitive sets and literature basis sets.

Algorithms↗
PubMed2005

Density-functional calculations of relativistic spin-orbit effects on nuclear magnetic shielding in paramagnetic molecules.

Terms arising from the relativistic spin-orbit effect on both hyperfine and Zeeman interactions are introduced to density-functional theory calculation of nuclear magnetic shielding in paramagnetic molecules. The theory is a generalization of the former nonrelativistic formulation for doublet systems and is consistent to O(alpha4), the fourth power of the fine structure constant, for the spin-orbit terms. The new temperature-dependent terms arise from the deviation of the electronic g tensor from the free-electron g value as well as spin-orbit corrections to hyperfine coupling tensor A, the latter introduced in the present work. In particular, the new contributions include a redefined isotropic pseudocontact contribution that consists of effects due to both the g tensor and spin-orbit corrections to hyperfine coupling. The implementation of the spin-orbit terms makes use of all-electron atomic mean-field operators and/or spin-orbit pseudopotentials. Sample results are given for group-9 metallocenes and a nitroxide radical. The new O(alpha4) corrections are found significant for the metallocene systems while they obtain small values for the nitroxide radical. For the isotropic shifts, none of the three beyond-leading-order hyperfine contributions are negligible.

Journal Article↗
PubMed2005

Perturbational calculations of parity-violating effects in nuclear-magnetic-resonance parameters.

We investigate the effects of the parity-violating electroweak interaction in the spectral parameters of nuclear magnetic resonance. Perturbational theory of parity-violating effects in the nuclear magnetic shielding is presented to the order of G(F)alpha, and in the indirect spin-spin coupling, to the order of G(F)alpha3. These leading-order parity-violating corrections are evaluated using analytical linear-response theory methods based on Hartree-Fock and density-functional theory reference states. Parity-violating contributions to spin-spin couplings are evaluated for the first time at the first-principles level. Calculations are carried out for two chiral halomethanes, bromochlorofluoromethane and bromofluoroiodomethane.

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↗
PubMed2004

Calculation of binary magnetic properties and potential energy curve in xenon dimer: second virial coefficient of (129)Xe nuclear shielding.

Quantum chemical calculations of the nuclear shielding tensor, the nuclear quadrupole coupling tensor, and the spin-rotation tensor are reported for the Xe dimer using ab initio quantum chemical methods. The binary chemical shift delta, the anisotropy of the shielding tensor Delta sigma, the nuclear quadrupole coupling tensor component along the internuclear axis chi( parallel ), and the spin-rotation constant C( perpendicular ) are presented as a function of internuclear distance. The basis set superposition error is approximately corrected for by using the counterpoise correction (CP) method. Electron correlation effects are systematically studied via the Hartree-Fock, complete active space self-consistent field, second-order Møller-Plesset many-body perturbation, and coupled-cluster singles and doubles (CCSD) theories, the last one without and with noniterative triples, at the nonrelativistic all-electron level. We also report a high-quality theoretical interatomic potential for the Xe dimer, gained using the relativistic effective potential/core polarization potential scheme. These calculations used valence basis set of cc-pVQZ quality supplemented with a set of midbond functions. The second virial coefficient of Xe nuclear shielding, which is probably the experimentally best-characterized intermolecular interaction effect in nuclear magnetic resonance spectroscopy, is computed as a function of temperature, and compared to experiment and earlier theoretical results. The best results for the second virial coefficient, obtained using the CCSD(CP) binary chemical shift curve and either our best theoretical potential or the empirical potentials from the literature, are in good agreement with experiment. Zero-point vibrational corrections of delta, Delta sigma, chi (parallel), and C (perpendicular) in the nu=0, J=0 rovibrational ground state of the xenon dimer are also reported.

Journal Article↗
PubMed2004

Nuclear magnetic shielding and quadrupole coupling tensors in liquid water: a combined molecular dynamics simulation and quantum chemical study.

Nuclear magnetic resonance (NMR) shielding tensors for the oxygen and hydrogen nuclei, as well as nuclear quadrupole coupling tensors for the oxygen and deuterium nuclei of water in the liquid and gaseous state, are calculated using Hartree-Fock and density functional theory methods, for snapshots sampled from Car-Parrinello molecular dynamics trajectories. Clusters representing local liquid structures and instantaneous configurations of a single molecule representing low-density gas are fed into a quantum chemical program for the calculation of the NMR tensors. The average isotropic and anisotropic tensorial properties of 400 samples in both states, averaged using a common Eckart coordinate frame, are calculated from the data. We report results for the gas-to-liquid chemical shifts of (17)O and (1)H nuclei, as well as the corresponding change in the nuclear quadrupole couplings of (17)O and (2)H. Full thermally averaged shielding and quadrupole coupling tensors are reported for the gaseous and liquid-state water, for the first time in the case of liquid. Electron correlation effects, the difference of classical vs quantum mechanical rovibrational averaging, and different methods of averaging anisotropic properties are discussed.

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

Relativistic spin-orbit effects on hyperfine coupling tensors by density-functional theory.

A second-order perturbation theory treatment of spin-orbit corrections to hyperfine coupling tensors has been implemented within a density-functional framework. The method uses the all-electron atomic mean-field approximation and/or spin-orbit pseudopotentials in incorporating one- and two-electron spin-orbit interaction within a first-principles framework. Validation of the approach on a set of main-group radicals and transition metal complexes indicates good agreement between all-electron and pseudopotential results for hyperfine coupling constants of the lighter nuclei in the system, except for cases in which scalar relativistic effects become important. The nonrelativistic Fermi contact part of the isotropic hyperfine coupling constants is not always accurately reproduced by the exchange-correlation functionals employed, particularly for the triplet and pi-type doublet radicals in the present work. For this reason, ab initio coupled-cluster singles and doubles with perturbative triples results for the first-order contributions have been combined in the validation calculations with the density-functional results for the second-order spin-orbit contributions. In the cases where spin-orbit corrections are of significant magnitude relative to the nonrelativistic first-order terms, they improve the agreement with experiment. Antisymmetric contributions to the hyperfine tensor arise from the spin-orbit contributions and are discussed for the IO2 radical, whereas rovibrational effects have been evaluated for RhC, NBr, and NI.

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

Relativistic spin-orbit coupling effects on secondary isotope shifts of (13)C nuclear shielding in CX(2) (X = O, S, Se, Te).

Rovibrational corrections, temperature dependence, and secondary isotope shifts of the (13)C nuclear shielding in CX(2) (X = O, S, Se, Te) are calculated taking into account the relativistic spin-orbit (SO) interaction. The SO effect is considered for the first time for the secondary isotope shifts. The nuclear shielding hypersurface in terms of nuclear displacements is calculated by using a density-functional theory method. Ab initio multiconfiguration self-consistent field calculations are done at the equilibrium geometry for comparison. (13)C NMR measurements are carried out for CS(2). The calculated results are compared with both present and earlier experimental data on CO(2), CS(2), and CSe(2). The heavy-atom SO effects on the rovibrational corrections of (13)C shielding are shown to be significant. For CSe(2) and CTe(2), reliable prediction of secondary isotope effects and their temperature dependence requires the inclusion of the SO corrections. In particular, earlier discrepancies of theory and experiment for CSe(2) are fully resolved by taking the SO interactions into account.

Journal Article↗