Showing 1–50 of 222 results for author: Koch, H

Towards polaritonic molecular orbitals for large molecular systems

Authors: Yassir El Moutaoukal, Rosario R. Riso, Matteo Castagnola, Henrik Koch

Abstract: A comprehensive theoretical understanding of electron-photon correlation is essential for describing the reshaping of molecular orbitals in quantum electrodynamics (QED) environments. The strong coupling QED Hartree-Fock (SC-QED-HF) theory tackles these aspects by providing consistent molecular orbitals in the strong coupling regime. The previous implementation, however, has significant convergenc… ▽ More A comprehensive theoretical understanding of electron-photon correlation is essential for describing the reshaping of molecular orbitals in quantum electrodynamics (QED) environments. The strong coupling QED Hartree-Fock (SC-QED-HF) theory tackles these aspects by providing consistent molecular orbitals in the strong coupling regime. The previous implementation, however, has significant convergence issues that limit the applicability. In this work we introduce two second-order algorithms that significantly reduce the computational requirements, thereby enhancing the modeling of large molecular systems in QED environments. Furthermore, the implementation will enable the development of correlated methods based on a reliable molecular orbital framework. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 29 pages and 6 figures

Analytical evaluation of ground state gradients in quantum electrodynamics coupled cluster theory

Authors: Marcus Takvam Lexander, Sara Angelico, Eirik Fadum Kjønstad, Henrik Koch

Abstract: Analytical gradients of potential energy surfaces play a central role in quantum chemistry, allowing for molecular geometry optimizations and molecular dynamics simulations. In strong coupling conditions, potential energy surfaces can account for strong interactions between matter and the quantized electromagnetic field. In this paper, we derive expressions for the ground state analytical gradient… ▽ More Analytical gradients of potential energy surfaces play a central role in quantum chemistry, allowing for molecular geometry optimizations and molecular dynamics simulations. In strong coupling conditions, potential energy surfaces can account for strong interactions between matter and the quantized electromagnetic field. In this paper, we derive expressions for the ground state analytical gradients in quantum electrodynamics coupled cluster theory. We also present a Cholesky-based implementation for the coupled cluster singles and doubles model. We report timings to show the performance of the implementation and present optimized geometries to highlight cavity-induced molecular orientation effects in strong coupling conditions. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 45 pages and 3 figures

Almost surely smoothed scattering for cubic NLS

Authors: Nicolas Burq, Herbert Koch, Nicola Visciglia, Nikolay Tzvetkov

Abstract: We consider cubic NLS in dimensions 2, 3, 4 and we prove that almost surely solutions with randomized initial data at low regularity scatter. Moreover, we establish some smoothing properties of the associated scattering operator and precise the rate of convergence. We consider cubic NLS in dimensions 2, 3, 4 and we prove that almost surely solutions with randomized initial data at low regularity scatter. Moreover, we establish some smoothing properties of the associated scattering operator and precise the rate of convergence. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 23 pages

MSC Class: 35Q55; 35R60

Strong coupling electron-photon dynamics: a real-time investigation of energy redistribution in molecular polaritons

Authors: Matteo Castagnola, Marcus T. Lexander, Enrico Ronca, Henrik Koch

Abstract: We analyze the real-time electron-photon dynamics in long-range polariton-mediated energy transfer using a real-time quantum electrodynamics coupled cluster (RT-QED-CC) model, which allows for spatial and temporal visualization of transport processes. We compute the time evolution of photonic and molecular observables, such as the dipole moment and the photon coordinate, following the excitation o… ▽ More We analyze the real-time electron-photon dynamics in long-range polariton-mediated energy transfer using a real-time quantum electrodynamics coupled cluster (RT-QED-CC) model, which allows for spatial and temporal visualization of transport processes. We compute the time evolution of photonic and molecular observables, such as the dipole moment and the photon coordinate, following the excitation of the system induced by short laser pulses. Our simulation highlights the different time scales of electrons and photons under light-matter strong coupling, the role of dark states, and the differences with the electronic (Förster and Dexter) energy exchange mechanisms. The developed method can simulate multiple high-intensity laser pulses while explicitly retaining electronic and electron-photon correlation and is thus suited for nonlinear optics and transient absorption spectroscopies of molecular polaritons. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 19 pages and 9 figures

Unexpected hydrogen dissociation in thymine: predictions from a novel coupled cluster theory

Authors: Eirik F. Kjønstad, O. Jonathan Fajen, Alexander C. Paul, Sara Angelico, Dennis Mayer, Markus Gühr, Thomas J. A. Wolf, Todd J. Martínez, Henrik Koch

Abstract: The fate of thymine upon excitation by ultraviolet radiation has been the subject of intense debate over the past three decades. Today, it is widely believed that its ultrafast excited state decay stems from a radiationless transition from the bright $ππ^*$ state to a dark $nπ^*$ state. However, conflicting theoretical predictions have made the experimental data difficult to interpret. Here we sim… ▽ More The fate of thymine upon excitation by ultraviolet radiation has been the subject of intense debate over the past three decades. Today, it is widely believed that its ultrafast excited state decay stems from a radiationless transition from the bright $ππ^*$ state to a dark $nπ^*$ state. However, conflicting theoretical predictions have made the experimental data difficult to interpret. Here we simulate the ultrafast dynamics in thymine at the highest level of theory to date, performing wavepacket dynamics with a new coupled cluster method. Our simulation confirms an ultrafast $ππ^*$ to $nπ^*$ transition ($τ = 41 \pm 14$ fs). Furthermore, the predicted oxygen-edge X-ray absorption spectra agree quantitatively with the experimental results. Our simulation also predicts an as-yet uncharacterized photochemical pathway: a $πσ^*$ channel that leads to hydrogen dissociation at one of the two N-H bonds in thymine. Similar behavior has been identified in other heteroaromatic compounds, including adenine, and several authors have speculated that a similar pathway may exist in thymine. However, this was never confirmed theoretically or experimentally. This prediction calls for renewed efforts to experimentally identify or exclude the presence of this channel. △ Less

Submitted 7 March, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Comments: 42 pages, 23 figures

Coupled cluster theory for nonadiabatic dynamics: nuclear gradients and nonadiabatic couplings in similarity constrained coupled cluster theory

Authors: Eirik F. Kjønstad, Sara Angelico, Henrik Koch

Abstract: Coupled cluster theory is one of the most accurate electronic structure methods for predicting ground and excited state chemistry. However, the presence of numerical artifacts at electronic degeneracies, such as complex energies, has made it difficult to apply it in nonadiabatic dynamics simulations. While it has already been shown that such numerical artifacts can be fully removed by using simila… ▽ More Coupled cluster theory is one of the most accurate electronic structure methods for predicting ground and excited state chemistry. However, the presence of numerical artifacts at electronic degeneracies, such as complex energies, has made it difficult to apply it in nonadiabatic dynamics simulations. While it has already been shown that such numerical artifacts can be fully removed by using similarity constrained coupled cluster (SCC) theory [J. Phys. Chem. Lett. 2017, 8, 19, 4801-4807], simulating dynamics requires efficient implementations of gradients and nonadiabatic couplings. Here, we present an implementation of nuclear gradients and nonadiabatic derivative couplings at the similarity constrained coupled cluster singles and doubles (SCCSD) level of theory, thereby making possible nonadiabatic dynamics simulations using a coupled cluster theory that provides a correct description of conical intersections between excited states. We present a few numerical examples that show good agreement with literature values and discuss some limitations of the method. △ Less

Submitted 15 May, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Comments: 40 pages, 3 figures, 2 tables

Theory of magnetic properties in QED environments: application to molecular aromaticity

Authors: Alberto Barlini, Andrea Bianchi, Enrico Ronca, Henrik Koch

Abstract: In this work, we present ab initio cavity QED methods which include interactions with a static magnetic field and nuclear spin degrees of freedom using different treatments of the quantum electromagnetic field. We derive explicit expressions for QED-HF magnetizability, nuclear shielding, and spin-spin coupling tensors. We apply this theory to explore the influence of the cavity field on the magnet… ▽ More In this work, we present ab initio cavity QED methods which include interactions with a static magnetic field and nuclear spin degrees of freedom using different treatments of the quantum electromagnetic field. We derive explicit expressions for QED-HF magnetizability, nuclear shielding, and spin-spin coupling tensors. We apply this theory to explore the influence of the cavity field on the magnetizability of saturated, unsaturated, and aromatic hydrocarbons, showing the effects of different polarization orientations and coupling strengths. We also examine how the cavity affects aromaticity descriptors, such as the nucleus-independent chemical shift and magnetizability exaltation. We employ these descriptors to study the trimerization reaction of acetylene to benzene. We show how the optical cavity induces modifications in the aromatic character of the transition state leading to variations in the activation energy of the reaction. Our findings shed light on the effects induced by the cavity on magnetic properties, especially in the context of aromatic molecules, providing valuable insights into understanding the interplay between the quantum electromagnetic field and molecules. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 62 pages and 10 figures

