Optimization of atomic density-fitting basis functions for molecular two-electron integral approximations.

A general procedure for the optimization of atomic density-fitting basis functions is designed with the balance between accuracy and numerical stability in mind. Given one-electron wavefunctions and energies, weights are assigned to the product densities, modeling their contribution to the exchange and second-order correlation energy, and a simple weighted error measure is minimized. Generally contracted Gaussian auxiliary basis sets are optimized to match the wavefunction basis sets [D. N. Laikov, Theor. Chem. Acc. 138, 40 (2019)] for all 102 elements in a scalar-relativistic approximation [D. N. Laikov, J. Chem. Phys. 150, 061103 (2019)].

Control of Photoisomerization of an Azoazacryptand by Anion Binding and Cucurbit[8]uril Encapsulation in an Aqueous Solution.

Control of isomerization of a receptor bearing multiple light-switchable subunits in a confined space is critical for the design of synthetic molecular machines. Toward this goal, a new azacryptand containing three azobenzene subunits has been developed, and its photoisomerization in an aqueous solution has been studied depending on anion coordination and recognition by a larger host-cucurbit[8]uril (CB[8]). The cryptand in its hexaprotonated form shows considerable affinity for fluoride and perchlorate, which in turn affects the isomer distribution of the receptor under UV-light irradiation, stabilizing the isomers of the cryptand with Z-configurations. CB[8] was found to be able to encapsulate the isomers of the cryptand by forming a Matryoshka-type complex. The irradiation of a 10:1 CB[8]-cryptand mixture has led to a selective conversion of the cryptand to the E,E,Z isomer inside CB[8]. It has been demonstrated that the addition of fluoride to the resulted complex induces the release of the cryptand as a major E,E,E isomer, while other studied anions were ineffective in this reaction. To our knowledge, this work presents a first example of a host-controlled photoisomerization of an anion receptor bearing multiple switching azobenzenes that model the function of naturally occurring chaperones.

Synthesis and biological activity of 5-aryliden-2-thiohydantoin S-aryl derivatives.

Three new and complementary approaches to S-arylation of 2-thiohydantoins have been developed: copper-catalyzed cross coupling with either arylboronic acids or aryl iodides under mild conditions, or direct nucleophilic substitution in activated aryl halides. For 38 diverse compounds, reaction yields for all three methods have been determined. Selected by molecular docking, they have been tested on androgen receptor activation, and p53-Mdm2 regulation, and A549, MCF7, VA13, HEK293T, PC3, LnCAP cell lines for cytotoxicity, Two of them turned out to be promising as androgen receptor activators (likely by allosteric regulation), and another one is shown to activate the p53 cascade. It is hoped that 2-thiohydantoin S-arylidenes are worth further studies as biologically active compounds.

Ullmann-type C-Se Cross-Coupling in the Hydantoin Family: Synthesis, Mechanistic Studies, and Tests of Biological Activity.

An attractive strategy for C-Se bond formation by Ullmann-type copper(I)-promoted cross-coupling is developed. A wide range of aryliodides reacts with various disubstituted 2-selenohydantoins under mild conditions and provides Se-arylated imidazolines in moderate to high yields. Computational mechanistic studies show the oxidative addition/intramolecular reductive elimination likely to be the lowest-energy pathway. Cytotoxic activity of all 43 reaction products has been tested in vitro against MCF7 and A549 cancer cell lines with VA13 and MCF10a control cells.

Simple exchange hole models for long-range-corrected density functionals.

Density functionals with a range-separated treatment of the exchange energy are known to improve upon their semilocal forerunners and fixed-fraction hybrids. The conversion of a given semilocal functional into its short-range analog is not straightforward, however, and not even unique, because the latter has a higher information content that has to be recovered in some way. Simple models of the spherically averaged exchange hole as an interpolation between the uniform electron gas limit and a few-term Hermite function are developed here for use with generalized-gradient approximations, so that the energy density of the error-function-weighted Coulomb interaction is given by explicit closed-form expressions in terms of elementary and error functions. For comparison, some new nonoscillatory models in the spirit of earlier works are also built and studied; the energy densities from both kinds of models match rather closely (within less than 5%), but the latter slightly mismatches (by about 1%) the exact uniform electron gas limit.

Additive atomic approximation for relativistic effects: A two-component Hamiltonian for molecular electronic structure calculations.

An approximate relativistic two-component Hamiltonian for use in molecular electronic structure calculations is derived in the form of a sum of fixed atom-centered kinetic and spin-orbit operators added to the non-relativistic Hamiltonian. Starting from the well-known zeroth-order regular approximation, further steps are taken to get rid of its nonlinearity in the potential, ending up with a simple formulation with easily computable integrals that can seamlessly work with any traditional electronic structure method. Molecular tests show a good accuracy of this approximation.

