Hans Jørgen Aagaard Jensen

Relativistic 4-component calculations of indirect nuclear spin-spin couplings in MH4(M=C,Si,Ge,Sn,Pb (CH3)3H

Relativistic four-component random phase approximation (RPA) calculations of indirect nuclear spin–spin coupling constants in MH4 (M=C, Si, Ge, Sn, Pb) and Pb(CH3)3H are presented. The need for tight s-functions also in relativistic four-component calculations is verified and explained, and the effect of omission of (SS–LL) and (SS–SS) two-electron integrals is investigated. Already in GeH4 we see a relativistic increase in the coupling constant by 12%, and for PbH4 the effect is a 156% increase for the one-bond coupling. Large relativistic effects are also computed for the two-bonds couplings. We find that the relativistic effects on the one-bond couplings are mainly due to scalar relativistic factors rather than spin–orbit corrections.

Large-scale parallel configuration interaction. II. Two- and four-component double-group general active space implementation with application to BiH.

We present a parallel implementation of a large-scale relativistic double-group configuration interaction (CI) program. It is applicable with a large variety of two- and four-component Hamiltonians. The parallel algorithm is based on a distributed data model in combination with a static load balancing scheme. The excellent scalability of our parallelization scheme is demonstrated in large-scale four-component multireference CI (MRCI) benchmark tests on two of the most common computer architectures, and we also discuss hardware-dependent aspects with respect to possible speedup limitations. With the new code we have been able to calculate accurate spectroscopic properties for the ground state and the first excited state of the BiH molecule using extensive basis sets. We focused, in particular, on an accurate description of the splitting of these two states which is caused by spin-orbit coupling. Our largest parallel MRCI calculation thereby comprised an expansion length of 2.7x10(9) Slater determinants.

Correlated calculations of indirect nuclear spin-spin coupling constants for XH4 (X = Si, Ge, and Sn)

Abstract We report multiconfigurational linear response calculations of the one-bond and geminal indirect nuclear spin-spin couplings constants in SiH4, GeH4, and SnH4. Results of a basis set study on SnH4 using 13 different basis sets on Sn are given, showing the effects on the individual contributions to the coupling constant. In addition to the more conventional choice of active orbital space (one correlating orbital for each strongly occupied valence orbital), calculations employing extended active spaces were carried out. The results are compared with available experimental data showing that correlation effects are very important. For the one-bond couplings good agreement with experiment, within 1–5%, is obtained, whereas the geminai proton-proton couplings show larger deviations. The calculations reported here are, to our knowledge, the first that take into consideration correlation effects in all four contributions to the coupling constants for molecules as large as GeH4 and SnH4. All X-H couplings are dominated by the Fermi contact term while the picture is more diverse for the geminal couplings.

A multiconfigurational hybrid density-functional theory

We propose a multiconfigurational hybrid density-functional theory which rigorously combines a multiconfiguration self-consistent-field calculation with a density-functional approximation based on a linear decomposition of the electron-electron interaction. This gives a straightforward extension of the usual hybrid approximations by essentially adding a fraction λ of exact static correlation in addition to the fraction λ of exact exchange. Test calculations on the cycloaddition reactions of ozone with ethylene or acetylene and the dissociation of diatomic molecules with the Perdew-Burke-Ernzerhof and Becke-Lee-Yang-Parr density functionals show that a good value of λ is 0.25, as in the usual hybrid approximations. The results suggest that the proposed multiconfigurational hybrid approximations can improve over usual density-functional calculations for situations with strong static correlation effects.

Gauge origin independent calculations of nuclear magnetic shieldings in relativistic four-component theory

The use of perturbation-dependent London atomic orbitals, also called gauge including atomic orbitals, has proven efficient for calculations of NMR shielding constants and other magnetic properties in the nonrelativistic framework. In this paper, the theory of London atomic orbitals for NMR shieldings is extended to the four-component relativistic framework and our implementation is described. The relevance of London atomic orbitals in four-component calculations as well as computational aspects are illustrated with test calculations on hydrogen iodide. We find that the use of London atomic orbitals is an efficient method for reliable calculations of NMR shielding constants with standard basis sets, also for four-component calculations with spin-orbit coupling effects included in the wave function optimization. Furthermore, we find that it is important that the small component basis functions fulfill the magnetic balance for accurate description of the diamagnetic shielding and that the role of London atomic orbitals in the relativistic domain is to provide atomic magnetic balance even in the molecular case, thus greatly improving basis set convergence. The Sternheim approximation, which calculates the diamagnetic contribution as an expectation value, leads to significant errors and is not recommended.

Relativistic effects in the intermolecular interaction-induced nuclear magnetic resonance parameters of xenon dimer

Relativistic effects on the (129)Xe nuclear magnetic resonance shielding and (131)Xe nuclear quadrupole coupling (NQC) tensors are examined in the weakly bound Xe(2) system at different levels of theory including the relativistic four-component Dirac-Hartree-Fock (DHF) method. The intermolecular interaction-induced binary chemical shift delta, the anisotropy of the shielding tensor Deltasigma, and the NQC constant along the internuclear axis chi( parallel) are calculated as a function of the internuclear distance. DHF shielding calculations are carried out using gauge-including atomic orbitals. For comparison, the full leading-order one-electron Breit-Pauli perturbation theory (BPPT) is applied using a common gauge origin. Electron correlation effects are studied at the nonrelativistic (NR) coupled-cluster singles and doubles with perturbational triples [CCSD(T)] level of theory. The fully relativistic second-order Moller-Plesset many-body perturbation (DMP2) theory is used to examine the cross coupling between correlation and relativity on NQC. The same is investigated for delta and Deltasigma by BPPT with a density functional theory model. A semiquantitative agreement between the BPPT and DHF binary property curves is obtained for delta and Deltasigma in Xe(2). For these properties, the currently most complete theoretical description is obtained by a piecewise approximation where the uncorrelated relativistic DHF results obtained close to the basis-set limit are corrected, on the one hand, for NR correlation effects and, on the other hand, for the BPPT-based cross coupling of relativity and correlation. For chi( parallel), the fully relativistic DMP2 results obtain a correction for NR correlation effects beyond MP2. The computed temperature dependence of the second virial coefficient of the (129)Xe nuclear shielding is compared to experiment in Xe gas. Our best results, obtained with the piecewise approximation for the binary chemical shift combined with the previously published state of the art theoretical potential energy curve for Xe(2), are in excellent agreement with the experiment for the first time.

Molecular-Level Insight into the Spectral Tuning Mechanism of the DsRed Chromophore.

We present a detailed study of the protein environmental effects on the one- and two-photon absorption (1PA and 2PA, respectively) properties of the S0-S1 transition in the DsRed protein using the polarizable embedding density functional theory formalism. We find that steric factors and chromophore-protein interactions act in concert to enhance the 2PA activity inside the protein while adversely blue-shifting the 1PA maximum. A two-state model reveals that the 2PA intensity gain is primarily governed by the increased change in the permanent dipole moment between the ground and the excited states acquired inside the protein. Our results indicate that this mainly is attributable to counter-directional contributions stemming from Lys163 and the conserved Arg95 with the former additionally identified as a key residue in the color tuning mechanism. The results provide new insight into the tuning mechanism of DsRed and suggest a possible strategy for simultaneous improvement of its 1PA and 2PA properties.

Density matrix renormalization group with efficient dynamical electron correlation through range separation

We present a new hybrid multiconfigurational method based on the concept of range-separation that combines the density matrix renormalization group approach with density functional theory. This new method is designed for the simultaneous description of dynamical and static electron-correlation effects in multiconfigurational electronic structure problems.

The Multi-Configuration Self-Consistent Field Method Within a Polarizable Embedded Framework

We present a detailed derivation of Multi-Configuration Self-Consistent Field (MCSCF) optimization and linear response equations within the polarizable embedding scheme: PE-MCSCF. The MCSCF model enables a proper description of multiconfigurational effects in reaction paths, spin systems, excited states, and other properties which cannot be described adequately with current implementations of polarizable embedding in density functional or coupled cluster theories. In the PE-MCSCF scheme the environment surrounding the central quantum mechanical system is represented by distributed multipole moments and anisotropic dipole-dipole polarizabilities. The PE-MCSCF model has been implemented in DALTON. As a preliminary application, the low lying valence states of acetone and uracil in water has been calculated using Complete Active Space Self-Consistent Field (CASSCF) wave functions. The dynamics of the water environment have been simulated using a series of snapshots generated from classical Molecular Dynamics. The calculated shifts from gas-phase to water display between good and excellent correlation with experiment and previous calculations. As an illustration of another area of potential applications we present calculations of electronic transitions in the transition metal complex, [Fe(NO)(CN)5](2-) in a micro-solvated environment. This system is highly multiconfigurational and the influence of solvation is significant.

Assessment of charge-transfer excitations with time-dependent, range-separated density functional theory based on long-range MP2 and multiconfigurational self- consistent field wave functions

Charge transfer excitations can be described within Time-Dependent Density Functional Theory (TD-DFT), not only by means of the Coulomb Attenuated Method (CAM) but also with a combination of wave function theory and TD-DFT based on range separation. The latter approach enables a rigorous formulation of multi-determinantal TD-DFT schemes where excitation classes, which are absent in conventional TD-DFT spectra (like for example double excitations), can be addressed. This paper investigates the combination of both the long-range Multi-Configuration Self-Consistent Field (MCSCF) and Second Order Polarization Propagator Approximation (SOPPA) ansätze with a short-range DFT (srDFT) description. We find that the combinations of SOPPA or MCSCF with TD-DFT yield better results than could be expected from the pure wave function schemes. For the Time-Dependent MCSCF short-range DFT ansatz (TD-MC-srDFT) excitation energies calculated over a larger benchmark set of molecules with predominantly single reference character yield good agreement with their reference values, and are in general comparable to the CAM-B3LYP functional. The SOPPA-srDFT scheme is tested for a subset of molecules used for benchmarking TD-MC-srDFT and performs slightly better against the reference data for this small subset. Beyond the proof-of-principle calculations comprising the first part of this contribution, we additionally studied the low-lying singlet excited states (S1 and S2) of the retinal chromophore. The chromophore displays multireference character in the ground state and both excited states exhibit considerable double excitation character, which in turn cannot be described within standard TD-DFT, due to the adiabatic approximation. However, a TD-MC-srDFT approach can account for the multireference character, and excitation energies are obtained with accuracy comparable to CASPT2, although using a much smaller active space.