Observation of electronic coherence created at conical intersections and its decoherence in aqueous solution

Authors: Yi-Ping Chang, Tadas Balciunas, Zhong Yin, Marin Sapunar, Bruno N. C. Tenorio, Alexander C. Paul, Shota Tsuru, Henrik Koch, Jean Pierre Wolf, Sonia Coriani, Hans Jakob Wörner

Abstract: Electronic quantum coherence is a prerequisite for charge migration in molecules and emerging molecular quantum technologies. It also underlies new forms of spectroscopy and can enhance molecular function. Whether conical intersections can create electronic coherence and for how long electronic coherence can survive in aqueous environments has remained elusive and controversial, respectively. Here… ▽ More Electronic quantum coherence is a prerequisite for charge migration in molecules and emerging molecular quantum technologies. It also underlies new forms of spectroscopy and can enhance molecular function. Whether conical intersections can create electronic coherence and for how long electronic coherence can survive in aqueous environments has remained elusive and controversial, respectively. Here, we use X-ray spectroscopy to realize a breakthrough on both of these topics. We find that electronic relaxation from the $^1$B$_{2\rm{u}}$($ππ^*$) state of pyrazine through conical intersections (CIs) induces previously unobserved electronic and vibrational quantum coherences between the $^1$B$_{3\rm{u}}$(n$π^*$) and $^1$A$_{\rm{u}}$(n$π^*$) states in the gas phase, which correspond to an electronic ring current. These coherences are entirely suppressed when pyrazine is dissolved in water. These observations, supported by the latest advances in multiconfigurational electronic-structure calculations and non-adiabatic dynamics simulations, confirm that CIs can create electronic coherences and that aqueous solvation can decohere them in less than 40 fs. This study opens the door to the investigation of the broad class of electronic coherences created during light-induced molecular dynamics and to quantify their susceptibility to aqueous solvation. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 54 pages, 28 figures

X-ray Absorption Spectra for Aqueous Ammonia and Ammonium: Quantum Mechanical versus Molecular Mechanical Embedding Schemes

Authors: Sarai Dery Folkestad, Alexander C. Paul, Regina Paul, Peter Reinholdt, Sonia Coriani, Michael Odelius, Henrik Koch

Abstract: The X-ray absorption (XA) spectra of aqueous ammonia and ammonium are computed using a combination of coupled cluster singles and doubles (CCSD) with different quantum mechanical and molecular mechanical embedding schemes. Specifically, we compare frozen Hartree--Fock (HF) density embedding, polarizable embedding (PE), and polarizable density embedding (PDE). Integrating CCSD with frozen HF densit… ▽ More The X-ray absorption (XA) spectra of aqueous ammonia and ammonium are computed using a combination of coupled cluster singles and doubles (CCSD) with different quantum mechanical and molecular mechanical embedding schemes. Specifically, we compare frozen Hartree--Fock (HF) density embedding, polarizable embedding (PE), and polarizable density embedding (PDE). Integrating CCSD with frozen HF density embedding is possible within the CC-in-HF framework, which circumvents the conventional system-size limitations of standard coupled cluster methods. We reveal similarities between PDE and frozen HF density descriptions, while PE spectra differ significantly. By including approximate triple excitations, we also investigate the effect of improving the electronic structure theory. The spectra computed using this approach show an improved intensity ratio compared to CCSD-in-HF. Charge transfer analysis of the excitations shows the local character of the pre-edge and main-edge, while the post-edge is formed by excitations delocalized over the first solvation shell and beyond. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 35 pages, 12 figures

The Hele-Shaw semi-flow

Authors: Thomas Alazard, Herbert Koch

Abstract: We prove that the Cauchy problem is well-posed in a strong sense and in a general setting. Our main result is the construction of an abstract semi-flow for the Hele-Shaw problem within general fluid domains (enabling, for instance, changes in the topology of the fluid domain) and which satisfies several properties: We provide simple comparison arguments, establish a new stability estimate and deri… ▽ More We prove that the Cauchy problem is well-posed in a strong sense and in a general setting. Our main result is the construction of an abstract semi-flow for the Hele-Shaw problem within general fluid domains (enabling, for instance, changes in the topology of the fluid domain) and which satisfies several properties: We provide simple comparison arguments, establish a new stability estimate and derive several consequences, including monotonicity and continuity results for the solutions, along with many Lyapunov functionals. We establish an eventual analytic regularity result for any arbitrary initial data. We also study numerous qualitative properties, including global regularity for initial data in sub-critical Sobolev spaces, well-posedness in a strong sense for initial data with barely a modulus of continuity, as well as waiting-time phenomena for Lipschitz solutions, in any dimension. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 62 pages

MSC Class: 76D27; 35R37

Collective Strong Coupling Modifies Aggregation and Solvation

Authors: Matteo Castagnola, Tor S. Haugland, Enrico Ronca, Henrik Koch, Christian Schäfer

Abstract: Intermolecular interactions are pivotal for aggregation, solvation, and crystallization. We demonstrate that the collective strong coupling of several molecules to a single optical mode results in notable changes in the molecular excitations around an impurity, e.g., in the first aggregation or solvation shell. A competition between short-range Coulombic and long-range photonic correlation inverts… ▽ More Intermolecular interactions are pivotal for aggregation, solvation, and crystallization. We demonstrate that the collective strong coupling of several molecules to a single optical mode results in notable changes in the molecular excitations around an impurity, e.g., in the first aggregation or solvation shell. A competition between short-range Coulombic and long-range photonic correlation inverts the local transition density in a polaritonic state, suggesting notable changes in the polarizability of the solvation shell. Our results provide an alternative perspective on recent work in polaritonic chemistry and pave the way for the rigorous treatment of cooperative effects in aggregation, solvation, and crystallization. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Journal ref: J. Phys. Chem. Lett. 2024, 15, 1428-1434

Model-independent extraction of form factors and $|V_{cb}|$ in $\overline{B} \rightarrow D \ell^- \overlineν_\ell$ with hadronic tagging at BaBar

Authors: BaBar Collaboration, J. P. Lees, V. Poireau, V. Tisserand, E. Grauges, A. Palano, G. Eigen, D. N. Brown, Yu. G. Kolomensky, M. Fritsch, H. Koch, R. Cheaib, C. Hearty, T. S. Mattison, J. A. McKenna, R. Y. So, V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, E. A. Kozyrev, E. A. Kravchenko, S. I. Serednyakov, Yu. I. Skovpen, E. P. Solodov, K. Yu. Todyshev , et al. (186 additional authors not shown)

Abstract: Using the entire BaBar $Υ(4S)$ data set, the first two-dimensional unbinned angular analysis of the semileptonic decay $\overline{B} \rightarrow D \ell^- \overlineν_\ell$ is performed, employing hadronic reconstruction of the tag-side $B$ meson from $Υ(4S)\to B\overline{B}$. Here, $\ell$ denotes the light charged leptons $e$ and $μ$. A novel data-driven signal-background separation procedure with… ▽ More Using the entire BaBar $Υ(4S)$ data set, the first two-dimensional unbinned angular analysis of the semileptonic decay $\overline{B} \rightarrow D \ell^- \overlineν_\ell$ is performed, employing hadronic reconstruction of the tag-side $B$ meson from $Υ(4S)\to B\overline{B}$. Here, $\ell$ denotes the light charged leptons $e$ and $μ$. A novel data-driven signal-background separation procedure with minimal dependence on simulation is developed. This procedure preserves all multi-dimensional correlations present in the data. The expected $\sin^2θ_\ell$ dependence of the differential decay rate in the Standard Model is demonstrated, where $θ_\ell$ is the lepton helicity angle. Including input from the latest lattice QCD calculations and previously available experimental data, the underlying form factors are extracted using both model-independent (BGL) and dependent (CLN) methods. Comparisons with lattice calculations show flavor SU(3) symmetry to be a good approximation in the $B_{(s)}\to D_{(s)}$ sector. Using the BGL results, the CKM matrix element $|V_{cb}|=(41.09\pm 1.16)\times 10^{-3}$ and the Standard Model prediction of the lepton-flavor universality violation variable $\mathcal{R}(D)=0.300\pm 0.004$, are extracted. The value of $|V_{cb}|$ from $\overline{B} \rightarrow D \ell^- \overlineν_\ell$ tends to be higher than that extracted using $\overline{B} \rightarrow D \ell^- \overlineν_\ell$. The Standard Model $\mathcal{R}(D)$ calculation is at a $1.97σ$ tension with the latest HFLAV experimental average. △ Less

