No Arabic abstract
A reliable and efficient computation of the entire single-particle spectrum of correlated molecules is an outstanding challenge in the field of quantum chemistry, with standard density functional theory approaches often giving an inadequate description of excitation energies and gaps. In this work, we expand upon a recently-introduced approach which relies on a fully self-consistent many-body perturbation theory, coupled to a non-perturbative truncation of the effective dynamics at each step. We show that this yields a low-scaling and accurate method across a diverse benchmark test set, capable of treating moderate levels of strong correlation effects, and detail an efficient implementation for applications up to $sim1000$ orbitals on parallel resources. We then use this method to characterise the spectral properties of the artemisinin anti-malarial drug molecule, resolving discrepancies in previous works concerning the active sites of the lowest energy fundamental excitations of the system.
The possibility of using time-resolved vibronic spectroscopy for spectral analysis of mixtures of chemical compounds with similar optical properties, when traditional methods are inefficient, is demonstrated by using the method of computer simulation. The analysis is carried out by the example of molecules of a series of polyenes (butadiene, hexatraene, octatetraene, decapentaene, and decatetraene), their various cis- and trans-rotational isomers, and phenyl-substituted polyenes. Ranges of relative concentrations of molecules similar in their spectral properties, where reliable interpretation of time-resolved spectra of mixtures and both qualitative and quantitative analyses are possible, are determined. The use of computer simulation methods for oprimizing full-scale experiments in femtosecond spectroscopy is shown to hold much promise.
Density matrix perturbation theory (DMPT) is known as a promising alternative to the Rayleigh-Schrodinger perturbation theory, in which the sum-over-state (SOS) is replaced by algorithms with perturbed density matrices as the input variables. In this article, we formulate and discuss three types of DMPT, with two of them based only on density matrices: the approach of Kussmann and Ochsenfeld [J. Chem. Phys.127, 054103 (2007)] is reformulated via the Sylvester equation, and the recursive DMPT of A.M.N. Niklasson and M. Challacombe [Phys. Rev. Lett. 92, 193001 (2004)] is extended to the hole-particle canonical purification (HPCP) from [J. Chem. Phys. 144, 091102 (2016)]. Comparison of the computational performances shows that the aformentioned methods outperform the standard SOS. The HPCP-DMPT demonstrates stable convergence profiles but at a higher computational cost when compared to the original recursive polynomial method
Approximate natural orbitals are investigated as a way to improve a Monte Carlo configuration interaction (MCCI) calculation. We introduce a way to approximate the natural orbitals in MCCI and test these and approximate natural orbitals from MP2 and QCISD in MCCI calculations of single-point energies. The efficiency and accuracy of approximate natural orbitals in MCCI potential curve calculations for the double hydrogen dissociation of water, the dissociation of carbon monoxide and the dissociation of the nitrogen molecule are then considered in comparison with standard MCCI when using full configuration interaction as a benchmark. We also use the method to produce a potential curve for water in an aug-cc-pVTZ basis. A new way to quantify the accuracy of a potential curve is put forward that takes into account all of the points and that the curve can be shifted by a constant. We adapt a second-order perturbation scheme to work with MCCI (MCCIPT2) and improve the efficiency of the removal of duplicate states in the method. MCCIPT2 is tested in the calculation of a potential curve for the dissociation of nitrogen using both Slater determinants and configuration state functions.
We identify the dominant computational cost within the recently introduced stochastic and internally contracted FCIQMC-NEVPT2 method for large active space sizes. This arises from the contribution to the four-body intermediates arising from low-excitation level sampled determinant pairs. We develop an effective way to mitigate this cost via an additional stochastic step within the sampling of the required NEVPT2 intermediates. We find this systematically improvable additional sampling can reduce simulation time by 80% without introducing appreciable error. This saving is expected to increase for larger active spaces. We combine this enhanced sampling scheme with full stochastic orbital optimization for the first time, and apply it to find FCIQMC-NEVPT2 energies for spin states of an iron porphyrin system within (24,24) active spaces with relatively meagre computational resources. This active space size can now be considered as routine for NEVPT2 calculations of strongly correlated molecular systems within this improved stochastic methodology.
A semiempirical parametric method is proposed for modeling three-dimensional (time-resolved) vibronic spectra of polyatomic molecules. The method is based on the use of the fragment approach in the formation of molecular models for excited electronic states and parametrization of these molecular fragments by modeling conventional (one-dimensional) absorption and fluorescence spectra of polyatomic molecules. All matrix elements that are required for calculation of the spectra can be found by the methods developed. The time dependencies of the populations of a great number (>10^3) of vibronic levels can be most conveniently found by using the iterative numerical method of integration of kinetic equations. Convenient numerical algorithms and specialized software for PC are developed. Computer experiments showed the possibility of the real-time modeling of three-dimensional spectra of polyatomic molecules containing several tens of atoms.