Home »


XChem code


What is XChem? XChem package compiles a set of computational tools able to provide a full quantum mechanical description of molecular ionization in the time domain by including both electronic and nuclear degrees of freedom. The used approach has provided reliable predictions and has opened the way to new applications in emerging scientific disciplines, as attochemistry, single molecule imaging, coherent attosecond control of chemical reactions, etc.

XChem is a solution for an all-electron ab-initio calculation of the electronic continuum of molecular systems. XChem combines the tools of quantum chemistry (as implemented in Molcas) and scattering theory to accurately account for electron correlation in the single-ionization continuum of atoms, small and medium-size molecules. The validity of the XChem approach has been demonstrated, among others, in publications (2), (3) and (4) listed below.

XChem Applications:It is designed to study molecular photoionization involving multichannel scattering problem.
It is used to

  • Represent accurately molecular autoionization and Auger decay beyond the Born-Oppenheimer approximation
  • Design interfaces compatible with available ab initio Quantum Chemistry packages for a widespread use in molecular ionization problems
  • Keep computational cost at the same level as that required to evaluate bound states at the same level of theory

What can XCHEM do? XChem can compute:

  • The Close Coupling Matrix (CCM) for a user-defined set of ionization channels (each defined as an ionized molecular state coupled to electrons of given angular momenta) and including short range states relevant to the problem at hand
  • Scattering states and scattering phases by asymptotic fitting to the analytical solution
  • (Partial) photoionization cross sections within perturbation theory
  • Lifetime and character of resonances embedded in the molecular continuum, either via analysis of the cross section or via inclusion of a complex absorbing potential in the CCM yielding complex eigen-energies
  • The electron dynamics during and after ionization, caused by and probed with ultrashort laser pulses, by solving the time dependent Schrödinger Equation using the CCM
  • The angular distribution of photo electrons (MFPADs)

What can XCHEM be used for? XChem is a valuable tool for:

  • The theoretical study of ultrafast processes in pump-probe experimental method, such as Attosecond Transient Absorption Spectroscopy (ATAS) and Reconstruction of Attosecond Beatings by Interference of Two-Photon Transitions (RABBITT)
  • The investigation of photoionization of complex molecules close to threshold, where electron correlation effects may dramatically affect the (partial) photoionization cross section (Fig. 1)
  • The study of ionization processes intrinsically dependent on electron correlation, like Autoionization and Auger Decay
  • The computation of potential energy surfaces for the investigation of molecular dynamics during and after fast photoionization events

Figure 1: Photoionization cross section of N2 obtained with XCHEM. Left: total cross section compared to experiment (J .Phys. Chem. Lett. 2018, 9, 756-762). Right: Partial cross sections for different symmetries of the molecule before and after ionization and the ejected electron (Phys. Rev. A 98, 033413).

Who is XCHEM for? XChem is a valuable tool for:

  • Researchers in (computational) quantum chemistry or molecular physics interested in studying electron dynamics in the ionization continuum of molecules (e.g., photoionization, charge migration, etc)
  •  Laboratories investigating ultrafast phenomena in many-electron atoms, small and medium size molecular systems

XChem Approach and Features:XChem, at its core, lies a close coupling expansion combined with the use of

  • Gaussian
  • and B-Spline hybrid basis sets.

This approach yields the scattering states of the molecular system via the eigenstates of the close coupling matrix (CCM). While useful in their own right, the full potential lies in using the close coupling matrix as a starting point for time dependent calculations. Doing so, one may explicitly model the interaction of molecules with ultrashort (attosecond) pulses. The extreme band widths of such pulses lead to the coherent excitation of multiple ionization channels, whose coupling (accurately described in XCHEM) gives rise to complex phenomena.

Figure 2: Illustration of the XCHEM basis architecture in benzene. Cyan: B-splines, mangenta: Gaussians at the center of mass of the molecule, blue and black: Gaussians at the atomic sites not overlapping with B-splines.

An attractive feature of XCHEM is that the architecture of the basis functions (Fig. 2) and the use of Molcas allow one to describe the electronic continuum of medium-size molecules at the same level of theory as multi-configurational SCF methods (CASSCF, RASSCF) do for the ground and the lowest excited states of such molecules. At present the largest systems treated have of the order of ten atoms.


How to download


Copyright notice

The Authors of the software are Luca Argenti, Josep Vicent Borras de Llano, Inés Corral Pérez, Jesús González-Vazquez, Markus Klinker, Carlos A. Marante Valdés and Fernando Martín García, herein referred to as XChem Authors.
XChem software has been registered in the regional intellectual property register (reg. num. 16/2018/6823 and 16/2022/2643). XChem brand name and logo has been registered in Europe (request num. 018057414) and US (request num. 018057414_01).