Submitted 25 November, 2023; originally announced November 2023.

PlotMap: Automated Layout Design for Building Game Worlds

Authors: Yi Wang, Jieliang Luo, Adam Gaier, Evan Atherton, Hilmar Koch

Abstract: World-building, the process of developing both the narrative and physical world of a game, plays a vital role in the game's experience. Critically acclaimed independent and AAA video games are praised for strong world building, with game maps that masterfully intertwine with and elevate the narrative, captivating players and leaving a lasting impression. However, designing game maps that support a… ▽ More World-building, the process of developing both the narrative and physical world of a game, plays a vital role in the game's experience. Critically acclaimed independent and AAA video games are praised for strong world building, with game maps that masterfully intertwine with and elevate the narrative, captivating players and leaving a lasting impression. However, designing game maps that support a desired narrative is challenging, as it requires satisfying complex constraints from various considerations. Most existing map generation methods focus on considerations about gameplay mechanics or map topography, while the need to support the story is typically neglected. As a result, extensive manual adjustment is still required to design a game world that facilitates particular stories. In this work, we approach this problem by introducing an extra layer of plot facility layout design that is independent of the underlying map generation method in a world-building pipeline. Concretely, we present a system that leverages Reinforcement Learning (RL) to automatically assign concrete locations on a game map to abstract locations mentioned in a given story (plot facilities), following spatial constraints derived from the story. A decision-making agent moves the plot facilities around, considering their relationship to the map and each other, to locations on the map that best satisfy the constraints of the story. Our system considers input from multiple modalities: map images as pixels, facility locations as real values, and story constraints expressed in natural language. We develop a method of generating datasets of facility layout tasks, create an RL environment to train and evaluate RL models, and further analyze the behaviors of the agents through a group of comprehensive experiments and ablation studies, aiming to provide insights for RL-based plot facility layout design. △ Less

Submitted 15 June, 2024; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: In 2019 IEEE Conference on Games (CoG)

Wellposedness for the KdV hierarchy

Authors: Friedrich Klaus, Herbert Koch, Baoping Liu

Abstract: We prove a version of wellposedness for all equations of the KdV hierarchy in $H^{-1}$. Ingredients are 1) The Miura map which allows to define the Gardner hierarchy through the generating function of the energies so that the $N$th Gardner equation is equivalent to the $N$th KdV equation. 2) A rigorous relation between the generating functions of the energies and the KdV resp. Gardner Hamilton… ▽ More We prove a version of wellposedness for all equations of the KdV hierarchy in $H^{-1}$. Ingredients are 1) The Miura map which allows to define the Gardner hierarchy through the generating function of the energies so that the $N$th Gardner equation is equivalent to the $N$th KdV equation. 2) A rigorous relation between the generating functions of the energies and the KdV resp. Gardner Hamiltonians. 3) Kato smoothing estimates for weak solutions and approximate flows. △ Less

Submitted 22 September, 2023; originally announced September 2023.

Comments: 139 pages

MSC Class: 37K10; 35Q53

Coupled cluster cavity Born-Oppenheimer approximation for electronic strong coupling

Authors: Sara Angelico, Tor S. Haugland, Enrico Ronca, Henrik Koch

Abstract: Chemical and photochemical reactivity, as well as supramolecular organization and several other molecular properties, can be modified by strong interactions between light and matter. Theoretical studies of these phenomena require the separation of the Schrödinger equation into different degrees of freedom as in the Born-Oppenheimer approximation. In this paper, we analyze the electron-photon Hamil… ▽ More Chemical and photochemical reactivity, as well as supramolecular organization and several other molecular properties, can be modified by strong interactions between light and matter. Theoretical studies of these phenomena require the separation of the Schrödinger equation into different degrees of freedom as in the Born-Oppenheimer approximation. In this paper, we analyze the electron-photon Hamiltonian within the cavity Born-Oppenheimer approximation (CBOA), where the electronic problem is solved for fixed nuclear positions and photonic parameters. Specifically, we focus on intermolecular interactions in representative dimer complexes. The CBOA potential energy surfaces are compared with those obtained using a polaritonic approach, where the photonic and electronic degrees of freedom are treated at the same level. This allows us to assess the role of electron-photon correlation and the accuracy of CBOA. △ Less

Submitted 2 November, 2023; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 9 pages, 8 figures

Understanding X-ray absorption in liquid water: triple excitations in multilevel coupled cluster theory

Authors: Sarai Dery Folkestad, Alexander C. Paul, Regina Paul, Sonia Coriani, Michael Odelius, Marcella Iannuzzi, Henrik Koch

Abstract: We present the first successful application of the coupled cluster approach to simulate the X-ray absorption (XA) spectrum of liquid water. The system size limitations of standard coupled cluster theory are overcome by employing a newly developed coupled cluster method for large molecular systems. This method combines coupled cluster singles, doubles, and perturbative triples in a multilevel frame… ▽ More We present the first successful application of the coupled cluster approach to simulate the X-ray absorption (XA) spectrum of liquid water. The system size limitations of standard coupled cluster theory are overcome by employing a newly developed coupled cluster method for large molecular systems. This method combines coupled cluster singles, doubles, and perturbative triples in a multilevel framework (MLCC3-in-HF) and is able to describe the delicate nature of intermolecular interactions in liquid water. Using molecular geometries from state-of-the-art path-integral molecular dynamics, we obtain excellent agreement with experimental spectra. Additionally, we show that an accurate description of the electronic structure within the first solvation shell is sufficient to model the XA spectrum of liquid water. Furthermore, we present a rigorous charge transfer analysis with unprecedented reliability, achieved through MLCC3-in-HF. This analysis aligns with previous studies regarding the character of the prominent features of the spectrum. △ Less

Submitted 19 December, 2023; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 32 pages, 4+12 figures

Asymptotic stability of the sine-Gordon kinks under perturbations in weighted Sobolev norms

Authors: Herbert Koch, Dongxiao Yu

Abstract: We study the asymptotic stability of the sine-Gordon kinks under small perturbations in weighted Sobolev norms. Our main tool is the Bäcklund transform which reduces the study of the asymptotic stability of the kinks to the study of the asymptotic decay of solutions near zero. Our results consist of two parts. First, we present a different proof of the local asymptotic stability result in arXiv:20… ▽ More We study the asymptotic stability of the sine-Gordon kinks under small perturbations in weighted Sobolev norms. Our main tool is the Bäcklund transform which reduces the study of the asymptotic stability of the kinks to the study of the asymptotic decay of solutions near zero. Our results consist of two parts. First, we present a different proof of the local asymptotic stability result in arXiv:2009.04260. In its proof, we apply a result obtained by the inverse scattering method on the local decay of the solutions with sufficiently small and localized initial data. Moreover, we prove an $L^\infty$-type asymptotic stability result which is similar to that in arXiv:2106.09605; the main difference is that we remove the assumptions on the spatial symmetry of the perturbations. In its proof, we apply a result obtained by the method of testing by wave packets on the pointwise decay of the solutions with small and localized data. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: 54 pages

Strong coupling to circularly polarized photons: towards field-induced enantioselectivity

Authors: Rosario R. Riso, Enrico Ronca, Henrik Koch

Abstract: The development of new methodologies for the selective synthesis of individual enantiomers is still one of the major challenges in synthetic chemistry. Many biomolecules, and therefore many pharmaceutical compounds, are indeed chiral. While the use of chiral reactants or catalysts has led to significant progress in the field of asymmetric synthesis, a systematic approach applicable to general reac… ▽ More The development of new methodologies for the selective synthesis of individual enantiomers is still one of the major challenges in synthetic chemistry. Many biomolecules, and therefore many pharmaceutical compounds, are indeed chiral. While the use of chiral reactants or catalysts has led to significant progress in the field of asymmetric synthesis, a systematic approach applicable to general reactions has still not been proposed. In this work, we show that strong coupling to circularly polarized fields represents a promising alternative for reaching highly selective asymmetric synthesis in a non-invasive fashion. We demonstrate that the field induces stereoselectivity in the early stages of the chemical reaction, by selecting an energetically preferred direction of approach for the reagents. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Comments: 15 pages, 5 figures

