Armoured: Atomic R-matrix Method For Relativistic Dynamics

Lead Research Organisation: Queen's University Belfast
Department Name: Sch of Mathematics and Physics

Abstract

In order to 'view' electrons in the atom, we need to be able to describe their motion on a time-scale comparable to the interactions themselves. This is akin to taking a photograph- the faster an object is moving, the shorter the exposure time required to capture it. The process of harmonic generation- where an electron in an atom is driven by a laser field in such a way as it re-emits high harmonics of the laser light- has provided a unique tool for probing atomic and molecular structure- allowing the imaging of molecular orbitals and even videos of chemical reactions taking place. While this field of attosecond physics (1 attosecond = 1 billionth of a billionth of a second) is well established, the theoretical description and computational models of the underlying processes are underdeveloped.

In order to fully describe the fundamental dynamics of laser driven electrons in the atom we need to consider not just the effect of the laser field, nor even simply the net effect of the electrons, we must also describe the complex interactions between the many electrons. In order to facilitate this description in a computationally tractable way we have developed time-dependent R-matrix theory (TDRM), and associated computer codes. The theory makes use of an R-matrix division of configuration space in order to simplify the calculations without neglecting important multielectron effects.

It is our intention to extend the TDRM technique to describe relativistic effects in ultrafast processes. The spin-orbit interaction, wherein an electron's intrinsic (spin) and orbital angular momenta interfere to change its behaviour, gives rise to fine structure splittings in atomic energy levels. The transitions of electrons between these levels occurs on a time-scale governed by the energy gap. For heavier elements such as krypton and xenon this time-scale is of the order of a few femtoseconds- and so dynamics can evolve within a laser pulse. By opening up new electron-emission channels, and changing the selection rules for electron transitions, the spin-orbit effect can significantly alter the fundamental processes of attosecond physics- ionisation, harmonic generation etc. Hence, interesting new dynamics can be observed on the atomic scale. By using the TDRM method we will be able to model these dynamics from first principles, thus giving an unparallelled insight into some of the most fundamental physical processes in atomic science.

Planned Impact

FUNDAMENTAL SCIENCE

Understanding of fundamental phenomena

The calculations we perform, and the scientific problems we will address will directly inform our understanding of the very fundamental fields of ultrafast electron dynamics, electron correlation effects and atomic structure.

Improved understanding in these areas will directly inform:

- Astrophysics: atomic opacities, scattering cross sections (CS)
- Materials science: electron transport
- Fusion energy: Scattering CS
- Laser physics: electron dynamics, ultrafast phenomena, high harmonic generation
- Radiation damage: electron dynamics in molecular systems
- Plasma physics: Scattering CS

Improve and develop experimental techniques

We will directly address high harmonic spectroscopy and attosecond transient absorption spectroscopy. Performing HHG calculations, and thereby uncovering possible improvements in efficiency or generality of the process impacts on

- XUV-initiated HHG
- Attosecond laser pulse generation
- Multidimensional spectroscopy
- Free-electron laser science

The proposed improvements to the computational approaches we propose will impact on ongoing research into fusion energy.

COMPUTATIONAL PHYSICS

Understand, improve and enhance computational methods

The specific area of computational physics we address will be ab-initio modelling of multi-electron systems. However, the techniques we use:

- novel mixed basis-set/finite difference techniques
- mixed (shared/distributed) memory parallelism
- novel hardware for HPC (deployment of RMT on Intel Xeon Phi architecture)

are more generally applicable in a variety of areas using HPC.

ECONOMIC AND SOCIETAL

Improved understanding of light mediated charge transfer, allows better understanding and control of

- vision
- photosynthesis
- nano-scale electronics
- high-speed electronics
- the manufacture and efficient operation of photovoltaic cells

with clear implications for medicine, human and plant biology, emerging technologies and renewable energy respectively.

Direct engagement with two ongoing projects in magnetic confinement fusion has important long term consequences for clean energy, and knock-on impacts for defence and security.

The high performance computing aspects of the work will inform the development and design of new HPC facilities. The testing of novel HPC technologies (at present the Intel Xeon Phi coprocessor) forms an important part of our engagement with the Irish Centre for High End Computing. This has far reaching implications for all the research areas which use HPC worldwide including

- meteorology
- medicine
- aerospace engineering
- renewable energy research
- defence

We will directly engage the public to improve understanding and awareness of fundamental science. In particular, we will create video introductions to the areas of

- ultrafast physics
- supercomputing for science

By directly engaging with three companies we will establish new connections and develop existing links between fundamental science and industry in the areas of

- plasma physics (via Quantemol)
- laser physics and metrology (via M2)
- computer component design (via Intel)

Improving the understanding of the SO interaction, will inform the improvement and potentially, the development of new industrial techniques. In particular the SO interaction is key in

