High-Harmonic Generation Spectroscopy of Gas-Phase Bromoform.

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.

Use of Gaussian-Type Functions for Describing Fast Ion-Matter Irradiation with Time-Dependent Density Functional Theory.

The electronic stopping power is an observable property that quantifies the ability of swift ions to penetrate matter to transfer energy to the electron cloud. The recent literature has proven the value of Real-Time Time-Dependent Density Functional Theory to accurately evaluate this property from first-principles, but questions remain regarding the capability of computer codes relying on atom-centered basis functions to capture the physics at play. In this Perspective, we draw attention to the fact that irradiation by swift ions triggers electron emission into the continuum, especially at the Bragg peak. We investigate the ability of Gaussian atomic orbitals (AOC), which were fitted to mimic continuum wave functions, to improve electronic stopping power predictions. AOC are added to standard correlation-consistent basis sets or STO minimal basis sets. Our benchmarks for water irradiation by fast protons clearly advocate for the use of AOC, especially near the Bragg peak. We show that AOC only need to be placed on the molecules struck by the ion. The number of AOC that are added to the usual basis set is relatively small compared to the total number of atomic orbitals, making the use of such a basis set an excellent choice from a computational cost point of view. The optimum basis set combination is applied for the calculation of the stopping power of a proton in water with encouraging agreement with experimental data.

In silico design of bio-marker detection fluorescent probes.

Fluorescent probes capable of sensing the biological medium are of utmost importance in medical diagnostics. However, the optical spectrum of such probes needs to be tuned with care for compatibility with living tissues. More specifically, fluorescent bioprobes must be adjusted so as to avoid light interference with pigments (e.g. hemoglobin), tissue photodamage, scattering of the emitted light, and autofluorescence. This leads to two important conditions on the optical spectrum of the probes. On the one hand, the emission wavelength must be in an optical window of 650 to 950 nm. On the other hand, the Stokes shift must be large, ideally greater than 150 nm. In this paper, we showcase the in-silico design of potential fluorescent biomarkers fulfilling these two conditions by means of heteroatomic substitution and conjugation on a 1,2,4-triazole core initially far away from biological standards.

Role of Inner Molecular Orbitals in High-Harmonic Generation Spectra of Aligned Uracil.

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.

Photoionization and core resonances from range-separated time-dependent density-functional theory for open-shell states: Example of the lithium atom.

We consider the calculations of photoionization spectra and core resonances of open-shell systems using range-separated time-dependent density-functional theory. Specifically, we use the time-dependent range-separated hybrid (TDRSH) scheme, combining a long-range Hartree-Fock exchange potential and kernel with a short-range potential and kernel from a local density-functional approximation, and the time-dependent locally range-separated hybrid (TDLRSH) scheme, which uses a local range-separation parameter. To efficiently perform the calculations, we formulate a spin-unrestricted linear-response Sternheimer approach in a non-orthogonal B-spline basis set using appropriate frequency-dependent boundary conditions. We illustrate this approach on the Li atom, which suggests that TDRSH and TDLRSH are adequate simple methods for estimating the single-electron photoionization spectra of open-shell systems.

Insight into the inclusion of heteroatom impurities in silicon structures.

The bonding properties of the tilt boundary in poly-silicon and the effect of interstitial impurities are investigated by first-principles. In order to obtain thorough information on the nature of chemical bondings in these solid systems, an accurate topological analysis is performed, through partitioning of the electron localization function. Although the mechanism of segregation of single light impurities, such as carbon, nitrogen, and oxygen, in Si-based systems is known, it is only in the presence of multiple segregations that the distinctive structures of the various interstitial impurities emerge. The structural analysis of the modified Si systems and the comparison with the corresponding molecular structure within these solid phases provide an adequate description of interesting properties, for which bond charges provide more insight than bond length. It is shown that, in the presence of isovalent carbon, all systems try to preserve the tetrahedral coordination; on the contrary, trivalent nitrogen induces a strong local distortion to fit in the tetrahedral Si matrix while oxygen is the impurity that segregates more easily and more regularly. This work shows that impurities lead to local distortions and shows how the electron distribution rearranges to smooth these distortions. Overall, it shows how the analysis of bonds and their correlation with energetics and electronic structure is of fundamental importance for the understanding of the defects induced properties and of the basic mechanisms that influence them.

