Hilbert space multireference coupled cluster tailored by matrix product states.

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants-the state universal one and the Brillouin-Wigner and Mukherjee's state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.

Three-body potential and third virial coefficients for helium including relativistic and nuclear-motion effects.

The non-additive three-body interaction potential for helium was computed using the coupled-cluster theory and the full configuration interaction method. The obtained potential comprises an improved nonrelativistic Born-Oppenheimer energy and the leading relativistic and nuclear-motion corrections. The mean absolute uncertainty of our calculations due to the incompleteness of the orbital basis set was determined employing complete-basis-set extrapolation techniques and was found to be 1.2%. For three helium atoms forming an equilateral triangle with the side length of 5.6 bohr - a geometry close to the minimum of the total potential energy surface - our three-body potential amounts to -90.6 mK, with an estimated uncertainty of 0.5 mK. An analytic function, developed to accurately fit the computed three-body interaction energies, was chosen to correctly describe the asymptotic behavior of the three-body potential for trimer configurations corresponding to both the three-atomic and the atom-diatom fragmentation channels. For large triangles with sides r12, r23, and r31, the potential takes correctly into account all angular terms decaying as r-l12 r-m23 r-n21 with l + m + n ≤ 14 for the nonrelativistic Born-Oppenheimer energy and l + m + n ≤ 9 for the post-Born-Oppenheimer corrections. We also developed a short-range analytic function describing the local behavior of the total uncertainty of the computed three-body interaction energies. Using both fits we calculated the third pressure and acoustic virial coefficients for helium and their uncertainties for a wide range of temperatures. The results of these calculations were compared with available experimental data and with previous theoretical determinations. The estimated uncertainties of present calculations are 3-5 times smaller than those reported in the best previous works.

Comparing the effects of different exogenous hormone combinations on seed-derived callus induction in Nicotiana tabacum.

Callus from Nicotiana tabacum is used as a model in plant developmental research. We tested several phytohormone (Indoleacetic acid - IAA; 2,4-Dichlorophenoxyacetic acid - 2,4-D; kinetin - KIN; 6-Benzylaminopurine - BAP) combinations to compare different approaches to callus induction directly from the seeds of Nicotiana tabacum. Callus formation was observed up to 4 weeks after sowing and the most effective were 0.5 mg/L of 2,4-D with 0.25 mg/L of BAP and 2 mg/L 2,4-D with 1 mg/L of BAP. The calli were green, photosynthetically active and after 6 weeks of growth, no stress symptoms (estimated on the basis of fluorescence of chlorophyll a in photosystem II) were noticed.

Ground state of the Fe(ii)-porphyrin model system corresponds to quintet: a DFT and DMRG-based tailored CC study.

Fe(ii)-porphyrins play an important role in many reactions relevant to material science and biological processes, due to their closely lying spin states. Although the prevalent opinion is that these systems posses the triplet ground state, the recent experiment on Fe(ii)-phthalocyanine under conditions matching those of an isolated molecule points toward the quintet ground state. We present a thorough DFT and DMRG-based tailored CC study of Fe(ii)-porphyrin model, in which we address all previously discussed correlation effects. We examine the importance of geometrical parameters, the Fe-N distances in particular, and conclude that the system possesses the quintet ground state.

Near-Linear Scaling in DMRG-Based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T).

We present a new implementation of density matrix renormalization group based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO). Compared to the previous local pair natural orbital (LPNO) version of the method, the new implementation is more accurate, offers more favorable scaling, and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen), and iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8 and 99.9% of the total canonical correlation energy which is an improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.

A water-soluble [60]fullerene-derivative stimulates chlorophyll accumulation and has no toxic effect on Chlamydomonas reinhardtii.

Chlamydomonas reinhardtii (WT 2137) P. A. Dang. (Volvocales, Chlorophyceae) is a green microalgae serving as a suitable model in scientific research and a promising industrial biotechnology platform for production of biofuel, hydrogen and recombinant proteins. Fullerenes (C60) are allotropic carbon nanoparticles discovered in 1985 and used in biomedical studies since the early 1990s, when water solubilization methodologies were developed. Recently, surface-modified hydroxylated derivatives of fullerenes were proven to enhance algal growth and drought tolerance in plants. Here, a novel type of water-soluble [60]fullerene derivative with 12 glycine residues (GF) has been synthesized and tested for acute toxicity (up to 50 µg/ml) and as a potential biostimulant of algal growth. The effects of GF on pigment composition and growth rate of Chlamydomonas reinhardtii were systematically investigated. Our results suggest that GF was not toxic, and no negative change in the pigment content and no stress symptoms were observed. No changes in the photosynthetic parameters based on the fluorescence of chlorophyll a in Photosystem II (NPQ, Fv/Fm, Fv/F0, PI and RC/ABS) were observed. The GF had no effect on cell size and growth rate. At a concentration of 20 µg/ml, GF stimulated chlorophyll accumulation in 3-day-old cultures.

Perturbative triples correction to domain-based local pair natural orbital variants of Mukherjee's state specific coupled cluster method.

In this article we report an implementation of the perturbative triples correction to Mukherjee's state-specific multireference coupled cluster method based on the domain-based pair natural orbital approach (DLPNO-MkCC). We tested the performance of DLPNO-MkCCSD(T) in calculations involving tetramethyleneethane and isomers of naphthynes. These tests show that more than 97% of triples energy was recovered with respect to the canonical MkCCSD(T) method, which together with the DLPNO-MkCCSD part accounts for about 99.70-99.85% of the total correlation energy. The applicability of the method was demonstrated on calculations of singlet-triplet gaps for several large systems: triangulene, dynemicin A, and a beryllium complex.

Domain-Based Local Pair Natural Orbital Version of Mukherjee's State-Specific Coupled Cluster Method.

This article reports development of a local variant of Mukherjee's state-specific multireference coupled cluster method based on the domain-based pair natural orbital approach (DLPNO-MkCC). The current implementation is restricted to connected single and double excitations and model space with up to biexcited references. The performance of the DLPNO-MkCCSD was tested on calculations of tetramethyleneethane. The results show that above 99.9% of the correlation energy was recovered, with respect to the conventional MkCC method. To demonstrate the applicability of the method to large systems, singlet-triplet gaps of triangulene and bis(1-(2,6-diisopropylphenyl)-3,3,5,5-tetramethylpyrrolidine-2-ylidene)beryllium complex were studied. For the last system (105 atoms), we were able to perform a calculation in cc-pVTZ with 2158 basis functions on a single CPU in less than 9 days.