- magnetocrystaline anisotropy, the key process in the manufacture of ferromagnetic materials
- giant magnetoresistance, the key process in the operation of hard-drive read-heads

Both the PDRA and PhD student will receive high-end technical training relevant to academia and industry, both in HPC techniques and fundamental science. Our record in preparing researchers for positions outside of academia is evidenced by the current high-level positions held by former CTAMOP members: Michael Lysaght at the Irish Centre for High-End Computing and Laura Moore at SAP.

Publications

10 25 50
 
Description We have been successful in the development of computational codes to enable the description of ultra-fast processes in heavy atoms. The Breit-Pauli atomic R-matrix codes have been redesigned to enable us to connect the output of the Breit-Pauli R-matrix codes to the current time-dependent R-matrix RMT codes. This required the implementation of new basis-set techniques within the R-matrix I style of R-matrix codes. New interface routines were developed so that the output from both R-matrix I and R-matrix II codes can be used within the time-dependent codes.

The data sets generated using the Breit-Pauli R-matrix approach have been used within the RMT codes to investigate ultra-fast dynamics induced by the relativistic spin-orbit interaction. In a test case for the codes, we demonstrated that we can select the excitation of specific autoionising states through selection of the polarisation direction of the probe pulse in an attosecond pump-probe experiment through exploiting the spin-orbit interaction.
This first successful application of the codes demonstrates that they are capable of elucidating ultra-fast dynamics in heavy multi-electron systems. This is an area of strong, topical experimental interest, but which is very difficult to study reliably with theory. We have further used these codes to demonstrate that ultra-short light pulses can aid the spin-polarisation of ejected electrons through exploitation of the different preferred emission angles for co- and counter-rotating electrons. The codes developed within this project therefore provide a unique scientific capability to provide theoretical support to state-of-the-art experimental interests.

The codes developed as part of this project have been published to make them available to the scientific community, with a full description of the codes and their operation. In addition, codes are made available through code repositories. These repositories are maintained to ensure that the latest updates to the codes are also accessible to the community.
Exploitation Route The development of the codes will allow greater accuracy in the description of plasma models. These models are of interest for the industrial use of plasmas, as well as all scientific research relating to plasmas, eg. those created in intense light fields, or in astrophysics.

Several current topics of experimental interest require information that can only be provided through codes with the capability of the new time-dependent codes. It is our intention to apply the current codes to a number of problems relevant to experiment to demonstrate the capability of the codes. These studies will act as a showcase on the type of calculations that the codes can perform, and hopefully lead to new joint experimental and theoretical studies.

One possible application area of the current codes is the production of spin-polarised electrons. Circularly polarised laser fields can be tuned such that the electron emission is dominated by electrons that move in the opposite direction around the atom. When this occurs in heavy atoms, one can preferentially generate electrons with a particular spin direction. Studies using the present code can identify schemes that will optimise the overall efficiency of the process.

It is difficult to forecast how progress in experiment will occur, and what innovations this will lead to, but the current codes have a unique capability that will facilitate progress in a range of areas.
Sectors Digital/Communication/Information Technologies (including Software),Other

 
Description This project modified the R-matrix I suite of codes to make the approach compatible with the R-matrix with time dependence codes. It also opened up a new area of ultra-fast research by allowing the study of dynamics driven by spin-orbit coupling. This may benefit future research involving heavy atoms, such as iodine, which can be used as a marker in various (bio)chemical studies. The code is made available to the community via software repositories, and as part of the atomic and molecular software gateway, AMOSgateway.org. This portal provides access to software developed by research groups around the world for use in research and education world-wide. The gateway has greatly simplified the training for young researchers to become familiar with the wide range of codes available to the community. High-quality training can now be provided at any location in the world without the need to set up extensive training suites locally. The success of the gateway is based on the availability of a good range of codes for different physics purposes, and the work within this project has helped to broaden the range.
First Year Of Impact 2022
Sector Digital/Communication/Information Technologies (including Software),Education
Impact Types Policy & public services

 
Description AQuA DIP: Advanced Quantum Approaches to Double Ionisation Processes
Amount £865,857 (GBP)
Funding ID EP/T019530/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 05/2020 
End 10/2024
 
Description PARAMOR- Platform And Resource for Atomic, Molecular and Optical Research
Amount £668,309 (GBP)
Funding ID EP/V05208X/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 07/2021 
End 06/2026
 
Description RECAP: R-matrix calculations for Electron Correlation in Attosecond Processes
Amount £33,547 (GBP)
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 01/2018 
End 12/2018
 
Description UK Atomic, Molecular and Optical physics R-matrix consortium (UK AMOR)
Amount £368,071 (GBP)
Funding ID EP/R029342/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 07/2018 
End 07/2023
 
