-
DMRG-tailored coupled cluster method in the 4c-relativistic domain: General implementation and application to the NUHFI and NUF$_3$ molecules
Authors:
Jakub Višňák,
Jan Brandejs,
Mihály Máté,
Lucas Visscher,
Örs Legeza,
Jiří Pittner
Abstract:
Heavy atom compounds represent a challenge for computational chemistry, due to the need for simultaneous treatment of relativistic and correlation effects. Often such systems exhibit also strong correlation which hampers the application of perturbation theory or single-reference coupled cluster (CC) methods. As a viable alternative, we have proposed to externally correct the CC method using the de…
▽ More
Heavy atom compounds represent a challenge for computational chemistry, due to the need for simultaneous treatment of relativistic and correlation effects. Often such systems exhibit also strong correlation which hampers the application of perturbation theory or single-reference coupled cluster (CC) methods. As a viable alternative, we have proposed to externally correct the CC method using the density matrix renormalization group (DMRG) wave functions, yielding the DMRG-tailored CC method. In a previous paper [J. Chem. Phys. {\bf 152}, 174107 (2020)] we have reported a first implementation of this method in the relativistic context, which was restricted to molecules with real double group symmetry. In this work we present a fully general implementation of the method, covering complex and quaternion double groups as well. The 4c-TCC method thus becomes applicable to polyatomic molecules including heavy atoms. For assessment of the method, we performed calculations of the chiral uranium compound NUHFI, which was previously studied in the context of the enhancement of parity violation effects. In particular, we performed calculations of a cut of the potential energy surface of this molecule along the dissociation of the N-U bond, where the system exhibits a strong multireference character. Since there are no experimental data for NUHFI, we have performed also an analogous study of the (more symmetric) NUF$_3$ molecule, where the vibrational frequency of the N-U bond can be compared with spectroscopic data.
△ Less
Submitted 30 April, 2024; v1 submitted 27 March, 2024;
originally announced March 2024.
-
Compressing multireference character of wave functions via fermionic mode optimization
Authors:
Mihály Máté,
Klára Petrov,
Szilárd Szalay,
Örs Legeza
Abstract:
In this work, we present a brief overview of the fermionic mode optimization within the framework of tensor network state methods [PRL 117,210402(2016)], and demonstrate that it has the potential to compress the multireference character of the wave functions after finding optimal molecular orbitals (modes), based on entanglement minimization. Numerical simulations have been performed for the nitro…
▽ More
In this work, we present a brief overview of the fermionic mode optimization within the framework of tensor network state methods [PRL 117,210402(2016)], and demonstrate that it has the potential to compress the multireference character of the wave functions after finding optimal molecular orbitals (modes), based on entanglement minimization. Numerical simulations have been performed for the nitrogen dimer in the cc-pVDZ basis for the equilibrium and for stretched geometries.
△ Less
Submitted 28 April, 2022;
originally announced April 2022.
-
Assessing the accuracy of tailored coupled cluster methods corrected by electronic wave functions of polynomial cost
Authors:
Aleksandra Leszczyk,
Mihály Máté,
Örs Legeza,
Katharina Boguslawski
Abstract:
Tailored coupled cluster theory represents a computationally inexpensive way to describe static and dynamical electron correlation effects. In this work, we scrutinize the performance of various tailored coupled cluster methods externally corrected by electronic wave functions of polynomial cost. Specifically, we focus on frozen-pair coupled-cluster (fpCC) methods, which are externally corrected b…
▽ More
Tailored coupled cluster theory represents a computationally inexpensive way to describe static and dynamical electron correlation effects. In this work, we scrutinize the performance of various tailored coupled cluster methods externally corrected by electronic wave functions of polynomial cost. Specifically, we focus on frozen-pair coupled-cluster (fpCC) methods, which are externally corrected by pair-coupled cluster doubles (pCCD), and coupled cluster theory tailored by matrix product state wave functions optimized by the density matrix renormalization group (DMRG) algorithm. As test system, we selected a set of various small- and medium-sized molecules containing diatomics (N$_2$, F$_2$, C$_2$, CN$^+$, BN, BO$^+$, and Cr$_2$) and molecules (ammonia, ethylene, cyclobutadiene, benzene) for which conventional single-reference coupled cluster singles and doubles (CCSD) is not able to produce accurate results for spectroscopic constants, potential energy surfaces, and barrier heights. Most importantly, DMRG-tailored and pCCD-tailored approaches yield similar errors in spectroscopic constants and potential energy surfaces compared to multireference and/or experimental reference data and generally outrank the conventional single-reference CCSD approach. Although fpCC methods provide a reliable description for the dissociation pathway of molecules featuring single and quadruple bonds, they fail in the description of triple or hextuple bond-breaking processes or avoided crossing regions.
△ Less
Submitted 23 March, 2021;
originally announced March 2021.
-
Fermionic systems for quantum information people
Authors:
Szilárd Szalay,
Zoltán Zimborás,
Mihály Máté,
Gergely Barcza,
Christian Schilling,
Örs Legeza
Abstract:
The operator algebra of fermionic modes is isomorphic to that of qubits, the difference between them is twofold: the embedding of subalgebras corresponding to mode subsets and multiqubit subsystems on the one hand, and the parity superselection in the fermionic case on the other. We discuss these two fundamental differences extensively, and illustrate these through the Jordan--Wigner representatio…
▽ More
The operator algebra of fermionic modes is isomorphic to that of qubits, the difference between them is twofold: the embedding of subalgebras corresponding to mode subsets and multiqubit subsystems on the one hand, and the parity superselection in the fermionic case on the other. We discuss these two fundamental differences extensively, and illustrate these through the Jordan--Wigner representation in a coherent, self-contained, pedagogical way, from the point of view of quantum information theory. Our perspective leads us to develop useful new tools for the treatment of fermionic systems, such as the fermionic (quasi-)tensor product, fermionic canonical embedding, fermionic partial trace, fermionic products of maps and fermionic embeddings of maps. We formulate these by direct, easily applicable formulas, without mode permutations, for arbitrary partitionings of the modes. It is also shown that fermionic reduced states can be calculated by the fermionic partial trace, containing the proper phase factors. We also consider variants of the notions of fermionic mode correlation and entanglement, which can be endowed with the usual, local operation based motivation, if the parity superselection rule is imposed. We also elucidate some other fundamental points, related to joint map extensions, which make the parity superselection inevitable in the description of fermionic systems.
△ Less
Submitted 10 December, 2021; v1 submitted 4 June, 2020;
originally announced June 2020.
-
How creating one additional well can generate Bose-Einstein condensation
Authors:
Mihály Máté,
Örs Legeza,
Rolf Schilling,
Mason Yousif,
Christian Schilling
Abstract:
The realization of Bose-Einstein condensation in ultracold trapped gases has led to a revival of interest in that fascinating quantum phenomenon. This experimental achievement necessitated both extremely low temperatures and sufficiently weak interactions. Particularly in reduced spatial dimensionality even an infinitesimal interaction immediately leads to a departure to quasi-condensation. We pro…
▽ More
The realization of Bose-Einstein condensation in ultracold trapped gases has led to a revival of interest in that fascinating quantum phenomenon. This experimental achievement necessitated both extremely low temperatures and sufficiently weak interactions. Particularly in reduced spatial dimensionality even an infinitesimal interaction immediately leads to a departure to quasi-condensation. We propose a system of strongly interacting bosons which overcomes those obstacles by exhibiting a number of intriguing related features: (i) The tuning of just a single control parameter drives a transition from quasi-condensation to complete condensation, (ii) the destructive influence of strong interactions is compensated by the respective increased mobility, (iii) topology plays a crucial role since a crossover from one- to `infinite'-dimensionality is simulated, (iv) a ground state gap opens which makes the condensation robust to thermal noise. Remarkably, all these features can be derived by analytical and exact numerical means despite the non-perturbative character of the system.
△ Less
Submitted 19 February, 2021; v1 submitted 23 February, 2020;
originally announced February 2020.
-
Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer
Authors:
Fabian M. Faulstich,
Mihály Máté,
Andre Laestadius,
Mihály András Csirik,
Libor Veis,
Andrej Antalik,
Jiří Brabec,
Reinhold Schneider,
Jiří Pittner,
Simen Kvaal,
Örs Legeza
Abstract:
In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate chemical properties of the used method, such as energy size extensivity and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use…
▽ More
In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate chemical properties of the used method, such as energy size extensivity and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of a so-called CAS-ext gap describing the basis splitting between the complete active space and the external part. Moreover, the behavior of the energy error as a function of the optimal basis splitting is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose the complete active space. Furthermore, we extend the mathematical analysis with a numerical study on the complete active space dependence of the error.
△ Less
Submitted 20 September, 2018;
originally announced September 2018.
-
Bit Patterned Magnetic Recording: Theory, Media Fabrication, and Recording Performance
Authors:
Thomas R. Albrecht,
Hitesh Arora,
Vipin Ayanoor-Vitikkate,
Jean-Marc Beaujour,
Daniel Bedau,
David Berman,
Alexei L. Bogdanov,
Yves-Andre Chapuis,
Julia Cushen,
Elizabeth E. Dobisz,
Gregory Doerk,
He Gao,
Michael Grobis,
Bruce Gurney,
Weldon Hanson,
Olav Hellwig,
Toshiki Hirano,
Pierre-Olivier Jubert,
Dan Kercher,
Jeffrey Lille,
Zuwei Liu,
C. Mathew Mate,
Yuri Obukhov,
Kanaiyalal C. Patel,
Kurt Rubin
, et al. (6 additional authors not shown)
Abstract:
Bit Patterned Media (BPM) for magnetic recording provide a route to densities $>1 Tb/in^2$ and circumvents many of the challenges associated with conventional granular media technology. Instead of recording a bit on an ensemble of random grains, BPM uses an array of lithographically defined isolated magnetic islands, each of which stores one bit. Fabrication of BPM is viewed as the greatest challe…
▽ More
Bit Patterned Media (BPM) for magnetic recording provide a route to densities $>1 Tb/in^2$ and circumvents many of the challenges associated with conventional granular media technology. Instead of recording a bit on an ensemble of random grains, BPM uses an array of lithographically defined isolated magnetic islands, each of which stores one bit. Fabrication of BPM is viewed as the greatest challenge for its commercialization. In this article we describe a BPM fabrication method which combines e-beam lithography, directed self-assembly of block copolymers, self-aligned double patterning, nanoimprint lithography, and ion milling to generate BPM based on CoCrPt alloys. This combination of fabrication technologies achieves feature sizes of $<10 nm$, significantly smaller than what conventional semiconductor nanofabrication methods can achieve. In contrast to earlier work which used hexagonal close-packed arrays of round islands, our latest approach creates BPM with rectangular bitcells, which are advantageous for integration with existing hard disk drive technology. The advantages of rectangular bits are analyzed from a theoretical and modeling point of view, and system integration requirements such as servo patterns, implementation of write synchronization, and providing for a stable head-disk interface are addressed in the context of experimental results. Optimization of magnetic alloy materials for thermal stability, writeability, and switching field distribution is discussed, and a new method for growing BPM islands on a patterned template is presented. New recording results at $1.6 Td/in^2$ (teradot/inch${}^2$, roughly equivalent to $1.3 Tb/in^2$) demonstrate a raw error rate $<10^{-2}$, which is consistent with the recording system requirements of modern hard drives. Extendibility of BPM to higher densities, and its eventual combination with energy assisted recording are explored.
△ Less
Submitted 19 March, 2015;
originally announced March 2015.