A new parametrizable model of molecular electronic structure.

A new electronic structure model is developed in which the ground state energy of a molecular system is given by a Hartree-Fock-like expression with parametrized one- and two-electron integrals over an extended (minimal + polarization) set of orthogonalized atom-centered basis functions, the variational equations being solved formally within the minimal basis but the effect of polarization functions being included in the spirit of second-order perturbation theory. It is designed to yield good dipole polarizabilities and improved intermolecular potentials with dispersion terms. The molecular integrals include up to three-center one-electron and two-center two-electron terms, all in simple analytical forms. A method to extract the effective one-electron Hamiltonian of nonlocal-exchange Kohn-Sham theory from the coupled-cluster one-electron density matrix is designed and used to get its matrix representation in a molecule-intrinsic minimal basis as an input to the parametrization procedure--making a direct link to the correlated wavefunction theory. The model has been trained for 15 elements (H, Li-F, Na-Cl, 720 parameters) on a set of 5581 molecules (including ions, transition states, and weakly bound complexes) whose first- and second-order properties were computed by the coupled-cluster theory as a reference, and a good agreement is seen. The model looks promising for the study of large molecular systems, it is believed to be an important step forward from the traditional semiempirical models towards higher accuracy at nearly as low a computational cost.

Communication: stabilization of radical anions with weakly bound electron in condensed media: a case study of diacetonyl radical anion.

The radical anion resulting from electron capture by diacetonyl molecule has been characterized by EPR and optical absorption spectroscopy in glassy ether matrices at 77 K. In non-polar alkane glasses this species was not observed under the same conditions, which confirms the crucial role of matrix interactions in stabilizing this species. Calculations at the MP2 level show the vertical detachment energy to increase gradually from roughly zero for a bare anion to ∼1 eV for the complex involving six ether molecules.

Neglect of four- and approximation of one-, two-, and three-center two-electron integrals in a symmetrically orthogonalized basis.

A new integral approximation for use in molecular electronic structure calculations is proposed as an alternative to the traditional neglect of diatomic differential overlap models. The similarity between the symmetrically orthogonalized and the original basis functions (assumed orthonormal within each atomic set but nonorthogonal between different centers) is used to construct a robust approximation for the two-electron integrals, with the error being quadratic in the deviation between the products of the functions. Invariance properties of this procedure are rigorously proved. Numerical studies on a representative set of molecules at valence-only minimal basis Hartree-Fock level show that the approximation introduces relatively small errors, encouraging its future application in the semiempirical field.

High-resolution electron spin resonance spectroscopy of XeF* in solid argon. The hyperfine structure constants as a probe of relativistic effects in the chemical bonding properties of a heavy noble gas atom.

Xenon fluoride radicals were generated by solid-state chemical reactions of mobile fluorine atoms with xenon atoms trapped in Ar matrix. Highly resolved electron spin resonance spectra of XeF* were obtained in the temperature range of 5-25 K and the anisotropic hyperfine parameters were determined for magnetic nuclei 19F, 129Xe, and 131Xe using naturally occurring and isotopically enriched xenon. Signs of parallel and perpendicular hyperfine components were established from analysis of temperature changes in the spectra and from numerical solutions of the spin Hamiltonian for two nonequivalent magnetic nuclei. Thus, the complete set of components of hyperfine- and g-factor tensors of XeF* were obtained: 19F (Aiso=435, Adip=1249 MHz) and 129Xe (Aiso=-1340, Adip=-485 MHz); g(parallel)=1.9822 and g(perpendicular)=2.0570. Comparison of the measured hyperfine parameters with those predicted by density-functional theory (DFT) calculations indicates, that relativistic DFT gives true electron spin distribution in the 2Sigma+ ground-state, whereas nonrelativistic theory underestimates dramatically the electron-nuclear contact Fermi interaction (Aiso) on the Xe atom. Analysis of the obtained magnetic-dipole interaction constants (Adip) shows that fluorine 2p and xenon 5p atomic orbitals make a major contribution to the spin density distribution in XeF*. Both relativistic and nonrelativistic calculations give close magnetic-dipole interaction constants, which are in agreement with the measured values. The other relativistic feature is considerable anisotropy of g-tensor, which results from spin-orbit interaction. The orbital contribution appears due to mixing of the ionic 2Pi states with the 2Sigma+ ground state, and the spin-orbit interaction plays a significant role in the chemical bonding of XeF*.