Rigid urea and self-healing thiourea ethanolamine monolayers

A series of long-tail alkyl ethanolamine analogs containing amide-, urea-, and thiourea moieties was synthesized and the behavior of the corresponding monolayers was assessed on the Langmuir–Pockels trough combined with grazing incidence X-ray diffraction experiments and complemented by computer simulations. All compounds form stable monolayers at the soft air/water interface. The phase behavior is dominated by strong intermolecular headgroup hydrogen bond networks. While the amide analog forms well-defined monolayer structures, the stronger hydrogen bonds in the urea analogs lead to the formation of small three-dimensional crystallites already during spreading due to concentration fluctuations. The hydrogen bonds in the thiourea case form a two-dimensional network, which ruptures temporarily during compression and is recovered in a self-healing process, while in the urea clusters the hydrogen bonds form a more planar framework with gliding planes keeping the structure intact during compression. Because the thiourea analogs are able to self-heal after rupture, such compounds could have interesting properties as tight, ordered, and self-healing monolayers

Software for the frontiers of quantum chemistry:An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022

Embedding potentials for excited states of embedded species

Frozen-Density-Embedding Theory (FDET) is a formalism to obtain the upper bound of the ground-state energy of the total system and the corresponding embedded wavefunction by means of Euler-Lagrange equations [T. A. Wesolowski, Phys. Rev. A77(1), 012504 (2008)]. FDET provides the expression for the embedding potential as a functional of the electron density of the embedded species, electron density of the environment, and the field generated by other charges in the environment. Under certain conditions, FDET leads to the exact ground-state energy and density of the whole system. Following Perdew-Levy theorem on stationary states of the ground-state energy functional, the other-than-ground-state stationary states of the FDET energy functional correspond to excited states. In the present work, we analyze such use of other-than-ground-state embedded wavefunctions obtained in practical calculations, i.e., when the FDET embedding potential is approximated. Three computational approaches based on FDET, that assure self-consistent excitation energy and embedded wavefunction dealing with the issue of orthogonality of embedded wavefunctions for different states in a different manner, are proposed and discussed

Development of Novel Computational Strategies to Match the Challenges of Supramolecular Chemistry, Biochemistry, and Materials Science

Recent formal developments and applications of the 'freeze-and-conquer' strategy proposed by Wesolowski and Warshel in 1993 to study large systems at quantum mechanical level are reviewed. This universal approach based on density functional theory allows one to link, via the orbital-free embedding potential, two parts of a larger system described at different levels of accuracy leading thus to significant savings in computational costs. As a result, applicability of conventional methods of quantum chemistry can be extended to even larger systems. It is shown that the 'freeze-and-thaw' approach applying the first-principles based approximation to the orbital-free embedding potential recently developed in our group provides a powerful and universal technique to study such embedded molecules (or molecular complexes), which are not linked with their microscopic environment by covalent bonds

Effect of the f-Orbital Delocalization on the Ligand-Field Splitting Energies in Lanthanide-Containing Elpasolites

The ligand-field induced splitting energies of f-levels in lanthanide-containing elpasolites are derived using the first-principles universal orbital-free embedding formalism [Wesolowski and Warshel, J. Phys. Chem. 1993, 97, 8050]. In our previous work concerning chloroelpasolite lattice (Cs2NaLnCl6), embedded orbitals and their energies were obtained using an additional assumption concerning the localization of embedded orbitals on preselected atoms leading to rather good ligand-field parameters. In this work, the validity of the localization assumption is examined by lifting it. In variational calculations, each component of the total electron density (this of the cation and that of the ligands) spreads over the whole system. It is found that the corresponding electron densities remain localized around the cation and the ligands, respectively. The calculated splitting energies of f-orbitals in chloroelpasolites are not affected noticeably in the whole lanthanide series. The same computational procedure is used also for other elpasolite lattices (Cs2NaLnX6, where X=F, Br, and I)materials which have not been fabricated or for which the ligand-field splitting parameters are not available

Ground States, Excited States and Metal-Ligand Bonding in Rare Earth Hexachloro Complexes: A DFT Based Ligand Field Study

Metal (4f)-ligand (Cl 3p) bonding in LnCl63- (Ln = Ce to Yb) complexes has been studied on the basis of 4f->4f and Cl,3p->4f charge-transfer spectra and on the analysis of these spectra within the valence bond configuration interaction model to show that mixing of Cl 3p into the Ln 4f ligand field orbitals does not exceed 1%. Contrary to this, Kohn-Sham formalism of density functional theory using currently available approximations to the exchange-correlation functional tends to strongly overestimate 4f-3p covalency, yielding, for YbCl63-, a much larger mixing of Cl 3p->4f charge transfer into the f13 ionic ground-state wave function. Thus, ligand field density functional theory, which was recently developed and applied with success to complexes of 3d metals in our group, yields anomalously large ligand field splittings for Ln, the discrepancy with experiment increasing from left to the right of the Ln 4f series. It is shown that eliminating artificial ligand-to-metal charge transfer in Kohn-Sham calculations by a procedure described in this work leads to energies of 4f-4f transitions in good agreement with experiment. We recall an earlier concept of Ballhausen and Dahl which describes ligand field in terms of a pseudopotential and give a thorough analysis of the contributions to the ligand field from electrostatics (crystal field) and exchange (Pauli) repulsion. The close relation of the present results with those obtained using the first-principles based and electron density dependent effective embedding potential is pointed out along with implications for applications to other systems

Extracting Information about Chemical Bonding from Molecular Electron Densities via Single Exponential Decay Detector (SEDD)

The recently introduced molecular descriptor (Single Exponential Decay Detector - SEDD) [P. de Silva, J. Korchowiec, T. A. Wesolowski, ChemPhysChem 2012, 13, 3462] is used to visualize bonding patterns in molecules. In each point of space SEDD is simply related to the electron density: SEDD(r) = ln[1/?2(?(??/?)2)2]. Either experimental or computed densities ?(r) can be used to evaluate SEDD. Here, maps of SEDD are obtained from theoretical densities and reveal such features as core electrons, chemical bonds, lone pairs and delocalization in aromatic systems. It is shown that SEDD provides fingerprints of aromaticity, which can be separated into geometric and electronic effects