BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Henrik Koch

Henrik Koch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phcond-mat.mtrl-sciquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
3topics
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

13 published item(s)

preprint2026arXiv

Polaritonic Bloch's Theorem beyond the Long-Wavelength Approximation

Cavity quantum electrodynamics provides a powerful tool to manipulate material properties, yet it remains a matter of debate whether and how quantized fields affect the periodicity of crystals. Here, we extend Bloch's theorem to crystals under strong light-matter coupling, revealing that polariton quasiparticles preserve lattice periodicity. We introduce a general framework to incorporate multimode cavity fields in a simple and tractable way, showing that additional modes contribute small energy corrections noticeable only at low frequencies. Within the single-photon approximation, these contributions reduce to a spatially uniform effective field in the crystal plane, providing a formal justification for the single-mode and long-wavelength approximations commonly used in molecular polaritonics. Together, these results establish a rigorous framework for describing polaritonic states in crystalline solids.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Efficient implementation of molecular CCSD gradients with Cholesky-decomposed electron repulsion integrals

We present an efficient implementation of ground and excited state CCSD gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition, like density-fitting, is an inner projection method, and thus similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. The computational time spent per optimization cycle are furthermore found to be significantly lower compared to other implementations based on an inner projection method. We illustrate the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the characteristic features of ionization in QED environments

The ionization of molecular systems is important in many chemical processes such as electron transfer and hot electron injection. Strong coupling between molecules and quantized fields (e.g. inside optical cavities) represents a new promising way to modify molecular properties in a non-invasive way. Recently, strong light-matter coupling has shown the potential to significantly improve the rates of hot electron driven processes, for instance in water splitting. In this paper, we demonstrate that inside an optical cavity the residual interaction between an outgoing free electron and the vacuum field is significant. We further show that, since the quantized field is also interacting with the ionized molecule, the free electron and the molecular system are correlated. We develop a theoretical framework to account for the field induced correlation and show that the interaction between the free electron and the field free electron-field interaction has sizeable effects on the ionization potential of typical organic molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Oscillator strengths in the framework of equation of motion multilevel CC3

We present an efficient implementation of the equation of motion oscillator strengths for the closed-shell multilevel coupled cluster singles and doubles with perturbative triples method (MLCC3) in the electronic structure program eT. The orbital space is split into an active part treated with CC3 and an inactive part computed at the coupled cluster singles and doubles (CCSD) level of theory. Asymptotically, the CC3 contribution scales as $O(n_\text{V} n^3_\text{v} n^3_\text{o})$ floating-point operations (FLOP), where $n_V$ is the total number of virtual orbitals while $n_\text{v}$ and $n_\text{o}$ are the number of active virtual and occupied orbitals, respectively. The CC3 contribution, thus, only scales linearly with the full system size and can become negligible compared to the cost of CCSD. We demonstrate the capabilities of our implementation by calculating the UV-VIS spectrum of azobenzene and a core excited state of betaine 30 with more than 1000 molecular orbitals.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Molecular orbital theory in cavity QED environments

Coupling between molecules and vacuum photon fields inside an optical cavity has proven to be an effective way to engineer molecular properties, in particular reactivity. To ease the rationalization of cavity induced effects we introduce an ab initio method leading to the first fully consistent molecular orbital theory for quantum electrodynamics environments. Our framework is non-perturbative and explains modifications of the electronic structure due to the interaction with the photon field. We show that the newly developed orbital theory can be used to predict cavity induced modifications of molecular reactivity and pinpoint classes of systems with significant cavity effects. We also investigate cavity-induced modifications of molecular reactivity in the vibrational strong coupling regime.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Simulating weak-field attosecond processes with a Lanczos reduced basis approach to time-dependent equation of motion coupled cluster theory

A time-dependent equation of motion coupled cluster singles and doubles (TD-EOM-CCSD) method is implemented, which uses a reduced basis calculated with the asymmetric band Lanczos algorithm. The approach is used to study weak-field processes in small molecules induced by ultrashort valence pump and core probe pulses. We assess the reliability of the procedure by comparing TD-EOM-CCSD absorption spectra to spectra obtained from the time-dependent coupled-cluster singles and doubles (TDCCSD) method and observe that spectral features can be reproduced for several molecules, at much lower computational times. We discuss how multiphoton absorption and symmetry can be handled in the method and general features of the core-valence separation (CVS) projection technique. We also model the transient absorption of an attosecond X-ray probe pulse by the glycine molecule.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A biorthonormal formalism for nonadiabatic coupled cluster dynamics

In coupled cluster methods, the electronic states are biorthonormal in the sense that the left states are orthonormal to the right states. Here we present an extension of this formalism to a left and right total molecular wave function. Starting from left and right Born-Huang expansions, we derive projected Schrödinger equations for the left and right nuclear wave functions. Observables may be extracted from the resulting wave function pair using standard expressions. The formalism is shown to be invariant under electronic basis transformations, such as normalization of the electronic states. Consequently, the nonadiabatic coupling elements can be expressed with biorthonormal wave functions. Calculating normalization factors that scale as full-CI is therefore not necessary, contrary to claims in the literature. For nuclear dynamics, we therefore need expressions for the vector and scalar couplings in the biorthonormal formalism. We derive these expressions using a Lagrangian formalism.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accelerated multimodel Newton-type algorithms for faster convergence of ground and excited state coupled cluster equations

