1
|
Woźniak AP, Moszyński R. Modeling of High-Harmonic Generation in the C 60 Fullerene Using Ab Initio, DFT-Based, and Semiempirical Methods. J Phys Chem A 2024; 128:2683-2702. [PMID: 38534023 PMCID: PMC11017253 DOI: 10.1021/acs.jpca.3c07865] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 03/04/2024] [Accepted: 03/13/2024] [Indexed: 03/28/2024]
Abstract
We report calculations of the high-harmonic generation spectra of the C60 fullerene molecule carried out by employing a diverse set of real-time time-dependent quantum chemical methods. All methodologies involve expanding the propagated electronic wave function in bases consisting of the ground and singly excited time-independent eigenstates obtained through the solution of the corresponding linear-response equations. We identify the correlation and exchange effect in the spectra by comparing the results from methods relying on the Hartree-Fock reference determinant with those obtained using approaches based on the density functional theory with different exchange-correlation functionals. The effect of the full random-phase approximation treatment of the excited electronic states is also analyzed and compared with the configuration interaction singles and the Tamm-Dancoff approximation. We also showcase the fact that the real-time extension of the semiempirical method INDO/S can be effectively applied for an approximate description of laser-driven dynamics in large systems.
Collapse
Affiliation(s)
| | - Robert Moszyński
- Faculty of Chemistry, University
of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
2
|
Morassut C, Ravindran A, Ciavardini A, Luppi E, De Ninno G, Coccia E. High-Harmonic Generation Spectroscopy of Gas-Phase Bromoform. J Phys Chem A 2024; 128:2015-2024. [PMID: 38469750 DOI: 10.1021/acs.jpca.3c07699] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/13/2024]
Abstract
High-Harmonic Generation (HHG) spectra of randomly aligned bromoform (CHBr3) molecules have been experimentally measured and theoretically simulated at various laser pulse intensities. From the experiments, we obtained a significant number of harmonics that goes beyond the cutoff limit predicted by the three-step model (3SM) with ionization from HOMO. To interpret the experiment, we resorted to real-time time-dependent configuration interaction with single excitations. We found that electronic bound states provide an appreciable contribution to the harmonics. More in detail, we analyzed the electron dynamics by decomposing the HHG signal in terms of single molecular-orbital contributions, to explain the appearance of harmonics around 20-30 eV beyond the expected cutoff due to HOMO. HHG spectra can be therefore explained by considering the contribution at high energy of HOMO-6 and HOMO-9, thus indicating a complex multiple-orbital strong-field dynamics. However, even though the presence of the bromoform cation should be not enough to produce such a signal, we could not exclude a priori that the origin of harmonics in the H29-H45 to be due to the cation, which has more energetic ionization channels.
Collapse
Affiliation(s)
- Chiara Morassut
- Laboratoire de Chimie Théorique, Sorbonne Université, CNRS, Paris F-75005, France
- Department of Chemical and Pharmaceutical Sciences, University of Trieste, via L. Giorgieri 1, 34127 Trieste, Italy
| | - Arun Ravindran
- Laboratory of Quantum Optics, University of Nova Gorica, Si-5270 Ajdovščina, Slovenija
| | - Alessandra Ciavardini
- Laboratory of Quantum Optics, University of Nova Gorica, Si-5270 Ajdovščina, Slovenija
| | - Eleonora Luppi
- Laboratoire de Chimie Théorique, Sorbonne Université, CNRS, Paris F-75005, France
| | - Giovanni De Ninno
- Laboratory of Quantum Optics, University of Nova Gorica, Si-5270 Ajdovščina, Slovenija
- Elettra-Sincrotrone Trieste S.C.p.A., S.S. 14-km 163.5 in AREA Science Park, 34149 Basovizza, Trieste, Italy
| | - Emanuele Coccia
- Department of Chemical and Pharmaceutical Sciences, University of Trieste, via L. Giorgieri 1, 34127 Trieste, Italy
| |
Collapse
|
3
|
Morassut C, Coccia E, Luppi E. Quantitative performance analysis and comparison of optimal-continuum Gaussian basis sets for high-harmonic generation spectra. J Chem Phys 2023; 159:124108. [PMID: 38127378 DOI: 10.1063/5.0153825] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2023] [Accepted: 08/24/2023] [Indexed: 12/23/2023] Open
Abstract
Quantum-chemistry methods in the time domain with Gaussian basis sets are increasingly used to compute high-harmonic generation (HHG) spectra of atomic and molecular systems. The quality of these approaches is limited by the accuracy of Gaussian basis sets to describe continuum energy states. In the literature, optimal-continuum Gaussian basis sets have been proposed: Kaufmann et al. [J. Phys. B: At., Mol. Opt. Phys. 22, 2223 (1989)], Woźniak et al. [J. Chem. Phys. 154, 094111 (2021)], Nestmann and Peyerimhoff [J. Phys. B: At., Mol. Opt. Phys. 23, L773 (1990)], Faure et al. [Comput. Phys. Commun. 144, 224 (2002)], and Krause et al. [J. Chem. Phys. 140, 174113 (2014)]. In this work, we have compared the performances of these basis sets to simulate HHG spectra of H atom at different laser intensities. We have also investigated different strategies to balance basis sets with these continuum functions, together with the role of angular momentum. To quantify the performance of the different basis sets, we introduce local and global HHG descriptors. Comparisons with the grid and exact calculations are also provided.
Collapse
Affiliation(s)
- C Morassut
- Laboratoire de Chimie Théorique, Sorbonne Université, CNRS, Paris F-75005, France
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università di Trieste, Via L. Giorgieri 1, 34127 Trieste, Italy
| | - E Coccia
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università di Trieste, Via L. Giorgieri 1, 34127 Trieste, Italy
| | - E Luppi
- Laboratoire de Chimie Théorique, Sorbonne Université, CNRS, Paris F-75005, France
| |
Collapse
|
4
|
Luppi E, Coccia E. Role of Inner Molecular Orbitals in High-Harmonic Generation Spectra of Aligned Uracil. J Phys Chem A 2023; 127:7335-7343. [PMID: 37640677 DOI: 10.1021/acs.jpca.3c03990] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/31/2023]
Abstract
In this work, we decompose the high-harmonic generation (HHG) signal of aligned gas-phase uracil into single molecular-orbital (MO) contributions. We compute HHG spectra for a pulse linearly polarized perpendicular to the molecular plane, with an intensity of 0.6 and 0.85 × 1014 W/cm2 and a wavelength of 800 nm. We use the real-time time-dependent Configuration Interaction with singles method, coupled to a Gaussian-based representation of the time-dependent wavefunction. The strong-field dynamics is affected by the energy of the ionization/recombination channels and by the coupling between the orbital symmetry and laser polarization. In the configuration studied here, we expect that π-type MOs favorably couple with the incoming pulse and play a substantial role in generating the HHG spectrum. Indeed, we show that HOMO, HOMO - 1, and HOMO - 4, which all are π-like, determine the intensity of harmonic peaks at different energies, while HOMO - 2 and HOMO - 3 provide a smaller contribution. It is worth mentioning that HOMO - 4 produces a stronger signal than that from HOMO - 1, even though the corresponding ionization energy, in an one-electron picture, is around 2.5 eV larger and more than 4 eV larger than the HOMO one.
Collapse
Affiliation(s)
- Eleonora Luppi
- Laboratoire de Chimie Théorique, Sorbonne Université, CNRS, Paris F-75005, France
| | - Emanuele Coccia
- Department of Chemical and Pharmaceutical Sciences, University of Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| |
Collapse
|
5
|
Monti M, Stener M, Coccia E. Electronic circular dichroism from real-time propagation in state space. J Chem Phys 2023; 158:084102. [PMID: 36859092 DOI: 10.1063/5.0136392] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
In this paper, we propose to compute the electronic circular dichroism (ECD) spectra of chiral molecules using a real-time propagation of the time-dependent Schrödinger equation (TDSE) in the space of electronic field-free eigenstates, by coupling TDSE with a given treatment of the electronic structure of the target. The time-dependent induced magnetic moment is used to compute the ECD spectrum from an explicit electric perturbation. The full matrix representing the transition magnetic moment in the space of electronic states is generated from that among pairs of molecular orbitals. In the present work, we show the ECD spectra of methyloxirane, of several conformers of L-alanine, and of the Λ-Co(acac)3 complex, computed from a singly excited ansatz of time-dependent density functional theory eigenstates. The time-domain ECD spectra properly reproduce the frequency-domain ones obtained in the linear-response regime and quantitatively agree with the available experimental data. Moreover, the time-domain approach to ECD allows us to naturally go beyond the ground-state rotationally averaged ECD spectrum, which is the standard outcome of the linear-response theory, e.g., by computing the ECD spectra from electronic excited states.
Collapse
Affiliation(s)
- M Monti
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università di Trieste, Via L. Giorgieri 1, 34127 Trieste, Italy
| | - M Stener
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università di Trieste, Via L. Giorgieri 1, 34127 Trieste, Italy
| | - E Coccia
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università di Trieste, Via L. Giorgieri 1, 34127 Trieste, Italy
| |
Collapse
|
6
|
Woźniak AP, Lewenstein M, Moszyński R. Exploring the attosecond laser-driven electron dynamics in the hydrogen molecule with different real-time time-dependent configuration interaction approaches. ADVANCES IN QUANTUM CHEMISTRY 2023. [DOI: 10.1016/bs.aiq.2023.01.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
|
7
|
A TD-CIS study of high-harmonic generation of uracil cation fragments. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
8
|
Zhu Y, Herbert JM. High harmonic spectra computed using time-dependent Kohn-Sham theory with Gaussian orbitals and a complex absorbing potential. J Chem Phys 2022; 156:204123. [PMID: 35649850 DOI: 10.1063/5.0079910] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
High harmonic spectra for H2 and H2 + are simulated by solving the time-dependent Kohn-Sham equation in the presence of a strong laser field using an atom-centered Gaussian representation of the density and a complex absorbing potential. The latter serves to mitigate artifacts associated with the finite extent of the basis functions, including spurious reflection of the outgoing electronic wave packet. Interference between the outgoing and reflected waves manifests as peak broadening in the spectrum as well as the appearance of spurious high-energy peaks after the harmonic progression has terminated. We demonstrate that well-resolved spectra can be obtained through the use of an atom-centered absorbing potential. As compared to grid-based algorithms, the present approach is more readily extensible to larger molecules.
Collapse
Affiliation(s)
- Ying Zhu
- Department of Chemistry and Biochemistry, and Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M Herbert
- Department of Chemistry and Biochemistry, and Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
9
|
Wozniak AP, Przybytek M, Lewenstein M, Moszynski R. Effects of electronic correlation on the high harmonic generation in helium: a time-dependent configuration interaction singles vs time-dependent full configuration interaction study. J Chem Phys 2022; 156:174106. [DOI: 10.1063/5.0087384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we investigate the effects of full electronic correlation on the high harmonic generation in the helium atom subjected to laser pulses of extremely high intensity. To do this, we perform real-time propagations of the helium atom wavefunction using quantum chemistry methods coupled to Gaussian basis sets. The calculations are done within the real-time time-dependent configuration interaction framework, at two levels of theory: time-dependent configuration interation with single excitations (TD-CIS, uncorrelated method) and time-dependent full configuration interaction (TD-FCI, fully correlated method). The electronic wavefunction is expanded in Dunning basis sets supplemented with functions adapted to describing highly excited and continuum states. We also compare the TD-CI results with grid-based propagations of the helium atom within the single-active-electron approximation. Our results show that when including the dynamical electron correlation, a noticeable improvement to the description of HHG can be achieved, in terms of e.g. a more constant intensity in the lower energy part of the harmonic plateau. However, such effects can be captured only if the basis set used suffices to reproduce the most basic features, such as the HHG cutoff position, at the uncorrelated level of theory.
Collapse
|
10
|
Witzorky C, Paramonov G, Bouakline F, Jaquet R, Saalfrank P, Klamroth T. Gaussian-Type Orbital Calculations for High Harmonic Generation in Vibrating Molecules: Benchmarks for H 2. J Chem Theory Comput 2021; 17:7353-7365. [PMID: 34747605 DOI: 10.1021/acs.jctc.1c00837] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The response of the hydrogen molecular ion, H2+, to few-cycle laser pulses of different intensities is simulated. To treat the coupled electron-nuclear motion, we use adiabatic potentials computed with Gaussian-type basis sets together with a heuristic ionization model for the electron and a grid representation for the nuclei. Using this mixed-basis approach, the time-dependent Schrödinger equation is solved, either within the Born-Oppenheimer approximation or with nonadiabatic couplings included. The dipole response spectra are compared to all-grid-based solutions for the three-body problem, which we take as a reference to benchmark the Gaussian-type basis set approaches. Also, calculations employing the fixed-nuclei approximation are performed, to quantify effects due to nuclear motion. For low intensities and small ionization probabilities, we get excellent agreement of the dynamics using Gaussian-type basis sets with the all-grid solutions. Our investigations suggest that high harmonic generation (HHG) and high-frequency response, in general, can be reliably modeled using Gaussian-type basis sets for the electrons for not too high harmonics. Further, nuclear motion destroys electronic coherences in the response spectra even on the time scale of about 30 fs and affects HHG intensities, which reflect the electron dynamics occurring on the attosecond time scale. For the present system, non-Born-Oppenheimer effects are small. The Gaussian-based, nonadiabatically coupled, time-dependent multisurface approach to treat quantum electron-nuclear motion beyond the non-Born-Oppenheimer approximation can be easily extended to approximate wavefunction methods, such as time-dependent configuration interaction singles (TD-CIS), for systems where no benchmarks are available.
Collapse
Affiliation(s)
- Christoph Witzorky
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| | - Guennaddi Paramonov
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| | - Foudhil Bouakline
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| | - Ralph Jaquet
- Theoretische Chemie, Universität Siegen, D-57068 Siegen, Germany
| | - Peter Saalfrank
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| | - Tillmann Klamroth
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| |
Collapse
|
11
|
Coccia E, Luppi E. Time-dependent ab initioapproaches for high-harmonic generation spectroscopy. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 34:073001. [PMID: 34731835 DOI: 10.1088/1361-648x/ac3608] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 11/03/2021] [Indexed: 06/13/2023]
Abstract
High-harmonic generation (HHG) is a nonlinear physical process used for the production of ultrashort pulses in XUV region, which are then used for investigating ultrafast phenomena in time-resolved spectroscopies. Moreover, HHG signal itself encodes information on electronic structure and dynamics of the target, possibly coupled to the nuclear degrees of freedom. Investigating HHG signal leads to HHG spectroscopy, which is applied to atoms, molecules, solids and recently also to liquids. Analysing the number of generated harmonics, their intensity and shape gives a detailed insight of, e.g., ionisation and recombination channels occurring in the strong-field dynamics. A number of valuable theoretical models has been developed over the years to explain and interpret HHG features, with the three-step model being the most known one. Originally, these models neglect the complexity of the propagating electronic wavefunction, by only using an approximated formulation of ground and continuum states. Many effects unravelled by HHG spectroscopy are instead due to electron correlation effects, quantum interference, and Rydberg-state contributions, which are all properly captured by anab initioelectronic-structure approach. In this review we have collected recent advances in modelling HHG by means ofab initiotime-dependent approaches relying on the propagation of the time-dependent Schrödinger equation (or derived equations) in presence of a very intense electromagnetic field. We limit ourselves to gas-phase atomic and molecular targets, and to solids. We focus on the various levels of theory employed for describing the electronic structure of the target, coupled with strong-field dynamics and ionisation approaches, and on the basis used to represent electronic states. Selected applications and perspectives for future developments are also given.
Collapse
Affiliation(s)
- Emanuele Coccia
- Dipartimento di Scienze Chimiche e Farmaceutiche, University of Trieste, via Giorgieri 1, 34127 Trieste, Italy
| | - Eleonora Luppi
- Sorbonne Universités, UPMC Univ Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005 Paris, France
- CNRS, UMR 7616, Laboratoire de Chimie Théorique, F-75005 Paris, France
| |
Collapse
|
12
|
Bedurke F, Klamroth T, Saalfrank P. Many-electron dynamics in laser-driven molecules: wavefunction theory vs. density functional theory. Phys Chem Chem Phys 2021; 23:13544-13560. [PMID: 34105544 DOI: 10.1039/d1cp01100f] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
With recent experimental advances in laser-driven electron dynamics in polyatomic molecules, the need arises for their reliable theoretical modelling. Among efficient, yet fairly accurate methods for many-electron dynamics are Time-Dependent Configuration Interaction Singles (TD-CIS) (a Wave Function Theory (WFT) method), and Real-Time Time-Dependent Density Functional Theory (RT-TD-DFT), respectively. Here we compare TD-CIS combined with extended Atomic Orbital (AO) bases, TD-CIS/AO, with RT-TD-DFT in a grid representation of the Kohn-Sham orbitals, RT-TD-DFT/Grid. Possible ionization losses are treated by complex absorbing potentials in energy space (for TD-CIS/AO) or real space (for RT-TD-DFT), respectively. The comparison is made for two test cases: (i) state-to-state transitions using resonant lasers (π-pulses), i.e., bound electron motion, and (ii) large-amplitude electron motion leading to High Harmonic Generation (HHG). Test systems are a H2 molecule and cis- and trans-1,2-dichlorethene, C2H2Cl2, (DCE). From time-dependent electronic energies, dipole moments and from HHG spectra, the following observations are made: first, for bound state-to-state transitions enforced by π-pulses, TD-CIS nicely accounts for the expected population inversion in contrast to RT-TD-DFT, in agreement with earlier findings. Secondly, when using laser pulses under non-resonant conditions, dipole moments and lower harmonics in HHG spectra are obtained by TD-CIS/AO which are in good agreement with those obtained with RT-TD-DFT/Grid. Deviations become larger for higher harmonics and at low laser intensities, i.e., for low-intensity HHG signals. We also carefully test effects of basis sets for TD-CIS/AO and grid size for RT-TD-DFT/Grid, different exchange-correlation functionals in RT-TD-DFT, and absorbing boundaries. Finally, for the present examples, TD-CIS/AO is observed to be at least an order of magnitude more computationally efficient than RT-TD-DFT/Grid.
Collapse
Affiliation(s)
- Florian Bedurke
- Universität Potsdam, Institut für Chemie, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany.
| | - Tillmann Klamroth
- Universität Potsdam, Institut für Chemie, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany.
| | - Peter Saalfrank
- Universität Potsdam, Institut für Chemie, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany.
| |
Collapse
|
13
|
Woźniak AP, Lesiuk M, Przybytek M, Efimov DK, Prauzner-Bechcicki JS, Mandrysz M, Ciappina M, Pisanty E, Zakrzewski J, Lewenstein M, Moszyński R. A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation. J Chem Phys 2021; 154:094111. [PMID: 33685145 DOI: 10.1063/5.0040879] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A precise understanding of mechanisms governing the dynamics of electrons in atoms and molecules subjected to intense laser fields has a key importance for the description of attosecond processes such as the high-harmonic generation and ionization. From the theoretical point of view, this is still a challenging task, as new approaches to solve the time-dependent Schrödinger equation with both good accuracy and efficiency are still emerging. Until recently, the purely numerical methods of real-time propagation of the wavefunction using finite grids have been frequently and successfully used to capture the electron dynamics in small one- or two-electron systems. However, as the main focus of attoscience shifts toward many-electron systems, such techniques are no longer effective and need to be replaced by more approximate but computationally efficient ones. In this paper, we explore the increasingly popular method of expanding the wavefunction of the examined system into a linear combination of atomic orbitals and present a novel systematic scheme for constructing an optimal Gaussian basis set suitable for the description of excited and continuum atomic or molecular states. We analyze the performance of the proposed basis sets by carrying out a series of time-dependent configuration interaction calculations for the hydrogen atom in fields of intensity varying from 5 × 1013 W/cm2 to 5 × 1014 W/cm2. We also compare the results with the data obtained using Gaussian basis sets proposed previously by other authors.
Collapse
Affiliation(s)
| | - Michał Lesiuk
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Przybytek
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Dmitry K Efimov
- Institute of Theoretical Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Jakub S Prauzner-Bechcicki
- Marian Smoluchowski Institute of Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Michał Mandrysz
- Institute of Theoretical Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Marcelo Ciappina
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860, Castelldefels, Barcelona, Spain
| | - Emilio Pisanty
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860, Castelldefels, Barcelona, Spain
| | - Jakub Zakrzewski
- Institute of Theoretical Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Maciej Lewenstein
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860, Castelldefels, Barcelona, Spain
| | - Robert Moszyński
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
14
|
Luppi E, Coccia E. Probing the molecular frame of uracil and thymine with high-harmonic generation spectroscopy. Phys Chem Chem Phys 2021; 23:3729-3738. [PMID: 33395454 DOI: 10.1039/d0cp05559j] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work, we present the computed high-harmonic generation (HHG) spectra of uracil and thymine molecules, by means of the real-time time-dependent formulation of Gaussian-based configuration interaction with single excitations (RT-TD-CIS). According to the experimental work [Hutchison et al., Comparison of high-order harmonic generation in uracil and thymine ablation plumes, Phys. Chem. Chem. Phys., 2013, 15, 12308], a pulse wavelength of 780 nm has been used, together with an intensity of 1014 W cm-2 and a pulse duration of 23 optical cycles. In order to examine the effect of pulse polarisation, rotationally averaged (to mimic the gas-phase sample) and single-polarisations have been computed for both molecules. Our results show that the HHG signal for both molecules possibly originates from different ionisation channels, involving HOMO, HOMO-1, HOMO-2 and HOMO-3 orbitals, which lie within 4 eV. We characterize the HHG spectrum of thymine, supporting the idea that the absence of the thymine signal in the original work does not depend on the single-molecule behaviour. Present results for uracil are consistent with the experimental data. Moreover, we have observed that states below and above the chosen ionisation threshold provide different contributions to the HHG spectrum in averaged and single-polarisation calculations.
Collapse
Affiliation(s)
- Eleonora Luppi
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France and CNRS, UMR 7616, Laboratoire de Chimie Théorique, F-75005 Paris, France
| | - Emanuele Coccia
- Dipartimento di Scienze Chimiche e Farmaceutiche, via Giorgieri 1, 34127 Trieste, Italy.
| |
Collapse
|
15
|
Pauletti CF, Coccia E, Luppi E. Role of exchange and correlation in high-harmonic generation spectra of H 2, N 2, and CO 2: Real-time time-dependent electronic-structure approaches. J Chem Phys 2021; 154:014101. [PMID: 33412879 DOI: 10.1063/5.0033072] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
This study arises from the attempt to answer the following question: how different descriptions of electronic exchange and correlation affect the high-harmonic generation (HHG) spectroscopy of H2, N2, and CO2 molecules? We compare HHG spectra for H2, N2, and CO2 with different ab initio electronic structure methods: real-time time-dependent configuration interaction and real-time time-dependent density functional theory (RT-TDDFT) using truncated basis sets composed of correlated wave functions expanded on Gaussian basis sets. In the framework of RT-TDDFT, we employ Perdew-Burke-Ernzerhof (PBE) and long-range corrected Perdew-Burke-Ernzerhof (LC-ωPBE) functionals. We study HHG spectroscopy by disentangling the effect of electronic exchange and correlation. We first analyze the electronic exchange alone, and in the case of RT-TDDFT with LC-ωPBE, we use ω = 0.3 and ω = 0.4 to tune the percentage of long-range Hartree-Fock exchange and short-range exchange PBE. Then, we added the correlation as described by the PBE functional. All the methods give very similar HHG spectra, and they seem not to be particularly sensitive to the different description of exchange and correlation or to the correct asymptotic behavior of the Coulomb potential. Despite this general trend, some differences are found in the region connecting the cutoff and the background. Here, the harmonics can be resolved with different accuracy depending on the theoretical schemes used. We believe that the investigation of the molecular continuum and its coupling with strong fields merits further theoretical investigations in the near future.
Collapse
Affiliation(s)
| | - Emanuele Coccia
- Dipartimento di Scienze Chimiche e Farmaceutiche, Via Giorgieri 1, Trieste Italy
| | - Eleonora Luppi
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
| |
Collapse
|
16
|
Affiliation(s)
- E. Coccia
- Dipartimento di Scienze Chimiche e Farmaceutiche, University of Trieste, Trieste, Italy
| |
Collapse
|