Methodological aspects in the calculation of parity-violating effects in nuclear magnetic resonance parameters

We examine the quantum chemical calculation of parity-violating (PV) electroweak contributions to the spectral parameters of nuclear magnetic resonance (NMR) from a methodological point of view. Nuclear magnetic shielding and indirect spin-spin coupling constants are considered and evaluated for three chiral molecules, H2O2, H2S2, and H2Se2. The effects of the choice of a one-particle basis set and the treatment of electron correlation, as well as the effects of special relativity, are studied. All of them are found to be relevant. The basis-set dependence is very pronounced, especially at the electron correlated ab initio levels of theory. Coupled-cluster and density-functional theory (DFT) results for PV contributions differ significantly from the Hartree-Fock data. DFT overestimates the PV effects, particularly with nonhybrid exchange-correlation functionals. Beginning from third-row elements, special relativity is of importance for the PV NMR properties, shown here by comparing perturbational one-component and various four-component calculations. In contrast to what is found for nuclear magnetic shielding, the choice of the model for nuclear charge distribution--point charge or extended (Gaussian)--has a significant impact on the PV contribution to the spin-spin coupling constants.     

An infinite-order two-component relativistic Hamiltonian by a simple one-step transformation

The authors report the implementation of a simple one-step method for obtaining an infinite-order two-component (IOTC) relativistic Hamiltonian using matrix algebra. They apply the IOTC Hamiltonian to calculations of excitation and ionization energies as well as electric and magnetic properties of the radon atom. The results are compared to corresponding calculations using identical basis sets and based on the four-component Dirac-Coulomb Hamiltonian as well as Douglas-Kroll-Hess and zeroth-order regular approximation Hamiltonians, all implemented in the DIRAC program package, thus allowing a comprehensive comparison of relativistic Hamiltonians within the finite basis approximation.     

Efficient synthesis of 4,7-diamino substituted 1,10-phenanthroline-2,9-dicarboxamides

A convenient and high-yielding multigram synthesis of the versatile intermediate 4,7-dichloro-1,10-phenanthroline-2,9-dicarboxylic acid is described. The intermediate is further efficiently derivatized to 4,7-diamino-1,10-phenanthroline-2,9-dicarboxamides with potential G-quadruplex stabilizing effects.     

NVU dynamics. I. Geodesic motion on the constant-potential-energy hypersurface

An algorithm is derived for computer simulation of geodesics on the constant-potential-energy hypersurface of a system of N classical particles. First, a basic time-reversible geodesic algorithm is derived by discretizing the geodesic stationarity condition and implementing the constant-potential-energy constraint via standard Lagrangian multipliers. The basic NVU algorithm is tested by single-precision computer simulations of the Lennard-Jones liquid. Excellent numerical stability is obtained if the force cutoff is smoothed and the two initial configurations have identical potential energy within machine precision. Nevertheless, just as for NVE algorithms, stabilizers are needed for very long runs in order to compensate for the accumulation of numerical errors that eventually lead to "entropic drift" of the potential energy towards higher values. A modification of the basic NVU algorithm is introduced that ensures potential-energy and step-length conservation; center-of-mass drift is also eliminated. Analytical arguments confirmed by simulations demonstrate that the modified NVU algorithm is absolutely stable. Finally, we present simulations showing that the NVU algorithm and the standard leap-frog NVE algorithm have identical radial distribution functions for the Lennard-Jones liquid.     

Connections on modules over singularities of finite CM representation type

Let A be a commutative k-algebra, where k is an algebraically closed field of characteristic 0, and let M be an A-module. We consider the following question: Under what conditions is it possible to find a connection on M?We consider the maximal Cohen-Macaulay (MCM) modules over complete CM algebras that are isolated singularities, and usually assume that the singularities have finite CM representation type. It is known that any MCM module over a simple singularity of dimension d≤2 admits an integrable connection. We prove that an MCM module over a simple singularity of dimension d≥3 admits a connection if and only if it is free. Among singularities of finite CM representation type, we find examples of curves with MCM modules that do not admit connections, and threefolds with non-free MCM modules that admit connections.Let A be a singularity not necessarily of finite CM representation type, and consider the condition that A is a Gorenstein curve or a -Gorenstein singularity of dimension d≥2. We show that this condition is sufficient for the canonical module ω_A to admit an integrable connection, and conjecture that it is also necessary. In support of the conjecture, we show that if A is a monomial curve singularity, then the canonical module ω_A admits an integrable connection if and only if A is Gorenstein.     

The frequency-dependent dipole polarizability of the mercury dimer from four-component relativistic density-functional theory

The frequency-dependent dipole polarizability of Hg(2) is calculated using response theory within four-component relativistic density-functional theory [using the local-density approximation (LDA) and the hybrid functional B3LYP] including corrections for the basis-set superposition error. The anisotropic component of the polarizability tensor agrees well with the values obtained from collision-induced Raman spectroscopy carried out at a wavelength of 488 nm. The values obtained from the two density functionals agree closely with the experimentally derived anisotropy component of the dipole polarizability, despite their rather large differences in the dimer potential-energy curves (LDA is strongly overbinding while B3LYP is purely repulsive). The first two refractivity virial coefficients for the generalized Clausius-Mossotti function are derived.