Title Cooper Minimum in PICS and HHG spectra in Ar+ 
Description Data for harmonic generation spectra and photoionisation cross section from neutral and singly ionized argon/ 
Type Of Material Database/Collection of data 
Year Produced 2018 
Provided To Others? Yes  
Impact none yet 
URL https://pure.qub.ac.uk/portal/en/datasets/cooper-minimum-in-pics-and-hhg-spectra-in-ar(df98f631-13ce...
 
Title Data for Resolving Ultra-Fast Spin-Orbit Dynamics in Heavy Many-Electron Atoms 
Description Data associated with the manuscript Resolving "Ultra-Fast Spin-Orbit Dynamics in Heavy Many-Electron Atoms" 
Type Of Material Database/Collection of data 
Year Produced 2019 
Provided To Others? Yes  
Impact Proposed Experimental scheme for resolving ultrafast spin-orbit dynamics in autoionisation 
URL https://pure.qub.ac.uk/portal/en/datasets/data-for-resolving-ultrafast-spinorbit-dynamics-in-heavy-m...
 
Title Dataset for "Enhancing spin polarization using ultrafast angular streaking" 
Description This dataset contains result pertaining to a forthcoming paper "Enhancing spin polarization using ultrafast angular streaking" to be published in Physical Review A. Data included are the text files used to generate the figures, as well as python plot scripts. 
Type Of Material Database/Collection of data 
Year Produced 2021 
Provided To Others? Yes  
Impact First demonstration of angular streaking technique to isolate spin-polarised electrons 
URL https://pure.qub.ac.uk/en/datasets/dataset-for-enhancing-spin-polarization-using-ultrafast-angular-s...
 
Title Electron correlation and short-range dynamics in attosecond angular streaking in the Fluorine anion. 
Description This dataset contains results from an article "Electron correlation and short-range dynamics in attosecond angular streaking" to be published in Physical Review A. Data included are the input files needed to run the calculations, and output text files and plot scripts used to generate the figures. In particular the data represent the radial electronic wave function of the Fluorine anion after irradiation by an 800nm laser pulse. These data may then be constructed into angularly resolved spectra using the plot scripts provided. As explained in the associated article, the results demonstrate a negative angular shift in the main photoelectron peak. This indicates that the angular shifts found in attosecond angular streaking experiments may not be attributed directly to the Coulombic interaction with the residual ion, because, in this particular instance, the residual 'ion' is neutral. 
Type Of Material Database/Collection of data 
Year Produced 2020 
Provided To Others? Yes  
Impact Indicates that the angular shifts found in attosecond angular streaking experiments may not be attributed directly to the Coulombic interaction with the residual ion, because, in this particular instance, the residual 'ion' is neutral. 
URL https://pure.qub.ac.uk/en/datasets/electron-correlation-and-shortrange-dynamics-in-attosecond-angula...
 
Title Electron vortices in He and F- 
Description These data are raw wave function data for He and F- interacting with arbitrarily polarised light fields. The data are produced by the R-matrix with time-dependence code. Instructions for post-processing can be found in the README file. 
Type Of Material Database/Collection of data 
Year Produced 2018 
Provided To Others? Yes  
Impact Several papers on electron vortices followed 
URL https://pure.qub.ac.uk/portal/en/datasets/electron-vortices-in-he-and-f(668e9c0f-a1f1-4d48-b982-56b3...
 
Title Extreme-ultraviolet-initiated High-harmonic Generation in Ar+ 
Description Data relevant to paper: Extreme-ultraviolet-initiated High-harmonic Generation in Ar+ published in physical review A (2018) This is the raw RMT data for calculations of the harmonic spectra, generated by Ar+, in both single (IR) and two-colour (XUV + IR) laser fields. The data are contained in a set of directories, named in accordance with their corresponding figures in the article. Each file, in turn, is named consistently with the label for the spectrum appearing in the appropriate figure (e.g., in the legend). The only exception to this pertains to Figure 8(a), for which the individual data files have been named according to the time delay between the XUV and IR pulses. Note that the latter are specified in units of the XUV pulse period. For precise details regarding the choice of pulse parameters, please refer to the main text and figure captions of the article. 
Type Of Material Database/Collection of data 
Year Produced 2018 
Provided To Others? Yes  
Impact None yet 
URL https://pure.qub.ac.uk/portal/en/datasets/extremeultravioletinitiated-highharmonic-generation-in-ar(...
 
Title Modeling tomographic measurements of photoelectron vortices in counter-rotating circularly polarized laser pulses 
Description This dataset contains result pertaining to a forthcoming paper "Modeling tomographic measurements of photoelectron vortices in counter-rotating circularly polarized laser pulses" to be published in Physical Review A. Data included are the input files, output files, and text files used to generate the figures. 
Type Of Material Database/Collection of data 
Year Produced 2019 
Provided To Others? Yes  
Impact First full ab-initio theoretical modelling of tomography in electron vortices for general. multielectron systems 
URL https://pure.qub.ac.uk/en/datasets/modeling-tomographic-measurements-of-photoelectron-vortices-in-co...
 