Understanding the polaritonic ground state in cavity quantum electrodynamics

Authors: Tor S. Haugland, John P. Philbin, Tushar K. Ghosh, Ming Chen, Henrik Koch, Prineha Narang

Abstract: Molecular polaritons arise when molecules interact so strongly with light that they become entangled with each other. This light-matter hybridization alters the chemical and physical properties of the molecular system and allows chemical reactions to be controlled without the use of external fields. We investigate the impact of strong light-matter coupling on the electronic structure using perturb… ▽ More Molecular polaritons arise when molecules interact so strongly with light that they become entangled with each other. This light-matter hybridization alters the chemical and physical properties of the molecular system and allows chemical reactions to be controlled without the use of external fields. We investigate the impact of strong light-matter coupling on the electronic structure using perturbative approaches and demonstrate that Rayleigh-Schrödinger perturbation theory can reproduce the ground state energies in optical cavities to comparable accuracy as ab initio cavity quantum electrodynamics methodologies for currently relevant coupling strengths. The method is effective in both low and high cavity frequency regimes and straightforward to implement via response functions. Furthermore, we establish simple relations between cavity-induced intermolecular forces and van der Waals forces. These findings provide valuable insight into the manipulation of ground-state polaritonic energy landscapes, shedding light on the systems and conditions in which modifications can be achieved. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 11 pages, 3 Figures

Coupled cluster simulation of impulsive stimulated X-ray Raman scattering

Authors: Alice Balbi, Andreas S. Skeidsvoll, Henrik Koch

Abstract: Time-dependent equation-of-motion coupled cluster (TD-EOM-CC) is used to simulate impulsive stimulated x-ray Raman scattering (ISXRS) of ultrashort laser pulses by neon, carbon monoxide, pyrrole, and p-aminophenol. The TD-EOM-CC equations are expressed in the basis of field-free EOM-CC states, where the calculation of the core-excited states is simplified through the use of the core-valence separa… ▽ More Time-dependent equation-of-motion coupled cluster (TD-EOM-CC) is used to simulate impulsive stimulated x-ray Raman scattering (ISXRS) of ultrashort laser pulses by neon, carbon monoxide, pyrrole, and p-aminophenol. The TD-EOM-CC equations are expressed in the basis of field-free EOM-CC states, where the calculation of the core-excited states is simplified through the use of the core-valence separation (CVS) approximation. The transfer of electronic population from the ground state to the core- and valence-excited states is calculated for different numbers of included core- and valence-excited states, as well as for electric field pulses with different polarizations and carrier frequencies. The results indicate that Gaussian pulses can transfer significant electronic populations to the valence states through the Raman process. The sensitivity of this population transfer to the model parameters is analyzed. The time-dependent electronic density for p-aminophenol is also showcased, supporting the interpretation that ISXRS involves localized core excitations and can be used to rapidly generate valence wavepackets. △ Less

Submitted 31 May, 2023; originally announced May 2023.

Comments: 10 pages, 5 figures

Polaritonic response theory for exact and approximate wave functions

Authors: Matteo Castagnola, Rosario R. Riso, Alberto Barlini, Enrico Ronca, Henrik Koch

Abstract: Polaritonic chemistry is an interdisciplinary emerging field that presents several challenges and opportunities in chemistry, physics, and engineering. A systematic review of polaritonic response theory is presented, following a chemical perspective based on molecular response theory. We provide the reader with a general strategy for developing response theory for ab initio cavity quantum electrod… ▽ More Polaritonic chemistry is an interdisciplinary emerging field that presents several challenges and opportunities in chemistry, physics, and engineering. A systematic review of polaritonic response theory is presented, following a chemical perspective based on molecular response theory. We provide the reader with a general strategy for developing response theory for ab initio cavity quantum electrodynamics (QED) methods and critically emphasize details that still need clarification and require cooperation between the physical and chemistry communities. We show that several well-established results can be applied to strong coupling light-matter systems, leading to novel perspectives on the computation of matter and photonic properties. The application of the Pauli-Fierz Hamiltonian to polaritons is discussed, focusing on the effects of describing operators in different mathematical representations. We thoroughly examine the most common approximations employed in ab initio QED, such as the dipole approximation. We introduce the polaritonic response equations for recently developed ab initio QED Hartree-Fock and QED coupled cluster methods. The discussion focuses on the similarities and differences from standard quantum chemistry methods, providing practical equations for computing the polaritonic properties. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 37 pages, 3 figures

Effective single mode methodology for strongly coupled multimode molecular-plasmon nanosystems

Authors: Marco Romanelli, Rosario R. Riso, Tor S. Haugland, Enrico Ronca, Stefano Corni, Henrik Koch

Abstract: Strong coupling between molecules and quantized fields has emerged as an effective methodology to engineer molecular properties. New hybrid states are formed when molecules interact with quantized fields. Since the properties of these states can be modulated by fine-tuning the field features, an exciting and new side of chemistry can be explored. In particular, significant modifications of the mol… ▽ More Strong coupling between molecules and quantized fields has emerged as an effective methodology to engineer molecular properties. New hybrid states are formed when molecules interact with quantized fields. Since the properties of these states can be modulated by fine-tuning the field features, an exciting and new side of chemistry can be explored. In particular, significant modifications of the molecular properties can be achieved in plasmonic nanocavities, where the field quantization volume is reduced to sub-nanometric volumes. Intriguing applications of nanoplasmonics include the possibility of coupling the plasmons with a single molecule, instrumental for sensing, high-resolution spectroscopy, and single-molecule imaging. In this work, we focus on phenomena where the simultaneous effects of multiple plasmonic modes are critical. We propose a theoretical methodology to account for many plasmonic modes simultaneously while retaining computational feasibility. Our approach is conceptually simple and allows us to accurately account for the multimode effects and rationalize the nature of the interaction between multiple plasmonic excitations and molecules. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Comments: 27 pages, 6 figures

Search for $B$ Mesogenesis at BABAR

Authors: BABAR Collaboration, J. P. Lees, V. Poireau, V. Tisserand, E. Grauges, A. Palano, G. Eigen, D. N. Brown, Yu. G. Kolomensky, M. Fritsch, H. Koch, R. Cheaib, C. Hearty, T. S. Mattison, J. A. McKenna, R. Y. So, V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev, E. A. Kozyrev, E. A. Kravchenko, A. P. Onuchin, S. I. Serednyakov, Yu. I. Skovpen , et al. (218 additional authors not shown)

Abstract: A new mechanism has been proposed to simultaneously explain the presence of dark matter and the matter-antimatter asymmetry in the universe. This scenario predicts exotic $B$ meson decays into a baryon and a dark sector anti-baryon ($ψ_D$) with branching fractions accessible at $B$ factories. We present a search for $B \rightarrow Λψ_D$ decays using data collected by the $BABAR$ experiment at SLAC… ▽ More A new mechanism has been proposed to simultaneously explain the presence of dark matter and the matter-antimatter asymmetry in the universe. This scenario predicts exotic $B$ meson decays into a baryon and a dark sector anti-baryon ($ψ_D$) with branching fractions accessible at $B$ factories. We present a search for $B \rightarrow Λψ_D$ decays using data collected by the $BABAR$ experiment at SLAC. This reaction is identified by fully reconstructing the accompanying $B$ meson and requiring the presence of a single $Λ$ baryon in the remaining particles. No significant signal is observed, and bounds on the $B \rightarrow Λψ_D$ branching fraction are derived in the range $0.13 - 5.2\times 10^{-5}$ for $1.0 < m_{ψ_D} < 4.2$ GeV/$c^{2}$. These results set strong constraints on the parameter space allowed by the theory. △ Less

Submitted 31 January, 2023; originally announced February 2023.

Journal ref: PHYS. REV. D 107, 092001 (2023)

Non-adiabatic derivative coupling elements for the coupled cluster singles and doubles model

Authors: Eirik F. Kjønstad, Henrik Koch

Abstract: We present an efficient implementation of analytical non-adiabatic derivative coupling elements for the coupled cluster singles and doubles model. The derivative coupling elements are evaluated in a biorthonormal formulation in which the nuclear derivative acts on the right electronic state, where this state is biorthonormal with respect to the set of left states. This stands in contrast to earlie… ▽ More We present an efficient implementation of analytical non-adiabatic derivative coupling elements for the coupled cluster singles and doubles model. The derivative coupling elements are evaluated in a biorthonormal formulation in which the nuclear derivative acts on the right electronic state, where this state is biorthonormal with respect to the set of left states. This stands in contrast to earlier implementations based on normalized states and a gradient formula for the derivative coupling. As an illustration of the implementation, we determine a minimum energy conical intersection between the nπ* and ππ* states in the nucleobase thymine. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 8 pages, 4 figures

