Ions have a profound effect on the geometrical structure of liquid water and an aqueous environment is known to change the electronic structure of ions.

We explore the range of applicability of the nuclear ensemble method (NEM) for quantitative simulations of absorption spectra and their temperature variations. We formulate a "good practice" for the NEM based on statistical theory. Special attention is paid to proper treatment of uncertainty estimation including the convergence with the number of samples, which is often neglected in the field. As a testbed, we have selected a well-known chromophore, (E)-azobenzene. We measured its temperature difference UV− vis absorption spectra in methanol, which displayed two dominant features: a moderate increase in the intensity of the nπ* band and a pronounced decrease in intensity of the low-energy part of the ππ* band. We attributed both features to increasing non-Condon effects with temperature. We show that the NEM based on the path integral molecular dynamics combined with range-separated hybrid functionals provides quantitatively accurate spectra and their differences. Experimentally, the depletion of the absorption in the ππ* band showed a characteristic vibrational progression that cannot be reproduced with the NEM. We show that hundreds of thousands of samples are necessary to achieve an accuracy sufficient for the unambiguous explanation of the observed temperature effects. We provide a detailed analysis of the temperature effects on the spectrum based on the harmonic model of the system combined with the NEM. We also rationalize the vibrational structure of the spectrum using the Franck−Condon principle.

Criegee intermediates (CIs) are increasingly recognized as important intermediates in atmospheric and combustion chemistry. Modelling in these fields requires reliable characterization of the CI's UV absorption cross sections. Different experimental approaches provided seemingly contradictory results both for the simplest CI as well as for higher analogues. To resolve these discrepancies, we model simultaneously the positions, widths, shapes and absolute intensities of the first absorption bands of CIs using the semiclassical reflection principle approach based on the path integral molecular dynamics. The UV spectra were efficiently processed via the kernel density estimation method. We show that the multireference character complicates the description of the system and the appropriate treatment of the electron correlation is vital. Multi-reference methods with dynamical correlation are problematic due to high sensitivity of the results with respect to the active space. Single-reference based methods become reliable once the triple excitations components are included. Interestingly, the very similar CC2 and ADC(2) methods provide highly conflicting results which are, however, reconciled at the CC3 and ADC(3) levels. To calculate the spectra of higher CIs, we introduce composite EOM-CCSD/T and CC2/3 methods providing excellent absorption cross sections at an acceptable computational cost. We provide the first absolute intensities for individual conformers of the CHCHOO molecule.