Source author record

Lucas Visscher

Lucas Visscher appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the number of correlated electrons and the size of the basis set is a bottleneck which limits applications of coupled cluster (CC) algorithms. This is particularly so for calculations based on 2- or 4-component relativistic Hamiltonians, which often employ uncontracted basis sets and lead to large virtual molecular orbital (VMO) spaces. This problem may be alleviated by employing a more compact set of virtual spinors than those provided by the canonical Hartree-Fock (HF) set, such as those based on natural orbitals (NOs). In this paper we describe the implementation of a module for generating NOs for correlated wavefunctions, and in particular MP2 frozen natural orbitals (MP2FNOs), as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [J. Pototschnig et. al., J. Chem. Theory Comput. 17, 5509, 2021]. Our implementation is capable of manipulating both complex and quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, so that the code also makes it possible to generate CCSD natural orbitals in AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties at CCSD level such as the electric field gradients at the nuclei (EFGs), electric dipole and quadrupole moments for hydrogen halides HX (X=F-Ts), and parity-violating energy differences (PV) for H$_2$Y$_2$ (Y=O-Se), we find that MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the periodic table and that, with VMO spaces truncated to around half the size of the full spaces ones, it is possible to obtain reliable estimates for both energies and all molecular properties considered.

preprint2022arXiv

Exploring the Statically Screened $G3W2$ Correction to the $GW$ Self-Energy: Charged Excitations and Total Energies of Finite Systems

Electron correlation in finite and extended systems is often described in an effective single-particle framework within the $GW$ approximation. Here, we use the statically screened second-order exchange contribution to the self-energy ($G3W2$) to calculate a perturbative correction to the $GW$ self-energy. We use this correction to calculate total correlation energies of atoms, relative energies, as well as charged excitations of a wide range of molecular systems. We show that the second-order correction improves correlation energies with respect to the RPA and also improves relative energies for many, but not all considered systems. While the full $G3W2$ contribution does not give consistent improvements over $GW$, taking the average of $GW$ and $GW + G3W2$ generally gives excellent results. Improvements over quasiparticle self-consistent $GW$, which we show to give very accurate charged excitations in small and medium molecules by itself, are only minor. $G_0W_0$ quasiparticle energies evaluated with eigenvalue and orbitals from range-separated hybrids, however, are tremendously improved upon: The second-order corrected $G_0W_0$ outperforms all existing $GW$ methods for the systems considered herein and also does not come with substantially increased computational cost compared to $G_0W_0$ for systems with up to 100 atoms.

preprint2022arXiv

Quasiparticle Self-Consistent $GW$-Bethe-Salpeter equation calculations for large chromophoric systems

The $GW$-Bethe-Salpeter Equation (BSE) method is promising for calculating the low-lying excited states of molecular systems. So far, it has only been applied to rather small molecules, and in the commonly implemented diagonal approximations to the electronic self-energy it depends on a mean-field starting point. We describe here an implementation of the self-consistent and starting-point independent quasiparticle self-consistent (qs$GW$)-BSE approach which is suitable for calculations on large molecules. We herein show that eigenvalue-only self-consistency leads to an unfaithful description of certain excitonic states for Chlorophyll dimers while the qs$GW$-BSE vertical excitation energies (VEE) are in excellent agreement with spectroscopic experiments for Chlorophyll monomers and dimers measured in the gas phase. On the other hand, VEEs from time-dependent density functional theory calculations tend to disagree with experimental values and using different range-separated hybrid (RSH) kernels changes the VEEs by up to 0.5 eV. We use the new qs$GW$-BSE implementation to calculate the lowest excitation energies of the six chromophores of the photosystem II (PSII) reaction center (RC) with nearly 2000 correlated electrons. Using more than 11000 (6000) basis functions, the calculation could be completed in less than 5 (2) days one a single modern compute node. In agreement with previous TD-DFT calculations using RSH kernels on models that do also not include environment effects, our qs$GW$-BSE calculations only yield states with local character in the low-energy spectrum of the hexameric complex. Earlier work with RSH kernels has demonstrated that the protein environment facilitates the experimentally observed interchromophoric charge transfer. Therefore, future research will need to combine correlation effects beyond TD-DFT with an explicit treatment of environment electrostatics.

preprint2021arXiv

Generalization of intrinsic orbitals to Kramers-paired quaternion spinors, molecular fragments and valence virtual spinors

Localization of molecular orbitals finds its importance in the representation of chemical bonding (and anti-bonding) and in the local correlation treatments beyond mean-field approximation. In this paper, we generalize the intrinsic atomic and bonding orbitals [G. Knizia, J. Chem. Theory Comput. 2013, 9, 11, 4834-4843] to relativistic applications using complex and quaternion spinors, as well as to molecular fragments instead of atomic fragments only. By performing a singular value decomposition, we show how localized valence virtual orbitals can be expressed in this intrinsic minimal basis. We demonstrate our method on systems of increasing complexity, starting from simple cases such as benzene, acrylic-acid and ferrocene molecules, and then demonstrating the use of molecular fragments and inclusion of relativistic effects for complexes containing heavy elements such as tellurium, iridium and astatine. The aforementioned scheme is implemented into a standalone program interfaced with several different quantum chemistry packages.

preprint2020arXiv

Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets

The superposition of atomic potentials (SAP) approach has recently been shown to be a simple and efficient way to initialize electronic structure calculations [S. Lehtola, J. Chem. Theory Comput. 15, 1593 (2019)]. Here, we study the differences between effective potentials from fully numerical density functional and optimized effective potential calculations for fixed configurations. We find that the differences are small, overall, and choose exchange-only potentials at the local density approximation level of theory computed on top of Hartree-Fock densities as a good compromise. The differences between potentials arising from different atomic configurations are also found to be small at this level of theory. Furthermore, we discuss the efficient Gaussian-basis implementation of SAP via error function fits to fully numerical atomic radial potentials. The guess obtained from the fitted potentials can be easily implemented in any Gaussian-basis quantum chemistry code in terms of two-electron integrals. Fits covering the whole periodic table from H to Og are reported for non-relativistic as well as fully relativistic four-component calculations that have been carried out with fully numerical approaches.

preprint2016arXiv

Tight-Binding Approximations to Time-Dependent Density Functional Theory - a fast approach for the calculation of electronically excited states

We propose a new method of calculating electronically excited states that combines a density functional theory (DFT) based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive time-dependent density functional theory (TD-DFT) calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB.

preprint2016arXiv

Vibrationally resolved UV/Vis spectroscopy with time-dependent density functional based tight binding

We report a time-dependent density functional based tight-binding (TD-DFTB) scheme for the calculation of UV/Vis spectra, explicitly taking into account the excitation of nuclear vibrations via the adiabatic Hessian Franck-Condon (AH|FC) method with a harmonic approximation for the nuclear wavefunction. The theory of vibrationally resolved UV/Vis spectroscopy is first summarized from the viewpoint of TD-DFTB. The method is benchmarked against time-dependent density functional theory (TD-DFT) calculations for strongly dipole allowed excitations in various aromatic and polar molecules. Using the recent 3ob:freq parameter set of Elstner's group, very good agreement with TD-DFT calculations using local functionals was achieved.

preprint2014arXiv

Efficient Calculation of Electronic Absorption Spectra by Means of Intensity-Selected TD-DFTB

During the last two decades density functional based linear response approaches have become the de facto standard for the calculation of optical properties of small and medium-sized molecules. At the heart of these methods is the solution of an eigenvalue equation in the space of single-orbital transitions, whose quickly increasing number makes such calculations costly if not infeasible for larger molecules. This is especially true for time-dependent density functional tight binding (TD-DFTB), where the evaluation of the matrix elements is inexpensive. For the relatively large systems that can be studied the solution of the eigenvalue equation therefore determines the cost of the calculation. We propose to do an oscillator strength based truncation of the single-orbital transition space to reduce the computational effort of TD-DFTB based absorption spectra calculations. We show that even a sizeable truncation does not destroy the principal features of the absorption spectrum, while naturally avoiding the unnecessary calculation of excitations with small oscillator strengths. We argue that the reduced computational cost of intensity-selected TD-DFTB together with its ease of use compared to other methods lowers the barrier of performing optical properties calculations of large molecules, and can serve to make such calculations possible in a wider array of applications.

preprint2013arXiv

An Accurate and Linear Scaling Method to Calculate Charge-Transfer Excitation Energies and Diabatic Couplings

Quantum--Mechanical methods that are both computationally fast and accurate are not yet available for electronic excitations having charge transfer character. In this work, we present a significant step forward towards this goal for those charge transfer excitations that take place between non-covalently bound molecules. In particular, we present a method that scales linearly with the number of non-covalently bound molecules in the system and is based on a two-pronged approach: The molecular electronic structure of broken-symmetry charge-localized states is obtained with the Frozen Density Embedding formulation of subsystem Density-Functional Theory; subsequently, in a post-SCF calculation, the full-electron Hamiltonian and overlap matrix elements among the charge-localized states are evaluated with an algorithm which takes full advantage of the subsystem DFT density partitioning technique. The method is benchmarked against Coupled-Cluster calculations and achieves chemical accuracy for the systems considered for intermolecular separations ranging from hydrogen-bond distances to tents of Ångstroms. Numerical examples are provided for molecular clusters comprised of up to 56 non-covalently bound molecules.

preprint2012arXiv

Relativistic quantum chemistry on quantum computers

Last years witnessed a remarkable interest in application of quantum computing for solving problems in quantum chemistry more efficiently than classical computers allow. Very recently, even first proof-of-principle experimental realizations have been reported. However, so far only the non-relativistic regime (i.e. Schroedinger equation) has been explored, while it is well known that relativistic effects can be very important in chemistry. In this letter we present the first quantum algorithm for relativistic computations of molecular energies. We show how to efficiently solve the eigenproblem of the Dirac-Coulomb Hamiltonian on a quantum computer and demonstrate the functionality of the proposed procedure by numerical simulations of computations of the spin-orbit splitting in the SbH molecule. Finally, we propose quantum circuits with 3 qubits and 9 or 10 CNOTs, which implement a proof-of-principle relativistic quantum chemical calculation for this molecule and might be suitable for an experimental realization.

preprint2010arXiv

Relativistic general-order coupled-cluster method for high-precision calculations: Application to Al+ atomic clock

We report the implementation of a general-order relativistic coupled-cluster method for performing high-precision calculations of atomic and molecular properties. As a first application, the static dipole polarizabilities of the ground and first excited states of Al+ have been determined to precisely estimate the uncertainty associated with the BBR shift of its clock frequency measurement. The obtained relative BBR shift is -3.66+-0.44 for the 3s^2 ^1S_0^0 --> 3s3p ^3P_0^0 transition in Al+ in contrast to the value obtained in the latest clock frequency measurement, -9+-3 [Phys. Rev. Lett. 104, 070802 (2010)]. The method developed in the present work can be employed to study a variety of subtle effects such as fundamental symmetry violations in atoms.