Search

We investigated some variational methods to compute a wavefunction based on antisymmetric product of geminals (APG). The Waring decomposition on the APG wavefunction leads a finite sum of antisymmetrized geminal power (AGP) wavefunctions, each for which the variational principle can be applied. We call this as AGP-CI method which provides a variational solution of the APG wavefunction efficiently. However, number of AGP wavefunctions in the exact AGP-CI formalism become exponentially large in case of many-electron systems. Therefore, we also investigate the low-rank APG wavefunction, in which the geminal matrices are factorized by the Schur decomposition. Interestingly, only a few non-zero eigenvalues (up to half number of electrons) were found from the Schur decomposition on the APG wavefunction. We developed some methods to approximate the APG wavefunction by lowering the ranks of geminal matrices, and demonstrate their performance. Our new geminal method based on the low-rank decomposition can drastically reduce the number of variational parameters, although there is no efficient algorithm so far, due to some mathematical complications., Comment: 21 pages, 7 figures

In the absence of laboratory data, state-of-the-art quantum chemical approaches can provide estimates of the binding energy (BE) of interstellar species with grains. Without BE values, contemporary astrochemical models are compelled to utilize wild guesses, often delivering misleading information. Here, we employed a fully quantum chemical approach to estimate the BE of seven diatomic radicals - CH, NH, OH, SH, CN, NS, and NO - that play a crucial role in shaping the interstellar chemical composition, using a suitable amorphous solid water model as a substrate since water is the principal constituent of interstellar ice in dense and shielded regions. While the BEs are compatible with physisorption, the binding of CH in some sites shows chemisorption, in which a chemical bond to an oxygen atom of a water molecule is formed. While no structural change has been observed for the CN radical, it is believed that the formation of a hemibonded system between the outer layer of the water cluster and the radical is the reason for the unusually large BE in one of the binding sites considered in our study. A significantly lower BE for NO, consistent with recent calculations, is obtained, which helps explain the recently observed HONO/NH$_2$OH and HONO/HNO ratios in the low-mass hot corino IRAS 16293-2422 B with chemical models., Comment: 11 pages, 6 figures, to be published in Astronomy & Astrophysics (A&A) Journal

Currently, there is a growing interest in the development of a new hierarchy of methods based on the concept of seniority, which has been introduced quite recently in quantum chemistry. Despite the enormous potential of these methods, the accurate description of both dynamical and static correlation effects within a single and in-principle-exact approach remains a challenge. In this work, we propose an alternative formulation of reduced density-matrix functional theory (RDMFT) where the (one-electron reduced) density matrix is mapped onto an ab initio seniority-zero wave function. In this theory, the exact natural orbitals and their occupancies are determined self-consistently from an effective seniority-zero calculation. The latter involves a universal higher-seniority density matrix functional for which an adiabatic connection (AC) formula is derived and implemented under specific constraints that are related to the density matrix. The pronounced curvature of the (constrained) AC integrand, which is numerically observed in prototypical hydrogen chains and the Helium dimer, indicates that a description of higher-seniority correlations within second-order perturbation theory is inadequate in this context. Applying multiple linear interpolations along the AC or connecting second-order perturbation theory to a full-seniority treatment via Pad\'{e} approximants are better strategies. Such information is expected to serve as a guide in the future design of higher-seniority density-matrix functional approximations.

It has been implicitly assumed that ices on grains in molecular clouds and proto planetary disks are formed by homogeneous layers regardless of their composition or crystallinity. To verify this assumption, we observed the H2O deposition onto refractory substrates and the crystallization of amorphous ices (H2O, CO2, and CO) using an ultra-high-vacuum transmission electron microscope. In the H2O-deposition experiments, we found that three-dimensional islands of crystalline ice (Ic) were formed at temperatures above 130 K. The crystallization experiments showed that uniform thin films of amorphous CO and H2O became three-dimensional islands of polyhedral crystals; amorphous CO2, on the other hand, became a thin film of nano crystalline CO2 covering the amorphous H2O. Our observations show that crystal morphologies strongly depend not only on the ice composition, but also on the substrate. Using experimental data concerning the crystallinity of deposited ices and the crystallization timescale of amorphous ices, we illustrated the criteria for ice crystallinity in space and outlined the macroscopic morphology of icy grains in molecular clouds as follows: amorphous H2O covered the refractory grain uniformly, CO2 nano-crystals were embedded in the amorphous H2O, and a polyhedral CO crystal was attached to the amorphous H2O. Furthermore, a change in the grain morphology in a proto-planetary disk is shown. These results have important implications for the chemical evolution of molecules, non-thermal desorption, collision of icy grains, and sintering., Comment: 60 pages, 18 figures, 3 tables, published in ApJ