Photoionization and core resonances from range-separated density-functional theory: General formalism and example of the beryllium atom.

We explore the merits of linear-response range-separated time-dependent density-functional theory (TDDFT) for the calculation of photoionization spectra. We consider two variants of range-separated TDDFT, namely, the time-dependent range-separated hybrid (TDRSH) scheme, which uses a global range-separation parameter, and the time-dependent locally range-separated hybrid (TDLRSH), which uses a local range-separation parameter, and compare with standard time-dependent local-density approximation (TDLDA) and time-dependent Hartree-Fock (TDHF). We show how to calculate photoionization spectra with these methods using the Sternheimer approach formulated in a non-orthogonal B-spline basis set with appropriate frequency-dependent boundary conditions. We illustrate these methods on the photoionization spectrum of the Be atom, focusing, in particular, on the core resonances. Both the TDRSH and TDLRSH photoionization spectra are found to constitute a large improvement over the TDLDA photoionization spectrum and a more modest improvement over the TDHF photoionization spectrum.

Time-dependentab initioapproaches for high-harmonic generation spectroscopy.

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.

The role of Si vacancies in the segregation of O, C, and N at silicon grain boundaries: An ab initio study.

Grain boundaries (GBs) are defects originating in multi-crystalline silicon during crystal growth for device Si solar cell fabrication. The presence of GBs changes the coordination of Si, making it advantageous for charge carriers to recombine, which brings a significant reduction of carrier lifetimes. Therefore, GBs can be highly detrimental for device performances. Furthermore, GBs easily form vacancies with deep defect electronic states and are also preferential segregation sites for various impurity species, such as C, N, and O. We studied from first principles the correlation between structural, energetics, and electronic properties of the Σ3{111} Si GB with and without vacancies, and the segregation of C, N, and O atoms. C and O atoms strongly increase their ability to segregate when vacancies are present. However, the electronic properties of the Σ3{111} Si GB are not affected by the presence of O, while they can strongly change in the case of C. For N atoms, it is not possible to find a clear trend in the energetics and electronic properties both with and without vacancies in the GB. In fact, as N is not isovalent with Si, as C and O, it is more flexible in finding new chemical arrangements in the GB structure. This implies a stronger difficulty in controlling the properties of the material in the presence of N impurity atoms compared to C and O impurities.

Probing the molecular frame of uracil and thymine with high-harmonic generation spectroscopy.

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.

Role of exchange and correlation in high-harmonic generation spectra of H2, N2, and CO2: Real-time time-dependent electronic-structure approaches.

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.

Long-range corrected exchange-correlation kernels to describe excitons in second-harmonic generation.

We investigate the role of excitons in second-harmonic generation (SHG) through the long-range corrected (LRC) exchange-correlation kernels: empirical LRC, Bootstrap, and jellium-with-a-gap model. We calculate the macroscopic second-order frequency-dependent susceptibility χ(2). We also present the frequency-dependent macroscopic dielectric function ϵM which is a fundamental quantity in the theoretical derivation of χ(2). We assess the role of the long-range kernels in describing excitons in materials with different symmetry types: cubic zincblende, hexagonal wurtzite, and tetragonal symmetry. Our studies indicate that excitons play an important role in χ(2) bringing a strong enhancement of the SHG signal. Moreover, we found that the SHG enhancement follows a simple trend determined by the magnitude of the long-range corrected α-parameter. This trend is material dependent.

Linear-response range-separated density-functional theory for atomic photoexcitation and photoionization spectra.