The research carried out by the authors and giving rise to XChem-v0 software has been conducted at the Universidad Autonoma de Madrid and has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013) / ERC grant agreement n° 290853 (ERC-AdG-XCHEM) and under the European Union’s Horizon 2020 research and innovation programme, grant agreement No780284 (ERC-PoC-Imaging-XChem). XChem-v0 software has been further developed with the support of the Spanish Ministry of Science and Innotavion, National Research Agency (AEI), and the Regional Funds of the European Union (FEDER), under projects ref. ref. FIS2016-77889-R and  BES-2017-081521, giving rise to XChem-v1. Protection and analysis of XChem-v1 exploitation pathways, and furhter developments, are been carried out under project ref. PDC2021-121073-I00 funded by MCIN/ AEI /10.13039/501100011033 and by the European Union “NextGenerationEU”/PRTR .


Relevant Publications

(1) M. Nisoli et al. “Attosecond dynamics in molecules”, Chem. Rev., 117 (16), 10760, 10825, (2017), https://pubs.acs.org/doi/10.1021/acs.chemrev.6b00453
(2) C. Marante, L. Argenti and F. Martı́n, “Hybrid Gaussian–B-spline basis for the electronic continuum: Photoionization of atomic hydrogen”, Phys. Rev. A 90, 012506, (2014), https://journals.aps.org/pra/abstract/10.1103/PhysRevA.90.012506
(3) C. Marante et al. “Hybrid-Basis Close-Coupling Interface to Quantum Chemistry Packages for the Treatment of Ionization Problems”, J. Chem. Theory Comput., 13, 499−514, (2017), https://pubs.acs.org/doi/abs/10.1021/acs.jctc.6b00907
(4) C. Marante et al. “Photoionization using the XCHEM approach: Total and partial cross sections of Ne and resonance parameters above the 2s 2 2 p 5 threshold”, Phys. Rev. A 96, 022507, (2017), https://journals.aps.org/pra/abstract/10.1103/PhysRevA.96.022507
(5) M. Klinker et al. “Partial cross sections and interfering resonances in photoionization of molecular nitrogen”, Phys. Rev. A 98, 033413, (2018), https://journals.aps.org/pra/abstract/10.1103/PhysRevA.98.033413
(6) M. Klinker et al.  “Electron Correlation in the Ionization Continuum of Molecules: Photoionization of N2 in the Vicinity of the Hopfield Series of Autoionizing States” J. Phys. Chem. Lett. 2018, 9, 756−762. DOI: 10.1021/acs.jpclett.7b03220.
(7) L. Barreau, C. L. M. Petersson, M. Klinker, A. Camper, C. Marante, T. Gorman, D. Kiesewetter, L. Argenti, P. Agostini, J. Gonzalez-Vazquez, L. P. Salieres, F. DiMauro, and F. Martín. “Disentangling spectral phases of interfering autoionizing states from attosecond interferometric measurements”, Phys. Rev. Lett., 2019, 122, 253203.
DOI: https://doi.org/10.1103/PhysRevLett.122.253203
(8) B. I. Schneider, K. Bartschat, O. Zatsarinny, I. Bray, A. Scrinzi, F. Martín, M. Klinker, J. Tennyson, J. D. Gornkiel, and S. Pamidighantam, “A science gateway for atomic and molecular physics”. Proceedings of PEARC’19: Practice and Experience in Advanced Research Computing: Rise of the Machines, arXiv preprint arXiv:2001.02286 (https://arxiv.org/abs/2001.02286)
(9) B. I. Schneider, K. Bartschat, O. Zatsarinny, I. Bray, A. Scrinzi, F. Martín, M. Klinker, J. Tennyson, J. D. Gornkiel, and S. Pamidighantam, “Atomic and Molecular Scattering Applications in an Apache Airavata Science Gateway”. PEARC’20: Practice and Experience in Advanced Research Computing: Rise of the Machines, July, 2020. Pages 270-277. DOI: https://doi.org/10.1145/3311790.3397342.
(10) SM Poullain, M Klinker, J González-Vázquez, F Martín, “Resonant photoionization of O 2 up to the fourth ionization threshold”. Phys. Chem. Chem. Phys., 2019,21, 16497-16504. DOI: https://doi.org/10.1039/C9CP02150G.

XChem code

FullFun Code

M3C Code

AttoChem COST Action


XChem project

Attotrend project



Our group belongs to the Erasmus Mundus European program, which offers opportunities to obtain a master to European and non-European students. More details in Jobs.


© 2012-2023 CampuS All Rights Reserved