It is speculated that there might be some linkage between interstellar aldehydes and their corresponding alcohols. Here, an observational study and astrochemical modeling are coupled together to illustrate the connection between them. The ALMA Cycle 4 data of a hot molecular core, G10.47+0.03 is utilized for this study. Various aldehydes (acetaldehyde, propanal, and glycolaldehyde), alcohols (methanol and ethylene glycol), and a ketone (acetone) are identified in this source. The excitation temperatures and the column densities of these species were derived via the rotation diagram method assuming LTE conditions. An extensive investigation is carried out to understand the formation of these species. Six pairs of aldehyde-alcohol: i) methanal and methanol; ii) ethanal and ethanol; iii) propanal and 1-propanol; iv) propenal and allyl alcohol; v) propynal and propargyl alcohol; vi) glycolaldehyde and ethylene glycol; vii) along with one pair of ketone-alcohol (acetone and isopropanol) and viii) ketene-alcohol (ethenone and vinyl alcohol) are considered for this study. Two successive hydrogenation reactions in the ice phase are examined to form these alcohols from aldehydes, ketone, and ketene, respectively. Quantum chemical methods are extensively executed to review the ice phase formation route and the kinetics of these species. Based on the obtained kinetic data, astrochemical modeling is employed to derive the abundances of these aldehydes, alcohols, ketone, and ketene in this source. It is seen that our model could successfully explain the observed abundances of various species in this hot molecular core., Comment: 34 pages, 11 figures, Accepted for the publication in the Astrophysical Journal

Phosphorus related species are not known to be as omnipresent in space as hydrogen, carbon, nitrogen, oxygen, and sulfur-bearing species. Astronomers spotted very few P-bearing molecules in the interstellar medium and circumstellar envelopes. Limited discovery of the P-bearing species imposes severe constraints in modeling the P-chemistry. In this paper, we carry out extensive chemical models to follow the fate of P-bearing species in diffuse clouds, photon-dominated or photodissociation regions (PDRs), and hot cores/corinos. We notice a curious correlation between the abundances of PO and PN and atomic nitrogen. Since N atoms are comparatively abundant in diffuse clouds and PDRs than in the hot core/corino region, PO/PN reflects < 1 in diffuse clouds, << 1 in PDRs, and > 1 in the late warm-up evolutionary phase of the hot core/corino regions. During the end of the post-warm-up phase, we obtain PO/PN > 1 for hot core and < 1 for its low mass analog. We employ a radiative transfer model to investigate the transitions of some of the P-bearing species in diffuse cloud and hot core regions and estimate the line profiles. Our study estimates the required integration time to observe these transitions with ground-based and space-based telescopes. We also carry out quantum chemical computation of the infrared features of PH3 along with various impurities. We notice that SO2 overlaps with the PH3 bending-scissoring modes around ~ (1000 - 1100) cm-1. We also find that the presence of CO2 can strongly influence the intensity of the stretching modes around ~ 2400 cm-1 of PH3 ., Comment: 44 pages, 28 figures, Accepted for the publication in The Astronomical Journal

Spectroscopic studies play a key role in the identification and analysis of interstellar ices and their structure. Some molecules have been identified within the interstellar ices either as pure, mixed, or even as layered structures. Absorption band features of water ice can significantly change with the presence of different types of impurities (CO, CO2, CH3OH, H2CO, etc.). In this work, we carried out a theoretical investigation to understand the behavior of water band frequency, and strength in the presence of impurities. The computational study has been supported and complemented by some infrared spectroscopy experiments aimed at verifying the effect of HCOOH, NH3 , and CH3 OH on the band profiles of pure H2O ice. Specifically, we explored the effect on the band strength of libration, bending, bulk stretching, and free-OH stretching modes. Computed band strength profiles have been compared with our new and existing experimental results, thus pointing out that vibrational modes of H2O and their intensities can change considerably in the presence of impurities at different concentrations. In most cases, the bulk stretching mode is the most affected vibration, while the bending is the least affected mode. HCOOH was found to have a strong influence on the libration, bending, and bulk stretching band profiles. In the case of NH3, the free-OH stretching band disappears when the impurity concentration becomes 50%. This work will ultimately aid a correct interpretation of future detailed spaceborne observations of interstellar ices by means of the upcoming JWST mission., Comment: Accepted for Publication in the ACS Earth and Space Chemistry Journal

