Benchmark theoretical study of the electric polarizabilities of naphthalene, anthracene, and tetracene

Accurate prediction of global-density-dependent range-separation parameters based on machine learning

In this work, we develop an accurate and efficient XGBoost machine learning model for predicting the global-density-dependent range-separation parameter, ωGDD, for long-range corrected functional (LRC)-ωPBE. This ωGDDML model has been built using a wide range of systems (11 466 complexes, ten different elements, and up to 139 heavy atoms) with fingerprints for the local atomic environment and histograms of distances for the long-range atomic correlation for mapping the quantum mechanical range-separation values. The promising performance on the testing set with 7046 complexes shows a mean absolute error of 0.001 117 a0-1 and only five systems (0.07%) with an absolute error larger than 0.01 a0-1, which indicates the good transferability of our ωGDDML model. In addition, the only required input to obtain ωGDDML is the Cartesian coordinates without electronic structure calculations, thereby enabling rapid predictions. LRC-ωPBE(ωGDDML) is used to predict polarizabilities for a series of oligomers, where polarizabilities are sensitive to the asymptotic density decay and are crucial in a variety of applications, including the calculations of dispersion corrections and refractive index, and surpasses the performance of all other popular density functionals except for the non-tuned LRC-ωPBE. Finally, LRC-ωPBE (ωGDDML) combined with (extended) symmetry-adapted perturbation theory is used in calculating noncovalent interactions to further show that the traditional ab initio system-specific tuning procedure can be bypassed. The present study not only provides an accurate and efficient way to determine the range-separation parameter for LRC-ωPBE but also shows the synergistic benefits of fusing the power of physically inspired density functional LRC-ωPBE and the data-driven ωGDDML model.

Solvation Effects on Polarizability of Aromatic Fluids

Elucidating solvation effects on polarizability in condensed phases is important for the description of the optical and dielectric behavior of high-refractive-index molecular materials. We study these effects via the polarizability model combining electronic, solvation, and vibrational contributions. The method is applied to well-characterized highly polarizable liquid precursors: benzene, naphthalene, and phenanthrene. We find that the solvation and vibrational terms are of opposite signs and cancel almost exactly for benzene, but for naphthalene and phenanthrene, a 2.5 and 5.0% decrease relative to the equilibrium electronic polarizability of the respective monomer, α₁^e, is predicted, respectively. The increase in electronic polarizability affects interaction polarizability of all contacts, which is the main reason for the increasing importance of solvation contribution. The calculated refractive indices agree very well with experiment for all three systems.

An Imbalance in the Force: The Need for Standardized Benchmarks for Molecular Simulation

Force fields (FFs) for molecular simulation have been under development for more than half a century. As with any predictive model, rigorous testing and comparisons of models critically depends on the availability of standardized data sets and benchmarks. While such benchmarks are rather common in the fields of quantum chemistry, this is not the case for empirical FFs. That is, few benchmarks are reused to evaluate FFs, and development teams rather use their own training and test sets. Here we present an overview of currently available tests and benchmarks for computational chemistry, focusing on organic compounds, including halogens and common ions, as FFs for these are the most common ones. We argue that many of the benchmark data sets from quantum chemistry can in fact be reused for evaluating FFs, but new gas phase data is still needed for compounds containing phosphorus and sulfur in different valence states. In addition, more nonequilibrium interaction energies and forces, as well as molecular properties such as electrostatic potentials around compounds, would be beneficial. For the condensed phases there is a large body of experimental data available, and tools to utilize these data in an automated fashion are under development. If FF developers, as well as researchers in artificial intelligence, would adopt a number of these data sets, it would become easier to compare the relative strengths and weaknesses of different models and to, eventually, restore the balance in the force.

Investigation of reactive properties, adsorption on fullerene, DFT, molecular dynamics simulation of an anthracene derivative targeting dihydrofolate reductase and human dUTPase

Anthracenes are aromatic compounds with flexible structure and reactivity which are of great interest to theoretical and experimental chemists. Theoretical investigations of 1,4-dihydroxy-5,8-bis[2-(2-hydroxyethylamino)ethylamino]anthracene-9,10-dione (Mitoxantrone) (DDEA) based on density functional theory, molecular dynamics and adsorption on fullerene are reported in the present research. The suitable situation for adsorption with fullerene (C60) is the cyclohex-2-ene-1,4-dione ring of DDEA. Selected quantum-molecular descriptors have been calculated to predict the most reactive sites of the DDEA molecule. Interactions of DDEA with water have been studied using MD simulations. MD simulations were also used to study solubility parameter, a significant quantity for the development of pharmaceutical formulations. The affinity of DDEA on human dihydrofolate reductase and deoxyuridine triphosphatase enzymes was investigated by MD simulation of the protein-ligand complex obtained by molecular docking study.Communicated by Ramaswamy H. Sarma.

Evaluating the Solvent Stark Effect from Temperature‐Dependent Solvatochromic Shifts of Anthracene

The solvent Stark effect on the spectral shifts of anthracene is studied with temperature‐dependent solvatochromic measurements. The Stark contribution Δv_Stark to the absorption shift Δv_p in polar solvents is measured to be Δv_Stark=(53±35) cm⁻¹, in reasonable agreement with dielectric continuum theory estimate of 28 cm⁻¹, whereas the major shift Δv_p∼300 cm⁻¹ presumably originates from the solute quadrupole. We pay attention to the accurate correction of Δv_p for the nonpolar contribution that is crucial when the shifts are modest in magnitude.