Entanglement Coupled Cluster Theory: Exact Spin-Adaptation

Authors: Sarai Dery Folkestad, Bendik Støa Sannes, Henrik Koch

Abstract: We present a novel approach to spin-adapted coupled cluster theory. This approach is based on the entanglement of an open-shell molecule with electrons in a non-interacting bath; together they form a closed-shell state. For the total system, electron correlation can be added with the standard spin-adapted closed-shell coupled cluster formalism. A projection operator, which enforces conditions on t… ▽ More We present a novel approach to spin-adapted coupled cluster theory. This approach is based on the entanglement of an open-shell molecule with electrons in a non-interacting bath; together they form a closed-shell state. For the total system, electron correlation can be added with the standard spin-adapted closed-shell coupled cluster formalism. A projection operator, which enforces conditions on the electrons in the bath, is used to obtain the desired state of the molecular system. This entanglement coupled cluster theory is detailed and proof-of-concept calculations for doublet states are reported. The approach is extendable to open-shell systems with other values of the total spin. △ Less

Submitted 19 May, 2023; v1 submitted 28 January, 2023; originally announced January 2023.

Comments: 17 pages, 3 figures and 7 tables

Comparing real-time coupled cluster methods through simulation of collective Rabi oscillations

Authors: Andreas S. Skeidsvoll, Henrik Koch

Abstract: The time-dependent equation-of-motion coupled cluster (TD-EOM-CC) and time-dependent coupled cluster (TDCC) methods are compared by simulating Rabi oscillations for different numbers of non-interacting atoms in a classical electromagnetic field. While the TD-EOM-CC simulations are numerically stable, the oscillating time-dependent energy scales unreasonably with the number of subsystems resonant w… ▽ More The time-dependent equation-of-motion coupled cluster (TD-EOM-CC) and time-dependent coupled cluster (TDCC) methods are compared by simulating Rabi oscillations for different numbers of non-interacting atoms in a classical electromagnetic field. While the TD-EOM-CC simulations are numerically stable, the oscillating time-dependent energy scales unreasonably with the number of subsystems resonant with the field. The TDCC simulations give the correct scaling of the time-dependent energy in the initial stages of the Rabi cycle, but the numerical solution breaks down when the multi-atom system approaches complete population inversion. We present a general theoretical framework in which the two methods can be described, where the cluster amplitude time derivatives are taken as auxiliary conditions, leading to a shifted time-dependent Hamiltonian matrix. In this framework, TDCC has a shifted Hamiltonian with a block upper triangular structure, explaining the correct scaling properties of the method. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Comments: 9 pages and 4 figures

Higher Regularity of the Free Boundary in a Semilinear System

Authors: Morteza Fotouhi, Herbert Koch

Abstract: In this paper we are concerned with higher regularity properties of the elliptic system \[ Δ\mathbf{u}= |\mathbf{u}|^{q-1}\mathbf{u}χ_{\{|\mathbf{u}|>0\}},\qquad\mathbf{u}=(u^1,\dots,u^m) \] for $0\leq q<1$. We show analyticity of the regular part of the free boundary $\partial\{|\mathbf{u}|>0\}$, analyticity of $|\mathbf{u}|^{\frac{1-q}2} $ and $ \frac{\mathbf{u}}{|\mathbf{u}|}$ up to the regular… ▽ More In this paper we are concerned with higher regularity properties of the elliptic system \[ Δ\mathbf{u}= |\mathbf{u}|^{q-1}\mathbf{u}χ_{\{|\mathbf{u}|>0\}},\qquad\mathbf{u}=(u^1,\dots,u^m) \] for $0\leq q<1$. We show analyticity of the regular part of the free boundary $\partial\{|\mathbf{u}|>0\}$, analyticity of $|\mathbf{u}|^{\frac{1-q}2} $ and $ \frac{\mathbf{u}}{|\mathbf{u}|}$ up to the regular part of the free boundary. Applying a variant of the partial hodograph-Legendre transformation and the implicit function theorem, we arrive at a degenerate equation, which introduces substantial challenges to be dealt with. Along the lines of our study, we also establish a Cauchy-Kowalevski type statement to show the local existence of solution when the free boundary and the restriction of $ \frac{\mathbf{u}}{|\mathbf{u}|} $ from both sides to the free boundary are given as analytic data. △ Less

Submitted 1 May, 2023; v1 submitted 26 December, 2022; originally announced December 2022.

Linear Response Properties of Solvated Systems: A Computational Study

Authors: Linda Goletto, Sara Gómez, Josefine H. Andersen, Henrik Koch, Tommaso Giovannini

Abstract: We present a computational study of static and dynamic linear polarizabilities in solution. We use different theoretical approaches to describe solvent effects, ranging from quantum mechanics/molecular mechanics (QM/MM) to quantum embedding approaches. In particular, we consider non-polarizable and polarizable QM/MM methods, the latter based on the fluctuating charge (FQ) force field. In addition,… ▽ More We present a computational study of static and dynamic linear polarizabilities in solution. We use different theoretical approaches to describe solvent effects, ranging from quantum mechanics/molecular mechanics (QM/MM) to quantum embedding approaches. In particular, we consider non-polarizable and polarizable QM/MM methods, the latter based on the fluctuating charge (FQ) force field. In addition, we use a quantum embedding method defined in the context of multilevel Hartree-Fock (MLHF), where the system is divided into active and inactive regions, and combine it with a third layer described by means of the FQ model. The multiscale approaches are then used as reference wave functions for equation-of-motion coupled cluster (EOM-CC) response properties, allowing for the account of electron correlation. The developed models are applied to the calculation of linear response properties of two organic moieties -- namely, para-nitroaniline and benzonitrile -- in non-aqueous solvents -- 1,4-dioxane, acetonitrile, and tetrahydrofuran. The computed polarizabilities are then discussed in terms of the physico-chemical solute-solvent interactions described by each method (electrostatic, polarization and Pauli repulsion), and finally compared with the available experimental references. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Comments: 13 pages, 9 figures

Molecular van der Waals fluids in cavity quantum electrodynamics

Authors: John P. Philbin, Tor S. Haugland, Tushar K. Ghosh, Enrico Ronca, Ming Chen, Prineha Narang, Henrik Koch

Abstract: Intermolecular van der Waals interactions are central to chemical and physical phenomena ranging from biomolecule binding to soft-matter phase transitions. However, there are currently very limited approaches to manipulate van der Waals interactions. In this work, we demonstrate that strong light-matter coupling can be used to tune van der Waals interactions, and, thus, control the thermodynamic p… ▽ More Intermolecular van der Waals interactions are central to chemical and physical phenomena ranging from biomolecule binding to soft-matter phase transitions. However, there are currently very limited approaches to manipulate van der Waals interactions. In this work, we demonstrate that strong light-matter coupling can be used to tune van der Waals interactions, and, thus, control the thermodynamic properties of many-molecule systems. Our analyses reveal orientation dependent single molecule energies and interaction energies for van der Waals molecules (for example, H$_{2}$). For example, we find intermolecular interactions that depend on the distance between the molecules $R$ as $R^{-3}$ and $R^{0}$. Moreover, we employ non-perturbative \textit{ab initio} cavity quantum electrodynamics calculations to develop machine learning-based interaction potentials for molecules inside optical cavities. By simulating systems ranging from $12$ H$_2$ to $144$ H$_2$ molecules, we demonstrate that strong light-matter coupling can tune the structural and thermodynamic properties of molecular fluids. In particular, we observe varying degrees of orientational order as a consequence of cavity-modified interactions, and we explain how quantum nuclear effects, light-matter coupling strengths, number of cavity modes, molecular anisotropies, and system size all impact the extent of orientational order. These simulations and analyses demonstrate both local and collective effects induced by strong light-matter coupling and open new paths for controlling the properties of molecular clusters. △ Less

Submitted 30 June, 2023; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 7 pages and 3 figures in main text - 16 pages and 13 figures supplemental material

Strong coupling in chiral cavities: nonperturbative framework for enantiomer discrimination

Authors: Rosario R. Riso, Laura Grazioli, Enrico Ronca, Tommaso Giovannini, Henrik Koch