Site-occupation embedding theory (SOET) is an in-principle-exact multi-determinantal extension of density-functional theory for model Hamiltonians. Various extensions of recent developments in SOET [Senjean et al., Phys. Rev. B 97, 235105 (2018)] are explored in this work. An important step forward is the generalization of the theory to multiple impurity sites. We also propose a new single-impurity density-functional approximation (DFA) where the density-functional impurity correlation energy of the two-level (2L) Hubbard system is combined with the Bethe ansatz local density approximation (BALDA) to the full correlation energy of the (infinite) Hubbard model. In order to test the new DFAs, the impurity-interacting wavefunction has been computed self-consistently with the density matrix renormalization group method (DMRG). Double occupation and per-site energy expressions have been derived and implemented in the one-dimensional case. A detailed analysis of the results is presented, with a particular focus on the errors induced either by the energy functionals solely or by the self-consistently converged densities. Among all the DFAs (including those previously proposed), the combined 2L-BALDA is the one that performs the best in all correlation and density regimes. Finally, extensions in new directions, like a partition-DFT-type reformulation of SOET, a projection-based SOET approach, or the combination of SOET with Green functions, are briefly discussed as a perspective., Comment: Version 1 contains errors in Appendix 6 (corresponding to the derivatives of the 2L-BALDA functional). Version 2 is correct

We propose a new simple computational model to estimate adsorption energies of atoms and molecules to low-temperature amorphous water ice, and we present the adsorption energies of carbon (3P), nitrogen (4S), and oxygen (3P) atoms based on quantum chemistry calculations. The adsorption energies were estimated to be 14100 +- 420 K for carbon, 400 +- 30 K for nitrogen, and 1440 +-160 K for oxygen. The adsorption energy of oxygen is well consistent with experimentally reported value. We found that the binding of a nitrogen atom is purely physisorption, while that of a carbon atom is chemisorption in which a chemical bond to an O atom of a water molecule is formed. That of an oxygen atom has a dual character both physisorption and chemisorption. The chemisorption of atomic carbon also implies a possibility of further chemical reactions to produce molecules bearing a C-O bond, while it may hinder the formation of methane on water ice via sequential hydrogenation of carbon atoms. These would be of a large impact to the chemical evolution of carbon species in interstellar environments. We also investigated effects of the newly calculated adsorption energies onto chemical compositions of cold dense molecular clouds with the aid of gas-ice astrochemical simulations. We found that abundances of major nitrogen-bearing molecules, such as N2 and NH3, are significantly altered by applying the calculated adsorption energy, because nitrogen atoms can thermally diffuse on surfaces even at 10 K., Comment: 14 pages, 9 figures, Accepted for publication in ApJ

Site-occupation embedding theory (SOET) is an alternative formulation of density-functional theory (DFT) for model Hamiltonians where the fully-interacting Hubbard problem is mapped, in principle exactly, onto an impurity-interacting (rather than a non-interacting) one. It provides a rigorous framework for combining wavefunction (or Green function) based methods with DFT. In this work, exact expressions for the per-site energy and double occupation of the uniform Hubbard model are derived in the context of SOET. As readily seen from these derivations, the so-called bath contribution to the per-site correlation energy is, in addition to the latter, the key density functional quantity to model in SOET. Various approximations based on Bethe ansatz and perturbative solutions to the Hubbard and single impurity Anderson models are constructed and tested on a one-dimensional ring. The self-consistent calculation of the embedded impurity wavefunction has been performed with the density matrix renormalization group method. It has been shown that promising results are obtained in specific regimes of correlation and density. Possible further developments have been proposed in order to provide reliable embedding functionals and potentials., Comment: Regular article with 14 pages including 6 figures