We investigate the performance of the range-separated hybrid (RSH) scheme, which combines long-range Hartree-Fock (HF) and a short-range density-functional approximation (DFA), for calculating the photoexcitation/photoionization spectra of the H and He atoms, using a B-spline basis set in order to correctly describe the continuum part of the spectra. The study of these simple systems allows us to quantify the influence on the spectra of the errors coming from the short-range exchange-correlation DFA and from the missing long-range correlation in the RSH scheme. We study the differences using the long-range HF exchange (nonlocal) potential and the long-range exact exchange (local) potential. Contrary to the former, the latter supports a series of Rydberg states and gives reasonable photoexcitation/photoionization spectra, even without applying linear-response theory. The most accurate spectra are obtained with the linear-response time-dependent RSH (TDRSH) scheme. In particular, for the He atom at the optimal value of the range-separation parameter, TDRSH gives slightly more accurate photoexcitation and photoionization spectra than the standard linear-response time-dependent HF. More generally, this work shows the potential of range-separated density-functional theory for calculating linear and nonlinear optical properties involving continuum states.

Optimal Basis Set for Electron Dynamics in Strong Laser Fields: The case of Molecular Ion H2.

A clear understanding of the mechanisms that control the electron dynamics in a strong laser field is still a challenge that requires interpretation by advanced theory. Development of accurate theoretical and computational methods, able to provide a precise treatment of the fundamental processes generated in the strong field regime, is therefore crucial. A central aspect is the choice of the basis for the wave function expansion. Accuracy in describing multiphoton processes is strictly related to the intrinsic properties of the basis, such as numerical convergence, computational cost, and representation of the continuum. By explicitly solving the 1D and 3D time-dependent Schrödinger equation for H2+ in the presence of an intense electric field, we explore the numerical performance of using a real-space grid, a B-spline basis, and a Gaussian basis (improved by optimal Gaussian functions for the continuum). We analyze the performance of the three bases for high-harmonic generation and above-threshold ionization for H2+. In particular, for high-harmonic generation, the capability of the basis to reproduce the two-center interference and the hyper-Raman phenomena is investigated.

Ab initio lifetime correction to scattering states for time-dependent electronic-structure calculations with incomplete basis sets.

We propose a method for obtaining effective lifetimes of scattering electronic states for avoiding the artificial confinement of the wave function due to the use of incomplete basis sets in time-dependent electronic-structure calculations of atoms and molecules. In this method, using a fitting procedure, the lifetimes are extracted from the spatial asymptotic decay of the approximate scattering wave functions obtained with a given basis set. The method is based on a rigorous analysis of the complex-energy solutions of the Schrödinger equation. It gives lifetimes adapted to any given basis set without using any empirical parameters. The method can be considered as an ab initio version of the heuristic lifetime model of Klinkusch et al. [J. Chem. Phys. 131, 114304 (2009)]. The method is validated on H and He atoms using Gaussian-type basis sets for the calculation of high-harmonic-generation spectra.

Self-consistent double-hybrid density-functional theory using the optimized-effective-potential method.

We introduce an orbital-optimized double-hybrid (DH) scheme using the optimized-effective-potential (OEP) method. The orbitals are optimized using a local potential corresponding to the complete exchange-correlation energy expression including the second-order Møller-Plesset correlation contribution. We have implemented a one-parameter version of this OEP-based self-consistent DH scheme using the BLYP density-functional approximation and compared it to the corresponding non-self-consistent DH scheme for calculations on a few closed-shell atoms and molecules. While the OEP-based self-consistency does not provide any improvement for the calculations of ground-state total energies and ionization potentials, it does improve the accuracy of electron affinities and restores the meaning of the LUMO orbital energy as being connected to a neutral excitation energy. Moreover, the OEP-based self-consistent DH scheme provides reasonably accurate exchange-correlation potentials and correlated densities.