Abstract: The development of efficient techniques to distinguish mirror images of chiral molecules (enantiomers) is very important in both chemistry and physics. Enantiomers share most molecular properties except, for instance, the absorption of circularly polarized light. Enantiomer purification is therefore a challenging task that requires specialized equipment. Strong coupling between quantized fields an… ▽ More The development of efficient techniques to distinguish mirror images of chiral molecules (enantiomers) is very important in both chemistry and physics. Enantiomers share most molecular properties except, for instance, the absorption of circularly polarized light. Enantiomer purification is therefore a challenging task that requires specialized equipment. Strong coupling between quantized fields and matter (e.g. in optical cavities) is a promising technique to modify molecular processes in a non-invasive way. The modulation of molecular properties is achieved by changing the field characteristics. In this work, we investigate whether strong coupling to circularly polarized electromagnetic fields is a viable way to discriminate chiral molecules. To this end, we develop a nonperturbative framework to calculate the behavior of molecules in chiral cavities. We show that in this setting the enantiomers have different energies -- that is, one is more stable than the other. The field-induced energy differences are also shown to give rise to enantiospecific signatures in rotational spectra. △ Less

Submitted 6 February, 2023; v1 submitted 5 September, 2022; originally announced September 2022.

Comments: 18 pages and 14 figures

Conserved energies for the one dimensional Gross-Pitaevskii equation: low regularity case

Authors: Herbert Koch, Xian Liao

Abstract: We construct a family of conserved energies for the one dimensional Gross-Pitaevskii equation, but in the low regularity case (in \cite{KL} we have constructed conserved energies in the high regularity situation). This can be done thanks to regularization procedures and a study of the topological structure of the finite-energy space. The asymptotic (regularised conserved) phase change on the real… ▽ More We construct a family of conserved energies for the one dimensional Gross-Pitaevskii equation, but in the low regularity case (in \cite{KL} we have constructed conserved energies in the high regularity situation). This can be done thanks to regularization procedures and a study of the topological structure of the finite-energy space. The asymptotic (regularised conserved) phase change on the real line with values in $ \R/2π\Z$ is studied. We also construct a conserved quantity, the renormalized momentum $H_1$ (see Theorem \ref{thm:E1}), on the universal covering space of the finite-energy space. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 45 pages

On the characteristic features of ionization in QED environments

Authors: Rosario R. Riso, Tor S. Haugland, Enrico Ronca, Henrik Koch

Abstract: 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 o… ▽ More 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. △ Less

Submitted 6 May, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: 8 pages, 4 figures

Universal supercritical behavior for some skew-product maps

Authors: Hans Koch

Abstract: We consider skew-product maps over circle rotations $x\mapsto x+α$ (mod 1) with factors that take values in SL(2,R) In numerical experiments with $α$ the inverse golden mean, Fibonacci iterates of almost Mathieu maps with rotation number $1/4$ and positive Lyapunov exponent exhibit asymptotic scaling behavior. We prove the existence of such asymptotic scaling for "periodic" rotation numbers and fo… ▽ More We consider skew-product maps over circle rotations $x\mapsto x+α$ (mod 1) with factors that take values in SL(2,R) In numerical experiments with $α$ the inverse golden mean, Fibonacci iterates of almost Mathieu maps with rotation number $1/4$ and positive Lyapunov exponent exhibit asymptotic scaling behavior. We prove the existence of such asymptotic scaling for "periodic" rotation numbers and for large Lyapunov exponent. The phenomenon is universal, in the sense that it holds for open sets of maps, with the scaling limit being independent of the maps. The set of maps with a given periodic rotation number is a real analytic codimension $1$ manifold in a suitable space of maps. △ Less

Submitted 3 March, 2022; originally announced March 2022.

MSC Class: 37E20 37E45 47B93 37A05

Oscillator strengths in the framework of equation of motion multilevel CC3

Authors: Alexander C. Paul, Sarai D. Folkestad, Rolf H. Myhre, Henrik Koch

Abstract: 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. Asym… ▽ More 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. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 38 page and 10 pages SI and 6 figures

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

Authors: Anna Kristina Schnack-Petersen, Henrik Koch, Sonia Coriani, Eirik F. Kjønstad

Abstract: 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… ▽ More 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. △ Less

Submitted 4 February, 2022; originally announced February 2022.

Comments: 12 pages, 1 figure

Extendible Functions and Local Root Numbers Remarks on a paper of R.P. Langlands

Authors: Helmut Koch, Ernst-Wilhelm Zink

Abstract: This paper refers to Langlands' big set of notes [L] devoted to the question if the (normalized) local Hecke-Tate root number $Δ=Δ(E,χ)$, where $E$ is a finite separable extension of a fixed non-archimedean local field $F$, and $χ$ a quasicharacter of $E^\times$, can be appropriately extended to a local $\varepsilon$-factor $\varepsilon_Δ=\varepsilon_Δ(E,ρ)$ for all virtual representations $ρ$ of… ▽ More This paper refers to Langlands' big set of notes [L] devoted to the question if the (normalized) local Hecke-Tate root number $Δ=Δ(E,χ)$, where $E$ is a finite separable extension of a fixed non-archimedean local field $F$, and $χ$ a quasicharacter of $E^\times$, can be appropriately extended to a local $\varepsilon$-factor $\varepsilon_Δ=\varepsilon_Δ(E,ρ)$ for all virtual representations $ρ$ of the corresponding Weil group $W_E.$ Whereas Deligne [D] has given a relatively short proof by using the global Artin-Weil L-functions, the proof of Langlands is purely local and splits into two parts: the {\bf algebraic part} to find a minimal set of relations for the functions $Δ$, such that the existence (and uniqueness) of $\varepsilon_Δ$ will follow from these relations; and the more extensive {\bf arithmetic part} to give a direct proof that all these relations are actually fulfilled. Our aim is to cover the algebraic part of Langlands' notes which can be done completely in the framework of representations of solvable profinite groups, where two modifications of Brauer's theorem play a prominent role. Introduction / 1. The notion of extendible functions / 2. The kernel of the Brauer map and its generating relations / 3. A criterion for the extendibility of functions / 4. Recovering Theorem 3.1 for the case of local root numbers / Appendices A1. Proving Brauer 3 and Brauer 4 / A2. On type-III-groups. △ Less

Submitted 25 June, 2022; v1 submitted 29 December, 2021; originally announced December 2021.

Comments: 50 pages

MSC Class: 11S37; 11S20; 11S40; 20C15

Molecular orbital theory in cavity QED environments

Authors: Rosario R. Riso, Tor S. Haugland, Enrico Ronca, Henrik Koch

Abstract: 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… ▽ More 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. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Comments: 23 pages, 7 figures

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

Authors: Andreas S. Skeidsvoll, Torsha Moitra, Alice Balbi, Alexander C. Paul, Sonia Coriani, Henrik Koch

Abstract: 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 sp… ▽ More 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. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 15 pages and 5 figures

Dirichlet problem for weakly harmonic maps with rough data

Authors: Gael Diebou Yomgne, Herbert Koch

Abstract: Weakly harmonic maps from a domain $Ω$ (the upper half-space $\Rd$ or a bounded $C^{1,α}$ domain, $α\in (0,1]$) into a smooth closed manifold are studied. Prescribing small Dirichlet data in either of the classes $L^{\infty}(\partialΩ)$ or $BMO(\partialΩ)$, we establish solvability of the resulting boundary value problems by means of a nonvariational method. As a by-product, solutions are shown to… ▽ More Weakly harmonic maps from a domain $Ω$ (the upper half-space $\Rd$ or a bounded $C^{1,α}$ domain, $α\in (0,1]$) into a smooth closed manifold are studied. Prescribing small Dirichlet data in either of the classes $L^{\infty}(\partialΩ)$ or $BMO(\partialΩ)$, we establish solvability of the resulting boundary value problems by means of a nonvariational method. As a by-product, solutions are shown to be locally smooth, $C^{\infty}_{loc}$. Moreover, we show that boundary data can be chosen large in the underlying topologies if $Ω$ is smooth and bounded by perturbing strictly stable smooth harmonic maps. △ Less

Submitted 8 October, 2021; originally announced October 2021.

MSC Class: 35B20; 35G60; 42B35; 58E20

Strong coupling between localized surface plasmons and molecules by coupled cluster theory

Authors: Jacopo Fregoni, Tor S. Haugland, Silvio Pipolo, Tommaso Giovannini, Henrik Koch, Stefano Corni