Title Rotational asymmetry in photodetachment from F- in circularly polarised laser fields 
Description These data are raw wave functions for F- following photo detachment by circularly polarised laser fields. The data are produced using the R-matrix with time-dependence code. Information on the datasets and their post processing may be found the the README file and channel data file. 
Type Of Material Database/Collection of data 
Year Produced 2019 
Provided To Others? Yes  
Impact None yet 
URL https://pure.qub.ac.uk/portal/en/datasets/rotational-asymmetry-in-photodetachment-from-f-in-circular...
 
Title Interface between Breit-Pauli R-matrix codes and time-dependent R-matrix code RMT 
Description The software allows relativistic atomic data to be used in time-dependent calculations on ultra-fast processes in atoms 
Type Of Technology Software 
Year Produced 2018 
Impact It is too early to discuss impacts of this piece of software, which completed development in 2018. 
 
Title RMT: R-matrix with time-dependence 
Description RMT is a program which solves the time-dependent Schrödinger equation for general, multielectron atoms, ions and molecules interacting with laser light. As such it can be used to model ionization (single-photon, multiphoton and strong-field), recollision (high-harmonic generation, strong-field rescattering) and, more generally, absorption or scattering processes with a full account of the multielectron correlation effects in a time-dependent manner. Calculations can be performed for targets interacting with ultrashort, intense laser pulses of long wavelength and arbitrary polarization. Calculations for atoms can optionally include the Breit-Pauli correction terms for the description of relativistic (in particular, spin-orbit) effects. 
Type Of Technology Software 
Year Produced 2019 
Open Source License? Yes  
Impact Several papers and data sets have been produced using calculations from RMT 
URL https://www.sciencedirect.com/science/article/abs/pii/S0010465519303856
 
Title Time-dependent R-matrix code for the description of atomic and molecular systems in arbitrarily polarised light including semi-relativistic systems 
Description A highly parallelised scientific code that enables the description of atomic and molecular systems in ultra-short light fields of arbitrary polarisation. This code has the capability to describe semi-relativistic systems accounting for the spin-orbit interaction. 
Type Of Technology Software 
Year Produced 2019 
Open Source License? Yes  
Impact This code has enabled scientific studies on the impact of spin-orbit dynamics on the atomic response to intense laser light. It has also enabled studies on ultrafast dynamics in arbitrary polarised light. The codes have also been applied to study ultrafast molecular dynamics. 
 
Title Time-independent R-matrix codes 
Description The project has enabled us to update the already existing R-matrix codes for photoionization and scattering, The improved efficiency of the codes has enabled us to carry out calculations for atomic systems which require larger basis sets to obtain reliable outcomes. 
Type Of Technology Software 
Year Produced 2018 
Open Source License? Yes  
Impact The software has enabled us to carry out calculations for atomic ions in the iron sequence, with improved accuracy for inclusion in plasma models. It also provides input data for the time-dependent R-matrix codes to enable the study of ultra-fast processes of heavy atoms in light fields driving by spin-orbit dynamics. 
 
Description Attochem conference 2020 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Postgraduate students
Results and Impact International conference on Attosecond dynamics in atoms and molecules. Dr. Andrew Brown gave a talk entitled "The RMT project" during which he reported on recent results obtained and tools developed under the auspices of the UK-AMOR consortium. The major impact was the arrangement of several new collaborative studies with audience members.
Year(s) Of Engagement Activity 2020
URL https://atom.ubbcluj.ro/attochem/#attochem-registration
 
Description Industry sponsored workshop for code users 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Postgraduate students
Results and Impact 32 researchers attended a two day workshop at UCL. The workshop was co-sponsored and organised by Quantemol Ltd. an industrial partner. The workshop was divided into several hands-on sessions where the computer codes developed as part of the grant were demonstrated. One session was devoted to training in version control software tools provided by a representative from the ARCHER supercomputing service. One session was used to discuss the scientific impact of the work.
Year(s) Of Engagement Activity 2019
URL https://rmadam2019.wordpress.com/home/programme/
 
Description UK-AMOR kick-off meeting 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Postgraduate students
Results and Impact This was a collaborative workshop involving many research groups within the UK, who are working with the R-matrix codes. An overview of the development work being undertaken was given, alongside a presentation of research enabled by these developments. Groups who use the R-matrix codes attended the meeting so that feedback on these developments could be obtained. The meeting also provided an overview of the current set-up for distribution of the codes to interested end-users.
Year(s) Of Engagement Activity 2018