1
|
Ozon M, Tumashevich K, Lin JJ, Prisle NL. Inversion model for extracting chemically resolved depth profiles across liquid interfaces of various configurations from XPS data: PROPHESY. JOURNAL OF SYNCHROTRON RADIATION 2023; 30:941-961. [PMID: 37610342 PMCID: PMC10481271 DOI: 10.1107/s1600577523006124] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Accepted: 07/12/2023] [Indexed: 08/24/2023]
Abstract
PROPHESY, a technique for the reconstruction of surface-depth profiles from X-ray photoelectron spectroscopy data, is introduced. The inversion methodology is based on a Bayesian framework and primal-dual convex optimization. The acquisition model is developed for several geometries representing different sample types: plane (bulk sample), cylinder (liquid microjet) and sphere (droplet). The methodology is tested and characterized with respect to simulated data as a proof of concept. Possible limitations of the method due to uncertainty in the attenuation length of the photo-emitted electron are illustrated.
Collapse
Affiliation(s)
- Matthew Ozon
- Center for Atmospheric Research, University of Oulu, PO Box 4500, Finland
| | | | - Jack J. Lin
- Center for Atmospheric Research, University of Oulu, PO Box 4500, Finland
| | - Nønne L. Prisle
- Center for Atmospheric Research, University of Oulu, PO Box 4500, Finland
| |
Collapse
|
2
|
Pathak H, Panyala A, Peng B, Bauman NP, Mutlu E, Rehr JJ, Vila FD, Kowalski K. Real-Time Equation-of-Motion Coupled-Cluster Cumulant Green's Function Method: Heterogeneous Parallel Implementation Based on the Tensor Algebra for Many-Body Methods Infrastructure. J Chem Theory Comput 2023; 19:2248-2257. [PMID: 37096369 DOI: 10.1021/acs.jctc.3c00045] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
We report the implementation of the real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function method [ J. Chem. Phys. 2020, 152, 174113] within the Tensor Algebra for Many-body Methods (TAMM) infrastructure. TAMM is a massively parallel heterogeneous tensor library designed for utilizing forthcoming exascale computing resources. The two-body electron repulsion matrix elements are Cholesky-decomposed, and we imposed spin-explicit forms of the various operators when evaluating the tensor contractions. Unlike our previous real algebra Tensor Contraction Engine (TCE) implementation, the TAMM implementation supports fully complex algebra. The RT-EOM-CC singles (S) and doubles (D) time-dependent amplitudes are propagated using a first-order Adams-Moulton method. This new implementation shows excellent scalability tested up to 500 GPUs using the Zn-porphyrin molecule with 655 basis functions, with parallel efficiencies above 90% up to 400 GPUs. The TAMM RT-EOM-CCSD was used to study core photoemission spectra in the formaldehyde and ethyl trifluoroacetate (ESCA) molecules. Simulations of the latter involve as many as 71 occupied and 649 virtual orbitals. The relative quasiparticle ionization energies and overall spectral functions agree well with available experimental results.
Collapse
Affiliation(s)
- Himadri Pathak
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Ajay Panyala
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Nicholas P Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - John J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Fernando D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
3
|
Li J, Jin Y, Rinke P, Yang W, Golze D. Benchmark of GW Methods for Core-Level Binding Energies. J Chem Theory Comput 2022; 18:7570-7585. [PMID: 36322136 PMCID: PMC9753590 DOI: 10.1021/acs.jctc.2c00617] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The GW approximation has recently gained increasing attention as a viable method for the computation of deep core-level binding energies as measured by X-ray photoelectron spectroscopy. We present a comprehensive benchmark study of different GW methodologies (starting point optimized, partial and full eigenvalue-self-consistent, Hedin shift, and renormalized singles) for molecular inner-shell excitations. We demonstrate that all methods yield a unique solution and apply them to the CORE65 benchmark set and ethyl trifluoroacetate. Three GW schemes clearly outperform the other methods for absolute core-level energies with a mean absolute error of 0.3 eV with respect to experiment. These are partial eigenvalue self-consistency, in which the eigenvalues are only updated in the Green's function, single-shot GW calculations based on an optimized hybrid functional starting point, and a Hedin shift in the Green's function. While all methods reproduce the experimental relative binding energies well, the eigenvalue self-consistent schemes and the Hedin shift yield with mean absolute errors <0.2 eV the best results.
Collapse
Affiliation(s)
- Jiachen Li
- Department
of Chemistry, Duke University, Durham, North Carolina27708, United States
| | - Ye Jin
- Department
of Chemistry, Duke University, Durham, North Carolina27708, United States
| | - Patrick Rinke
- Department
of Applied Physics, Aalto University, Otakaari 1, FI-02150Espoo, Finland
| | - Weitao Yang
- Department
of Chemistry, Duke University, Durham, North Carolina27708, United States
| | - Dorothea Golze
- Department
of Applied Physics, Aalto University, Otakaari 1, FI-02150Espoo, Finland,Faculty
of Chemistry and Food Chemistry, Technische
Universität Dresden, 01062Dresden, Germany,
| |
Collapse
|
4
|
Plekan O, Grazioli C, Coreno M, Di Fraia M, Prince KC, Richter R, Ponzi A. Investigation of quinoline derivatives by photoemission spectroscopy and theoretical calculations. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111657] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
5
|
Vila FD, Kowalski K, Peng B, Kas JJ, Rehr JJ. Real-Time Equation-of-Motion CCSD Cumulant Green's Function. J Chem Theory Comput 2022; 18:1799-1807. [PMID: 35157796 DOI: 10.1021/acs.jctc.1c01179] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Many-body excitations in X-ray photoemission spectra have been difficult to simulate from first principles. We have recently developed a cumulant-based one-electron Green's function method using the real-time coupled-cluster-singles equation-of-motion approach (RT-EOM-CCS) that provides a general framework for treating these problems. Here we extend this approach to include double excitations in the ground-state energy and in the coupled cluster amplitudes, which have been implemented using subroutines generated by the Tensor Contraction Engine (TCE). As in the case of the singles approximation, RT-EOM-CCSD yields a nonperturbative cumulant form of the Green's function in terms of the time-dependent cluster amplitudes, adding nonlinear corrections to the traditional cumulant forms. The extended approach is applied to the core-hole spectral function for small molecular systems. We find that, when core-optimized basis sets are used, the doubles contributions reduce the mean absolute errors in the core binding energies of the 10e systems from 0.8 to 0.3 eV. They also significantly improve the quasiparticle-satellite gap by reducing its overestimation from about 3-5 to about 0-1 eV in CH4, NH3, and H2O, and also improving the overall shape of the satellite features. Finally, we demonstrate the application of the new implementation to the larger, classical XPS ESCA series of molecules and show that the singles approximation can be paired with a modest basis set to study carbon speciation.
Collapse
Affiliation(s)
- F D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - K Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - B Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - J J Kas
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - J J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
6
|
Pelimanni E, Saak CM, Michailoudi G, Prisle N, Huttula M, Patanen M. Solvent and cosolute dependence of Mg surface enrichment in submicron aerosol particles. Phys Chem Chem Phys 2022; 24:2934-2943. [PMID: 35060587 PMCID: PMC8809137 DOI: 10.1039/d1cp04953d] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Accepted: 01/14/2022] [Indexed: 01/17/2023]
Abstract
The formation of multicomponent aerosol particles from precursor solution droplets often involves segregation and surface enrichment of the different solutes, resulting in non-homogeneous particle structures and diverse morphologies. In particular, these effects can have a significant influence on the chemical composition of the particle-vapor interface. In this work, we investigate the bulk/surface partitioning of inorganic ions, Na+, Mg2 +, Ca2 +, Cl- and Br-, in atomiser-generated submicron aerosols using synchrotron radiation based X-ray photoelectron spectroscopy (XPS). Specifically, the chemical compositions of the outermost few nm thick surface layers of non-supported MgCl2/CaCl2 and NaBr/MgBr2 particles are determined. It is found that in MgCl2/CaCl2 particles, the relative abundance of the two species in the particle surface correlates well with their mixing ratio in the parent aqueous solution. In stark contrast, extreme surface enrichment of Mg2 + is observed in NaBr/MgBr2 particles formed from both aqueous and organic solution droplets, indicative of core-shell structures. Structural properties and hydration state of the particles are discussed.
Collapse
Affiliation(s)
- Eetu Pelimanni
- Nano and Molecular Systems Research Unit, Faculty of Science, University of Oulu, Box 3000, 90014, Finland.
| | - Clara-Magdalena Saak
- Department of Physics and Astronomy, Uppsala University, Box 516, 75120 Uppsala, Sweden
- University of Vienna, Department of Physical Chemistry, Währinger Straße 42, 1090 Vienna, Austria
| | - Georgia Michailoudi
- Nano and Molecular Systems Research Unit, Faculty of Science, University of Oulu, Box 3000, 90014, Finland.
| | - Nønne Prisle
- Nano and Molecular Systems Research Unit, Faculty of Science, University of Oulu, Box 3000, 90014, Finland.
- Center for Atmospheric Research, Faculty of Information Technology and Electrical Engineering, University of Oulu, P. O. Box 4500, 90014, Finland
| | - Marko Huttula
- Nano and Molecular Systems Research Unit, Faculty of Science, University of Oulu, Box 3000, 90014, Finland.
| | - Minna Patanen
- Nano and Molecular Systems Research Unit, Faculty of Science, University of Oulu, Box 3000, 90014, Finland.
| |
Collapse
|
7
|
Mejia-Rodriguez D, Kunitsa A, Aprà E, Govind N. Scalable Molecular GW Calculations: Valence and Core Spectra. J Chem Theory Comput 2021; 17:7504-7517. [PMID: 34855381 DOI: 10.1021/acs.jctc.1c00738] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a scalable implementation of the GW approximation using Gaussian atomic orbitals to study the valence and core ionization spectroscopies of molecules. The implementation of the standard spectral decomposition approach to the screened-Coulomb interaction, as well as a contour-deformation method, is described. We have implemented both of these approaches using the robust variational fitting approximation to the four-center electron repulsion integrals. We have utilized the MINRES solver with the contour-deformation approach to reduce the computational scaling by 1 order of magnitude. A complex heuristic in the quasiparticle equation solver further allows a speed-up of the computation of core and semicore ionization energies. Benchmark tests using the GW100 and CORE65 data sets and the carbon 1s binding energy of the well-studied ethyl trifluoroacetate, or ESCA molecule, were performed to validate the accuracy of our implementation. We also demonstrate and discuss the parallel performance and computational scaling of our implementation using a range of water clusters of increasing size.
Collapse
Affiliation(s)
- Daniel Mejia-Rodriguez
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Alexander Kunitsa
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Edoardo Aprà
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|