Abstract: Plasmonic nanocavities enable the confinement of molecules and electromagnetic fields within nano-metric volumes. As a consequence, the molecules experience a remarkably strong interaction with the electromagnetic field, to such an extent that the quantum states of the system become hybrids between light and matter: polaritons. Here we present a non-perturbative method to simulate the emerging pro… ▽ More Plasmonic nanocavities enable the confinement of molecules and electromagnetic fields within nano-metric volumes. As a consequence, the molecules experience a remarkably strong interaction with the electromagnetic field, to such an extent that the quantum states of the system become hybrids between light and matter: polaritons. Here we present a non-perturbative method to simulate the emerging properties of such polaritons: it combines a high-level quantum chemical description of the molecule with a quantized description of the localized surface plasmons in the nanocavity. We apply the method to molecules of realistic complexity in a typical plasmonic nanocavity, featuring also a subnanometric asperity (picocavity). Our results disclose the effects of the mutual polarization and correlation of plasmons and molecular excitations, disregarded so far. They also quantify to what extent the molecular charge density can be manipulated by nanocavities, and stand as benchmarks to guide the development of methods for molecular polaritonics. △ Less

Submitted 12 July, 2021; originally announced July 2021.

Comments: 13 Pages Main text, 4 Figures and 11 pages Supporting Information

Asymptotic scaling and universality for skew products with factors in SL(2,R)

Authors: Hans Koch

Abstract: We consider skew-product maps over circle rotations $x\mapsto x+α$ (mod 1) with factors that take values in SL(2,R). This includes maps of almost Mathieu type. In numerical experiments, with $α$ the inverse golden mean, Fibonacci iterates of maps from the almost Mathieu family exhibit asymptotic scaling behavior that is reminiscent of critical phase transitions. In a restricted setup that is chara… ▽ More We consider skew-product maps over circle rotations $x\mapsto x+α$ (mod 1) with factors that take values in SL(2,R). This includes maps of almost Mathieu type. In numerical experiments, with $α$ the inverse golden mean, Fibonacci iterates of maps from the almost Mathieu family exhibit asymptotic scaling behavior that is reminiscent of critical phase transitions. In a restricted setup that is characterized by a symmetry, we prove that critical behavior indeed occurs and is universal in an open neighborhood of the almost Mathieu family. This behavior is governed by a periodic orbit of a renormalization transformation. An extension of this transformation is shown to have a second periodic orbit as well, and we present some evidence that this orbit attracts supercritical almost Mathieu maps. △ Less

Submitted 28 April, 2021; originally announced April 2021.

Comments: For the source code of our programs and possible updates, see https://web.ma.utexas.edu/users/koch/papers/skewunivers/

MSC Class: 37E20

Broadband coupling of fast electrons to high-Q whispering-gallery mode resonators

Authors: Niklas Müller, Vincent Hock, Holger Koch, Nora Bach, Christopher Rathje, Sascha Schäfer

Abstract: Transmission electron microscopy is an excellent experimental tool to study the interaction of free electrons with nanoscale light fields. However, up to now, applying electron microscopy to quantum optical investigations was hampered by the lack of experimental platforms which allow a strong coupling between fast electrons and high-quality resonators. Here, as a first step, we demonstrate the bro… ▽ More Transmission electron microscopy is an excellent experimental tool to study the interaction of free electrons with nanoscale light fields. However, up to now, applying electron microscopy to quantum optical investigations was hampered by the lack of experimental platforms which allow a strong coupling between fast electrons and high-quality resonators. Here, as a first step, we demonstrate the broad-band excitation of optical whispering-gallery modes in silica microresonators by fast electrons. In the emitted coherent cathodoluminescence spectrum, a comb of equidistant peaks is observed, resulting in cavity quality factors larger than 700. These results enable the study of quantum optical phenomena in electron microscopy with potential applications in quantum electron-light metrology. △ Less

Submitted 24 March, 2021; originally announced March 2021.

Comments: 19 pages, 4 figures, 52 references

Mapping of the magnetic field to correct systematic effects in a neutron electric dipole moment experiment

Authors: C. Abel, N. J. Ayres, G. Ban, G. Bison, K. Bodek, V. Bondar, E. Chanel, P. -J. Chiu, B. Clément, C. B. Crawford, M. Daum, S. Emmenegger, L. Ferraris-Bouchez, M. Fertl, P. Flaux, A. Fratangelo, W. C. Griffith, Z. D. Grujić, P. G. Harris, L. Hayen, N. Hild, M. Kasprzak, K. Kirch, P. Knowles, H. -C. Koch , et al. (28 additional authors not shown)

Abstract: Experiments dedicated to the measurement of the electric dipole moment of the neutron require outstanding control of the magnetic field uniformity. The neutron electric dipole moment (nEDM) experiment at the Paul Scherrer Institute uses a 199Hg co-magnetometer to precisely monitor magnetic field variations. This co-magnetometer, in the presence of field non-uniformity, is responsible for the large… ▽ More Experiments dedicated to the measurement of the electric dipole moment of the neutron require outstanding control of the magnetic field uniformity. The neutron electric dipole moment (nEDM) experiment at the Paul Scherrer Institute uses a 199Hg co-magnetometer to precisely monitor magnetic field variations. This co-magnetometer, in the presence of field non-uniformity, is responsible for the largest systematic effect of this measurement. To evaluate and correct that effect, offline measurements of the field non-uniformity were performed during mapping campaigns in 2013, 2014 and 2017. We present the results of these campaigns, and the improvement the correction of this effect brings to the neutron electric dipole moment measurement. △ Less

Submitted 3 May, 2022; v1 submitted 16 March, 2021; originally announced March 2021.

PANDA Phase One

Authors: G. Barucca, F. Davì, G. Lancioni, P. Mengucci, L. Montalto, P. P. Natali, N. Paone, D. Rinaldi, L. Scalise, B. Krusche, M. Steinacher, Z. Liu, C. Liu, B. Liu, X. Shen, S. Sun, G. Zhao, J. Zhao, M. Albrecht, W. Alkakhi, S. Bökelmann, S. Coen, F. Feldbauer, M. Fink, J. Frech , et al. (399 additional authors not shown)

Abstract: The Facility for Antiproton and Ion Research (FAIR) in Darmstadt, Germany, provides unique possibilities for a new generation of hadron-, nuclear- and atomic physics experiments. The future antiProton ANnihilations at DArmstadt (PANDA or $\overline{\rm P}$ANDA) experiment at FAIR will offer a broad physics programme, covering different aspects of the strong interaction. Understanding the latter in… ▽ More The Facility for Antiproton and Ion Research (FAIR) in Darmstadt, Germany, provides unique possibilities for a new generation of hadron-, nuclear- and atomic physics experiments. The future antiProton ANnihilations at DArmstadt (PANDA or $\overline{\rm P}$ANDA) experiment at FAIR will offer a broad physics programme, covering different aspects of the strong interaction. Understanding the latter in the non-perturbative regime remains one of the greatest challenges in contemporary physics. The antiproton-nucleon interaction studied with PANDA provides crucial tests in this area. Furthermore, the high-intensity, low-energy domain of PANDA allows for searches for physics beyond the Standard Model, e.g. through high precision symmetry tests. This paper takes into account a staged approach for the detector setup and for the delivered luminosity from the accelerator. The available detector setup at the time of the delivery of the first antiproton beams in the HESR storage ring is referred to as the \textit{Phase One} setup. The physics programme that is achievable during Phase One is outlined in this paper. △ Less

Submitted 9 June, 2021; v1 submitted 28 January, 2021; originally announced January 2021.

Comments: 35 pages, 15 figures

Journal ref: European Physics Journal A57, 44 (2021)

Combining Multilevel Hartree Fock and Multilevel Coupled Cluster with Molecular Mechanics: a Study of Electronic Excitations in Solutions

Authors: Linda Goletto, Tommaso Giovannini, Sarai D. Folkestad, Henrik Koch

Abstract: We investigate the coupling of different quantum-embedding approaches with a third molecular-mechanics layer, which can be either polarizable or non-polarizable. In particular, such a coupling is discussed for the multilevel families of methods, in which the system is divided into an active and an inactive orbital spaces. The computational cost of the resulting three layers approaches is reduced b… ▽ More We investigate the coupling of different quantum-embedding approaches with a third molecular-mechanics layer, which can be either polarizable or non-polarizable. In particular, such a coupling is discussed for the multilevel families of methods, in which the system is divided into an active and an inactive orbital spaces. The computational cost of the resulting three layers approaches is reduced by treating the long-range interactions at the classical level.The methods developed are tested against the calculation of the excitation energies of molecular systems in aqueous solution, for which an atomistic description of the environment is crucial to correctly reproduce the specific solute-solvent interactions, such as hydrogen bonding. In particular, we present the results obtained for three different moieties, acrolein, pyridine and para-nitroaniline, showing that an almost perfect agreement with experimental data can be achieved when the relevant physico-chemical interactions are included in the modeling of the condensed phase. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: 13 pages and 10 figures