Accurate molecular polarizabilities with coupled cluster theory and machine learning

The molecular dipole polarizability describes the tendency of a molecule to change its dipole moment in response to an applied electric field. This quantity governs key intra- and intermolecular interactions, such as induction and dispersion; plays a vital role in determining the spectroscopic signatures of molecules; and is an essential ingredient in polarizable force fields. Compared with other ground-state properties, an accurate prediction of the molecular polarizability is considerably more difficult, as this response quantity is quite sensitive to the underlying electronic structure description. In this work, we present highly accurate quantum mechanical calculations of the static dipole polarizability tensors of 7,211 small organic molecules computed using linear response coupled cluster singles and doubles theory (LR-CCSD). Using a symmetry-adapted machine-learning approach, we demonstrate that it is possible to predict the LR-CCSD molecular polarizabilities of these small molecules with an error that is an order of magnitude smaller than that of hybrid density functional theory (DFT) at a negligible computational cost. The resultant model is robust and transferable, yielding molecular polarizabilities for a diverse set of 52 larger molecules (including challenging conjugated systems, carbohydrates, small drugs, amino acids, nucleobases, and hydrocarbon isomers) at an accuracy that exceeds that of hybrid DFT. The atom-centered decomposition implicit in our machine-learning approach offers some insight into the shortcomings of DFT in the prediction of this fundamental quantity of interest.

Polarizable Drude Model with s-Type Gaussian or Slater Charge Density for General Molecular Mechanics Force Fields

Gas-phase electric properties of molecules can be computed routinely using wave function methods or density functional theory (DFT). However, these methods remain computationally expensive for high-throughput screening of the vast chemical space of virtual compounds. Therefore, empirical force fields are a more practical choice in many cases, particularly since force field methods allow one to routinely predict the physicochemical properties in the condensed phases. This work presents Drude polarizable models, to increase the physical realism in empirical force fields, where the core particle is treated as a point charge and the Drude particle is treated either as a 1 s-Gaussian or a ns-Slater ( n = 1, 2, 3) charge density. Systematic parametrization to large high-quality quantum chemistry data obtained from the open access Alexandria Library ( https://doi.org/10.5281/zenodo.1004711 ) ensures the transferability of these parameters. The dipole moments and isotropic polarizabilities of the isolated molecules predicted by the proposed Drude models are in agreement with experiment with accuracy similar to DFT calculations at the B3LYP/aug-cc-pVTZ level of theory. The results show that the inclusion of explicit polarization into the models reduces the root-mean-square deviation with respect to DFT calculations of the predicted dipole moments of 152 dimers and clusters by more than 50%. Finally, we show that the accuracy of the electrostatic interaction energy of the water dimers can be improved systematically by the introduction of polarizable smeared charges as a model for charge penetration.

The impact of chemical structure and molecular packing on the electronic polarisation of fullerene arrays

Electronic polarisation contributes to the electronic landscape as seen by separating charges in organic materials. The nature of electronic polarisation depends on the polarisability, density, and arrangement of polarisable molecules. In this paper, we introduce a microscopic, coarse-grained model in which we treat each molecule as a polarisable site, and use an array of such polarisable dipoles to calculate the electric field and associated energy of any arrangement of charges in the medium. The model incorporates chemical structure via the molecular polarisability and molecular packing patterns via the structure of the array. We use this model to calculate energies of charge pairs undergoing separation in finite fullerene lattices of different chemical and crystal structures. The effective dielectric constants that we estimate from this approach are in good quantitative agreement with those measured experimentally in C₆₀ and phenyl-C₆₁-butyric acid methyl ester (PCBM) films, but we find significant differences in dielectric constant depending on packing and on direction of separation, which we rationalise in terms of density of polarisable fullerene cages in regions of high field. In general, we find lattices containing molecules of more isotropic polarisability tensors exhibit higher dielectric constants. By exploring several model systems we conclude that differences in molecular polarisability (and therefore, chemical structure) appear to be less important than differences in molecular packing and separation direction in determining the energetic landscape for charge separation. We note that the results are relevant for finite lattices, but not necessarily for infinite systems. We propose that the model could be used to design molecular systems for effective electronic screening.

Theoretical study of molecular vibrations in electron momentum spectroscopy experiments on furan: an analytical versus a molecular dynamical approach

The influence of thermally induced nuclear dynamics (molecular vibrations) in the initial electronic ground state on the valence orbital momentum profiles of furan has been theoretically investigated using two different approaches. The first of these approaches employs the principles of Born-Oppenheimer molecular dynamics, whereas the so-called harmonic analytical quantum mechanical approach resorts to an analytical decomposition of contributions arising from quantized harmonic vibrational eigenstates. In spite of their intrinsic differences, the two approaches enable consistent insights into the electron momentum distributions inferred from new measurements employing electron momentum spectroscopy and an electron impact energy of 1.2 keV. Both approaches point out in particular an appreciable influence of a few specific molecular vibrations of A1 symmetry on the 9a1 momentum profile, which can be unravelled from considerations on the symmetry characteristics of orbitals and their energy spacing.

Models of charge pair generation in organic solar cells

A critical perspective on modelling of charge generation in organic photovoltaics, focussing on interfacial electronic states, electrostatics, and dynamic processes.