Current descriptions of the ab initio DMRG algorithm use two superficially different languages: an older language of the renormalization group and renormalized operators, and a more recent language of matrix product states and matrix product operators. The same algorithm can appear dramatically different when written in the two different vocabularies. In this work, we carefully describe the translation between the two languages in several contexts. First, we describe how to efficiently implement the ab-initio DMRG sweep using a matrix product operator based code, and the equivalence to the original renormalized operator implementation. Next we describe how to implement the general matrix product operator/matrix product state algebra within a pure renormalized operator-based DMRG code. Finally, we discuss two improvements of the ab initio DMRG sweep algorithm motivated by matrix product operator language: Hamiltonian compression, and a sum over operators representation that allows for perfect computational parallelism. The connections and correspondences described here serve to link the future developments with the past, and are important in the efficient implementation of continuing advances in ab initio DMRG and related algorithms., Comment: 35 pages, 10 figures

Linear response theory for the density matrix renormalization group (DMRG-LRT) was first presented in terms of the DMRG renormalization projectors [Dorando et al., J. Chem. Phys. 130, 184111 (2009)]. Later, with an understanding of the manifold structure of the matrix product state (MPS) ansatz, which lies at the basis of the DMRG algorithm, a way was found to construct the linear response space for general choices of the MPS gauge in terms of the tangent space vectors [Haegeman et al., Phys. Rev. Lett. 107, 070601 (2011)]. These two developments led to the formulation of the Tamm-Dancoff and random phase approximations (TDA and RPA) for MPS. This work describes how these LRTs may be efficiently implemented through minor modifications of the DMRG sweep algorithm, at a computational cost which scales the same as the ground-state DMRG algorithm. In fact, the mixed canonical MPS form implicit to the DMRG sweep is essential for efficient implementation of the RPA, due to the structure of the second-order tangent space. We present ab initio DMRG-TDA results for excited states of polyenes, the water molecule, and a [2Fe-2S] iron-sulfur cluster., Comment: 15 pages, 7 tables, 10 figures

The similarities between Hartree-Fock (HF) theory and the density-matrix renormalization group (DMRG) are explored. Both methods can be formulated as the variational optimization of a wave-function ansatz. Linearization of the time-dependent variational principle near a variational minimum allows to derive the random phase approximation (RPA). We show that the non-redundant parametrization of the matrix product state (MPS) tangent space [J. Haegeman et al., Phys. Rev. Lett. 107, 070601 (2011)] leads to the Thouless theorem for MPS, i.e. an explicit non-redundant parametrization of the entire MPS manifold, starting from a specific MPS reference. Excitation operators are identified, which extends the analogy between HF and DMRG to the Tamm-Dancoff approximation (TDA), the configuration interaction (CI) expansion, and coupled cluster theory. For a small one-dimensional Hubbard chain, we use a CI-MPS ansatz with single and double excitations to improve on the ground state and to calculate low-lying excitation energies. For a symmetry-broken ground state of this model, we show that RPA-MPS allows to retrieve the Goldstone mode. We also discuss calculations of the RPA-MPS correlation energy. With the long-range quantum chemical Pariser-Parr-Pople Hamiltonian, low-lying TDA-MPS and RPA-MPS excitation energies for polyenes are obtained., Comment: 16 pages, 3 figures and 1 table

We investigate tree tensor network states for quantum chemistry. Tree tensor network states represent one of the simplest generalizations of matrix product states and the density matrix renormalization group. While matrix product states encode a one-dimensional entanglement structure, tree tensor network states encode a tree entanglement structure, allowing for a more flexible description of general molecules. We describe an optimal tree tensor network state algorithm for quantum chemistry. We introduce the concept of half-renormalization which greatly improves the efficiency of the calculations. Using our efficient formulation we demonstrate the strengths and weaknesses of tree tensor network states versus matrix product states. We carry out benchmark calculations both on tree systems (hydrogen trees and \pi-conjugated dendrimers) as well as non-tree molecules (hydrogen chains, nitrogen dimer, and chromium dimer). In general, tree tensor network states require much fewer renormalized states to achieve the same accuracy as matrix product states. In non-tree molecules, whether this translates into a computational savings is system dependent, due to the higher prefactor and computational scaling associated with tree algorithms. In tree like molecules, tree network states are easily superior to matrix product states. As an ilustration, our largest dendrimer calculation with tree tensor network states correlates 110 electrons in 110 active orbitals., Comment: 15 pages, 19 figures