On instability mechanisms for inverse problems

Authors: Herbert Koch, Angkana Rüland, Mikko Salo

Abstract: In this article we present three robust instability mechanisms for linear and nonlinear inverse problems. All of these are based on strong compression properties (in the sense of singular value or entropy number bounds) which we deduce through either strong global smoothing, only weak global smoothing or microlocal smoothing for the corresponding forward operators, respectively. As applications we… ▽ More In this article we present three robust instability mechanisms for linear and nonlinear inverse problems. All of these are based on strong compression properties (in the sense of singular value or entropy number bounds) which we deduce through either strong global smoothing, only weak global smoothing or microlocal smoothing for the corresponding forward operators, respectively. As applications we for instance present new instability arguments for unique continuation, for the backward heat equation and for linear and nonlinear Calderón type problems in general geometries, possibly in the presence of rough coefficients. Our instability mechanisms could also be of interest in the context of control theory, providing estimates on the cost of (approximate) controllability in rather general settings. △ Less

Submitted 2 December, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: 93 pages, 1 figure, comments welcome, corrected some typos

Journal ref: Ars Inveniendi Analytica (2021), Paper No. 7, 93 pp

Intermolecular interactions in optical cavities: an ab initio QED study

Authors: Tor S. Haugland, Christian Schäfer, Enrico Ronca, Angel Rubio, Henrik Koch

Abstract: Intermolecular bonds are weak compared to covalent bonds, but they are strong enough to influence the properties of large molecular systems. In this work, we investigate how strong light-matter coupling inside an optical cavity can modify these intermolecular forces. We perform a detailed comparison between currently available ab initio electron-photon methodologies. The electromagnetic field insi… ▽ More Intermolecular bonds are weak compared to covalent bonds, but they are strong enough to influence the properties of large molecular systems. In this work, we investigate how strong light-matter coupling inside an optical cavity can modify these intermolecular forces. We perform a detailed comparison between currently available ab initio electron-photon methodologies. The electromagnetic field inside the cavity can modulate the ground state properties of weakly bound complexes. Controlling the field polarization, the interactions can be stabilized or destabilized, and electron densities, dipole moments, and polarizabilities can be altered. We demonstrate that electron-photon correlation is fundamental to describe intermolecular interactions in strong light-matter coupling. This work proposes optical cavities as a novel tool to manipulate and control ground state properties, solvent effects, and intermolecular interactions for molecules and materials. △ Less

Submitted 4 December, 2020; v1 submitted 2 December, 2020; originally announced December 2020.

Comments: 11 pages and 9 figures

Transient Resonant Auger-Meitner Spectra of Photoexcited Thymine

Authors: Thomas J. A. Wolf, Alexander C. Paul, Sarai D. Folkestad, Rolf H. Myhre, James P. Cryan, Nora Berrah, Phil H. Bucksbaum, Sonia Coriani, Giacomo Coslovich, Raimund Feifel, Todd J. Martinez, Stefan P. Moeller, Melanie Mucke, Razib Obaid, Oksana Plekan, Richard J. Squibb, Henrik Koch, Markus Gühr

Abstract: We present the first investigation of excited state dynamics by resonant Auger-Meitner spectroscopy (also known as resonant Auger spectroscopy) using the nucleobase thymine as an example. Thymine is photoexcited in the UV and probed with X-ray photon energies at and below the oxygen K-edge. After initial photoexcitation to a ππ* excited state, thymine is known to undergo internal conversion to an… ▽ More We present the first investigation of excited state dynamics by resonant Auger-Meitner spectroscopy (also known as resonant Auger spectroscopy) using the nucleobase thymine as an example. Thymine is photoexcited in the UV and probed with X-ray photon energies at and below the oxygen K-edge. After initial photoexcitation to a ππ* excited state, thymine is known to undergo internal conversion to an nπ* excited state with a strong resonance at the oxygen K-edge, red-shifted from the ground state π* resonances of thymine (see our previous study Wolf et al., Nat. Commun., 2017, 8, 29). We resolve and compare the Auger-Meitner electron spectra associated both with the excited state and ground state resonances, and distinguish participator and spectator decay contributions. Furthermore, we observe simultaneously with the decay of the nπ* state signatures the appearance of additional resonant Auger-Meitner contributions at photon energies between the nπ* state and the ground state resonances. We assign these contributions to population transfer from the nπ* state to a ππ* triplet state via intersystem crossing on the picosecond timescale based on simulations of the X-ray absorption spectra in the vibrationally hot triplet state. Moreover, we identify signatures from the initially excited ππ* singlet state which we have not observed in our previous study. △ Less

Submitted 29 September, 2020; originally announced September 2020.

A Hopf bifurcation in the planar Navier-Stokes equations

Authors: Gianni Arioli, Hans Koch

Abstract: We consider the Navier-Stokes equation for an incompressible viscous fluid on a square, satisfying Navier boundary conditions and being subjected to a time-independent force. As the kinematic viscosity is varied, a branch of stationary solutions is shown to undergo a Hopf bifurcation, where a periodic cycle branches from the stationary solution. Our proof is constructive and uses computer-assisted… ▽ More We consider the Navier-Stokes equation for an incompressible viscous fluid on a square, satisfying Navier boundary conditions and being subjected to a time-independent force. As the kinematic viscosity is varied, a branch of stationary solutions is shown to undergo a Hopf bifurcation, where a periodic cycle branches from the stationary solution. Our proof is constructive and uses computer-assisted estimates. △ Less

Submitted 27 September, 2020; originally announced September 2020.

Comments: 16 pages, 1 figure

A search for neutron to mirror-neutron oscillations

Authors: C. Abel, N. J. Ayres, G. Ban, G. Bison, K. Bodek, V. Bondar, E. Chanel, P. -J. Chiu, C. Crawford, M. Daum, R. T. Dinani, S. Emmenegger, P. Flaux, L. Ferraris-Bouchez, W. C. Griffith, Z. D. Grujic, N. Hild, K. Kirch, H. -C. Koch, P. A. Koss, A. Kozela, J. Krempel, B. Lauss, T. Lefort, A. Leredde , et al. (18 additional authors not shown)

Abstract: It has been proposed that there could be a mirror copy of the standard model particles, restoring the parity symmetry in the weak interaction on the global level. Oscillations between a neutral standard model particle, such as the neutron, and its mirror counterpart could potentially answer various standing issues in physics today. Astrophysical studies and terrestrial experiments led by ultracold… ▽ More It has been proposed that there could be a mirror copy of the standard model particles, restoring the parity symmetry in the weak interaction on the global level. Oscillations between a neutral standard model particle, such as the neutron, and its mirror counterpart could potentially answer various standing issues in physics today. Astrophysical studies and terrestrial experiments led by ultracold neutron storage measurements have investigated neutron to mirror-neutron oscillations and imposed constraints on the theoretical parameters. Recently, further analysis of these ultracold neutron storage experiments has yielded statistically significant anomalous signals that may be interpreted as neutron to mirror-neutron oscillations, assuming nonzero mirror magnetic fields. The neutron electric dipole moment collaboration performed a dedicated search at the Paul Scherrer Institute and found no evidence of neutron to mirror-neutron oscillations. Thereby, the following new lower limits on the oscillation time were obtained: $τ_{nn'} > 352~$s at $B'=0$ (95% C.L.), $τ_{nn'} > 6~\text{s}$ for all $0.4~μ\text{T}<B'<25.7~μ\text{T}$ (95% C.L.), and $τ_{nn'}/\sqrt{\cosβ}>9~\text{s}$ for all $5.0~μ\text{T}<B'<25.4~μ\text{T}$ (95% C.L.), where $β$ is the fixed angle between the applied magnetic field and the local mirror magnetic field which is assumed to be bound to the Earth. These new constraints are the best measured so far around $B'\sim10~μ$T, and $B'\sim20~μ$T. △ Less

Submitted 19 November, 2020; v1 submitted 23 September, 2020; originally announced September 2020.

Comments: Submitted to Physics Letters B, 9 pages, 3 figures