Toward an Understanding of the Oxidation Process of Methionine Enkephalin: A Combined Electrochemistry, Quantum Chemistry and Quantum Chemical Topology Analysis.

Recent experimental results about the oxidation of methionine enkephalin by ·OH radicals indicated an intramolecular electron transfer between the C-terminal methionine radical cation and the tyrosine N-terminus too fast to be observed. We have investigated the thermodynamic possibility of this intramolecular electron transfer by calculating the one-electron redox potentials of both residues for several conformations of the peptide, extracted from the experimental data of the Protein Data Bank (1PLW). Using a QM/MM approach, we show that the redox potential of the Met(•+)/Met couple is higher than that of the TyrOH(•+)/TyrOH one (tyrosine is denoted as TyrOH) for all conformations. The intramolecular electron transfer between both residues (from TyrOH to Met(•+)) is thus always thermodynamically allowed. Previously, we had performed topological studies on the intramolecular electron transfer which predicted this charge transfer. A study by cyclic voltammetry pointed out that the wave belonging to methionine is not present when methionine enkephalin is oxidized and only the direct involvement of the tyrosine residue is observed.

Basis convergence of range-separated density-functional theory.

Range-separated density-functional theory (DFT) is an alternative approach to Kohn-Sham density-functional theory. The strategy of range-separated density-functional theory consists in separating the Coulomb electron-electron interaction into long-range and short-range components and treating the long-range part by an explicit many-body wave-function method and the short-range part by a density-functional approximation. Among the advantages of using many-body methods for the long-range part of the electron-electron interaction is that they are much less sensitive to the one-electron atomic basis compared to the case of the standard Coulomb interaction. Here, we provide a detailed study of the basis convergence of range-separated density-functional theory. We study the convergence of the partial-wave expansion of the long-range wave function near the electron-electron coalescence. We show that the rate of convergence is exponential with respect to the maximal angular momentum L for the long-range wave function, whereas it is polynomial for the case of the Coulomb interaction. We also study the convergence of the long-range second-order Møller-Plesset correlation energy of four systems (He, Ne, N2, and H2O) with cardinal number X of the Dunning basis sets cc - p(C)V XZ and find that the error in the correlation energy is best fitted by an exponential in X. This leads us to propose a three-point complete-basis-set extrapolation scheme for range-separated density-functional theory based on an exponential formula.

Topological analyses of time-dependent electronic structures: application to electron-transfers in methionine enkephalin.

We have studied electron transfers (ET) between electron donors and acceptors, taking as illustrative example the case of ET in methionine enkephalin. Recent pulse and gamma radiolysis experiments suggested that an ultrafast ET takes place from the C-terminal tyrosine residue to the N-terminal, oxidized, methionine residue. According to standard theoretical frameworks like the Marcus theory, ET can be decomposed into two successive steps: i) the achievement through thermal fluctuations, of a set of nuclear coordinates associated with degeneracy of the two electronic states, ii) the electron tunneling from the donor molecular orbital to the acceptor molecular orbital. Here, we focus on the analysis of the time-dependent electronic dynamics during the tunneling event. This is done by extending the approaches based on the topological analyses of stationary electronic density and of the electron localization function (ELF) to the time-dependent domain. Furthermore, we analyzed isosurfaces of the divergence of the current density, showing the paths that are followed by the tunneling electron from the donor to the acceptor. We show how these functions can be calculated with constrained density functional theory. Beyond this work, the topological tools used here can open up new opportunities for the electronic description in the time-dependent domain.

Defects and strain enhancements of second-harmonic generation in Si/Ge superlattices.

Starting from experimental findings and interface growth problems in Si/Ge superlattices, we have investigated through ab initio methods the concurrent and competitive behavior of strain and defects in the second-harmonic generation process. Interpreting the second-harmonic intensities as a function of the different nature and percentage of defects together with the strain induced at the interface between Si and Ge, we found a way to tune and enhance the second-harmonic generation response of these systems.