Researcher profile

Lucas Visscher

Lucas Visscher contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 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.