We introduce a multimodel approach to solve coupled cluster equations, employing a quasi Newton algorithm for the ground state and an Olsen algorithm for the excited states. In these algorithms, both of which can be viewed as Newton algorithms, the Jacobian matrix of a lower level coupled cluster model is used in Newton equations associated with the target model. Improvements in convergence then implies savings for sufficiently large molecular systems, since the computational cost of macroiterations scales more steeply with system size than the cost of microiterations. The multimodel approach is suitable when there is a lower level Jacobian matrix that is much more accurate than the zeroth order approximation. Applying the approach to the CC3 equations, using the CCSD approximation of the Jacobian, we show that the time spent to determine the ground and valence excited states can be significantly reduced. We also find improved convergence for core excited states, indicating that similar savings will be obtained with an explicit implementation of the core-valence separated CCSD Jacobian transformation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy-Based Molecular Orbital Localization in a Specific Spatial Region

We present a novel energy-based localization procedure able to localize molecular orbitals into specific spatial regions. The method is applied to several cases including both conjugated and non-conjugated systems. The obtained localized molecular orbitals are used in a multiscale framework based on the multilevel Hartree-Fock approach. An almost perfect agreement with reference values is achieved for both ground state properties, such as dipole moments, and local excitation energies calculated at the coupled cluster level. The proposed approach is useful to extend the application range of high level electron correlation methods. In fact, the reduced number of molecular orbitals can lead to a large reduction in the computational cost of correlated calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Equation-of-motion MLCCSD and CCSD-in-HF oscillator strengths and their application to core excitations

We present an implementation of equation-of-motion oscillator strengths for the multilevel CCSD (MLCCSD) model where CCS is used as the lower level method (CCS/CCSD). In this model, the double excitations of the cluster operator are restricted to an active orbital space, whereas the single excitations are unrestricted. Calculated nitrogen K-edge spectra of adenosine, adenosine triphosphate (ATP), and an ATP-water system are used to demonstrate the performance of the model. Projected atomic orbitals (PAOs) are used to partition the virtual space into active and inactive orbital sets. Cholesky decomposition of the Hartree-Fock density is used to partition the occupied orbitals. This Cholesky-PAO partitioning is cheap, scaling as $\mathcal{O}(N^3)$, and is suitable for the calculation of core excitations which are localized in character. By restricting the single excitations of the cluster operator to the active space, as well as the double excitations, the CCSD-in-HF model is obtained. A comparison of the two models---MLCCSD and CCSD-in-HF---is presented for the core excitation spectra of the adenosine and ATP systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

eT 1.0: an open source electronic structure program with emphasis on coupled cluster and multilevel methods

The eT program is an open source electronic structure package with emphasis on coupled cluster and multilevel methods. It includes efficient spin adapted implementations of ground and excited singlet states, as well as equation of motion oscillator strengths, for CCS, CC2, CCSD, and CC3. Furthermore, eT provides unique capabilities such as multilevel Hartree-Fock and multilevel CC2, real-time propagation for CCS and CCSD, and efficient CC3 oscillator strengths. With a coupled cluster code based on an efficient Cholesky decomposition algorithm for the electronic repulsion integrals, eT has similar advantages as codes using density fitting, but with strict error control. Here we present the main features of the program and demonstrate its performance through example calculations. Because of its availability, performance, and unique capabilities, we expect eT to become a valuable resource to the electronic structure community.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multilevel Density Functional Theory

We introduce a novel density-based multilevel approach in density functional theory. In this multilevel density functional theory (MLDFT), the system is partitioned in an active and an inactive fragment, and all interactions are retained between the two parts. In MLDFT, the Kohn-Sham equations are solved in the MO basis for the active part only, while keeping the inactive density frozen. This results in a reduction of computational cost. We outline the theory and implementation, and discuss applications to aqueous solutions of methyloxirane and glycidol.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time-dependent coupled cluster theory for ultrafast transient absorption spectroscopy

We present a spin-adapted time-dependent coupled cluster singles and doubles model for the molecular response to a sequence of ultrashort laser pulses. The implementation is used to calculate the electronic response to a valence-exciting pump pulse, and a subsequent core-exciting probe pulse. We assess the accuracy of the integration procedures used in solving the dynamic coupled cluster equations, in order to find a compromise between computational cost and accuracy. The transient absorption spectrum of lithium fluoride is calculated for various delays of the probe pulse with respect to the pump pulse. We observe that the transient probe absorption oscillates with the pump-probe delay, an effect that is attributed to the interference of states in the pump-induced superposition.

0 citations0 reviewsNo code linkedOpen access