An official website of the United States government
We discuss an extension of the Single-Active-Electron (SAE) approximation in atoms by allowing the model potential to depend on the angular-momentum quantum number ℓ. We refer to this extension as the ℓ-SAE approximation. The main ideas behind ℓ-SAE are illustrated using the helium atom as a benchmark system. We show that introducing ℓ-dependent potentials improves the accuracy of key quantities in atomic structure computed from the Time-Independent Schrödinger Equation (TISE), including energies, oscillator strengths, and static and dynamic polarizabilities, compared to the standard SAE approach. Additionally, we demonstrate that the ℓ-SAE approximation is suitable for quantum simulations of light−atom interactions described by the Time-Dependent Schrödinger Equation (TDSE). As an illustration, we simulate High-order Harmonic Generation (HHG) and the three-sideband (3SB) version of the Reconstruction of Attosecond Beating by Interference of Two-photon Transitions (RABBITT) technique, achieving enhanced accuracy comparable to that obtained in all-electron calculations. One of the main advantages of the ℓ-SAE approach is that existing SAE codes can be easily adapted to handle ℓ-dependent potentials without any additional computational cost.
more »
« less
Single-active electron calculations of high-order harmonic generation from valence shells in atoms for quantitative comparison with TDDFT calculations
Abstract We present a reproducible ab-initio method to produce benchmark tests between time-dependent Schrödinger equation (TDSE) in the single-active-electron approximation (SAE) and time-dependent density functional theory (TDDFT) in the highly nonlinear multiphoton and tunneling regime of strong-field physics. To this end we compare results for high-order harmonic generation from valence shells in atoms using the SAE-TDSE approach and TDDFT calculations. As key to the benchmark comparison we obtain an analytic form of SAE potentials based on density functional theory, which we applied for different atoms and ions. The ionization energies of atomic ground and excited states, as well as the energies of inner shells, for the SAE potentials agree well with experimental data. Using these potentials we find remarkable agreement between the results of the two independent numerical approaches (TDDFT and SAE-TDSE) for the high-order harmonic yields in helium, demonstrating the accuracy of the SAE potentials as well as the predictive power of SAE-TDSE and TDDFT calculations for the nonperturbative and highly nonlinear strong-field process of high harmonic generation in the ultraviolet and visible wavelength regime. Finally, as another application of the SAE potentials, high harmonic spectra from outer and inner valence shells are calculated and it is shown that unphysical artifacts in the SAE-spectra from the individual shells are removed once all the amplitudes are considered.
more »
« less
- Award ID(s):
- 1734006
- PAR ID:
- 10360438
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Journal of Physics Communications
- Volume:
- 4
- Issue:
- 6
- ISSN:
- 2399-6528
- Page Range / eLocation ID:
- Article No. 065011
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We applied localized orbital scaling correction (LOSC) in Bethe–Salpeter equation (BSE) to predict accurate excitation energies for molecules. LOSC systematically eliminates the delocalization error in the density functional approximation and is capable of approximating quasiparticle (QP) energies with accuracy similar to or better than GW Green’s function approach and with much less computational cost. The QP energies from LOSC, instead of commonly used G 0 W 0 and ev GW, are directly used in BSE. We show that the BSE/LOSC approach greatly outperforms the commonly used BSE/ G 0 W 0 approach for predicting excitations with different characters. For the calculations of Truhlar–Gagliardi test set containing valence, charge transfer, and Rydberg excitations, BSE/LOSC with the Tamm–Dancoff approximation provides a comparable accuracy to time-dependent density functional theory (TDDFT) and BSE/ev GW. For the calculations of Stein CT test set and Rydberg excitations of atoms, BSE/LOSC considerably outperforms both BSE/ G 0 W 0 and TDDFT approaches with a reduced starting point dependence. BSE/LOSC is, thus, a promising and efficient approach to calculate excitation energies for molecular systems.more » « less
-
The quantum-electrodynamical time-dependent density functional theory equations are solved by time propagating the wave function on a tensor product of a Fock-space and real-space grid. Applications for molecules in cavities show the accuracy of the approach. Examples include the coupling strength and light frequency dependence of the energies, wave functions, optical absorption spectra, and Rabi splitting magnitudes in cavities, as well as a description of high harmonic generation in cavities.more » « less
-
The QMol-grid package provides a suite of routines for performing quantum-mechanical simulations in atomic and molecular systems, currently implemented in one spatial dimension. It supports ground- and excited-state calculations for the Schrödinger equation, density-functional theory, and Hartree–Fock levels of theory as well as propagators for field-free and field-driven time-dependent Schrödinger equation (TDSE) and real-time time-dependent density-functional theory (TDDFT), using symplectic-split schemes. The package is written using MATLAB’s object-oriented features and handle classes. It is designed to facilitate access to the wave function(s) (TDSE) and the Kohn–Sham orbitals (TDDFT) within MATLAB’s environment.more » « less
-
Abstract Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.more » « less
