Subscriber access provided by Kansas State University Libraries

Article

Electron Scattering from Pyridine Agnieszka Sieradzka, Francisco Blanco, Martina C. Fuss, Zdenek Masin, Jimena Gorfinkiel, and Gustavo Garcia J. Phys. Chem. A, Just Accepted Manuscript • Publication Date (Web): 09 Jun 2014 Downloaded from http://pubs.acs.org on June 16, 2014

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Electron Scattering from Pyridine Agnieszka Sieradzka,† F. Blanco,‡ M. C. Fuss,¶ Zdenˇek Maˇs´ın,† J. D. Gorfinkiel,∗,† and G. Garc´ıa¶,§ Department of Physical Sciences, The Open University, Walton Hall, Milton Keynes, MK7 6AA, United Kingdom, Departamento de F´ısica At´omica, Molecular y Nuclear, Universidad Complutense de Madrid, 28040 Madrid, Spain, Instituto de F´ısica Fundamental, Consejo Superior de Investigaciones Cient´ıficas, Serrano 113-bis, 28006 Madrid, Spain, and Centre for Medical Radiation Physics, University of Wollongong, 2522 NSW, Australia E-mail: [email protected].

Abstract We have calculated cross sections for elastic and inelastic electron scattering from pyridine in the energy range 1 eV to 1 keV. The R-matrix and IAM-SCAR methods have been used for low and higher collision energies respectively. Agreement with available theoretical data is good. We have also looked at the formation of shape resonances and compared our results with existing experimental data. We compare the results with data for electron scattering with pyrimidine. ∗

To whom correspondence should be addressed Department of Physical Sciences, The Open University, Walton Hall, Milton Keynes, MK7 6AA, United Kingdom ‡ Departamento de F´ısica At´ omica, Molecular y Nuclear, Universidad Complutense de Madrid, 28040 Madrid, Spain ¶ Instituto de F´ısica Fundamental, Consejo Superior de Investigaciones Cient´ıficas, Serrano 113-bis, 28006 Madrid, Spain § Centre for Medical Radiation Physics, University of Wollongong, 2522 NSW, Australia †

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Keywords 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

electron-molecule collisions, cross sections, resonances, organic compounds

Introduction Understanding and characterizing electron scattering from molecules is crucial in a number of applied fields (e.g. plasma processes in both technological and astrophysical environments, radiation damage in biological matter, etc.). Experimental and theoretical/computational techniques continue to be developed and improved in order to provide more accurate quantitative information and detailed insight into the molecular processes induced by electron collisions. 1 Increased interest in understanding, at a microscopic level, the effect of radiation on biological systems has led to the study of electron scattering from biological molecules: 2 nucleobases, deoxyribose, amino acids, etc.. Much of this work has focused on the dissociative electron attachment (DEA) process known to lead to DNA-strand breaks. 3,4 Low energy experimental work has looked at the anions produced in DEA from DNA components. Theoretical work has been mostly restricted to studying resonance (or temporary negative ion) formation, with a few studies investigating the dissociation step in one dimension (e.g. Gallup and Fabrikant 5 ). At the same time, the accurate determination of cross sections for all the possible processes (elastic scattering, target excitation, DEA, ionization, etc.) for a broad energy range is also receiving increased attention. This is due to the need for these data in single-track structure simulations used for modelling radiation-induced damage in biological matter. 6 These simulations require quantitative data for all the possible processes that take place over the whole range of available kinetic energies. Key among these processes are electron collisions, as large quantities of electrons are produced by the primary radiation. 7 Many recent studies, particularly theoretical ones, have focused on model molecules: for 2

ACS Paragon Plus Environment

Page 2 of 26

Page 3 of 26

The Journal of Physical Chemistry

example, tetrahydrofuran was initially studied instead of deoxyribose 8–12 because it has fewer 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

electrons (particularly helpful for ab initio calculations) and higher symmetry. Similarly, pyrimidine (C4 H4 N2 ) and its isomers (pyrazine and pyridazine) have attracted significant attention as models for the pyrimidinic nucleobases: thymine, cytosine, and uracil (see Sanz et al. 13 and references therein). These studies provide useful information that can be extrapolated to the scattering processes with the molecules of interest. 14 Differences in the collision processes and cross sections for the diazines are mainly due to their different dipole moment (other properties, like polarizability and ionization energy are very similar for the three of them); low energy resonances are very similar and can be traced back to benzene. It seemed therefore of interest to study pyridine, C5 H5 N, that possesses similar physicochemical properties to the diazines, in order to ascertain how the single (instead of double) substitution of a carbon atom in the benzene ring affects the collision processes and its theoretical description. Low energy electron scattering from pyridine has been studied computationally 15 using the Schwinger multichannel method. The authors reported elastic cross sections as well as the positions of the low-lying π ∗ resonances. Modelli and Burrow 16 and Nenner and Schulz 17 performed electron transmission spectroscopy experiments (ETS) on substituted pyridines and azabenzenes, among them pyridine, respectively. This allowed them to characterize three low-lying (below 5 eV) π ∗ resonances in all these targets.

The R-matrix method The R-matrix method 18 and its application to electron-molecule scattering 19 have recently been described in detail. We will nonetheless briefly summarize the method here. The R-matrix calculations reported in this work have been performed within the fixed-nuclei approximation. The key idea of the R-matrix method is the division of configuration space into two

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Page 4 of 26

regions, an inner and an outer region. These are separated by a sphere of radius r = a 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

centred on the centre of mass of the system. The R-matrix sphere must be large enough to contain the charge density of the target states of interest and that associated to the target orbitals used to build the L2 functions (see below). The more complex inner region problem is solved first. In this region, the wave function for the (N + 1) electron system can be written as a linear combination:

ΨΓk (X1 ....XN+1 )

= A +

nc n X X

Φi (X1 ....XN ; ˆ rN+1 ; σN +1 ) ×

i=1 j=1 m X

χi (X1 ....XN+1 )bik ,

uij (rN +1 ) aijk rN +1 (1)

i=1

where Xi represent the space and spin coordinates of electron i and σN+1 indicates the spin coordinates of the (N + 1)th electron. The operator A ensures the antisymmetrization of the whole wavefunction. The functions Φi (XN ; ˆ rN+1 ; σN +1 ) are built as products of the target wavefunctions Φi of each of the n states included in the calculation and the angular (spherical harmonics) and spin functions of the scattering electron. The functions uij describe the radial behaviour of the scattering electron and are generated as linear combinations of Gaussian functions centred on the centre of mass of the molecule. The L2 -integrable functions χi are built as products of target orbitals and they are crucial for the representation of shortrange correlation and polarization effects. The spin-space symmetry of the wavefunctions is denoted by Γ. In order to obtain the coefficients aijk and bik , the matrix obtained from the sum of the (N + 1) Hamiltonian and the Bloch operator 18 is diagonalized. This generates a set of ΨΓk and their associated eigenvalues, Ek , that allow us to describe the (N + 1) system in the inner region. In the outer region, where exchange between the scattering and target electrons can be neglected, a single centre expansion of the electron-molecule interaction potential is used. 4

ACS Paragon Plus Environment

Page 5 of 26

The Journal of Physical Chemistry

The radial functions describing the behaviour of the scattering electron are determined by 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

solving a set of coupled differential equations. This is done by propagating the R-matrix, constructed using inner region (ΨΓk and Ek ) and target (to define the channels) information. This propagation is carried out to a radius large enough so that an asymptotic expansion for the radial wavefunctions of the scattering electron in each channel can be used. K-matrices are determined by matching these to known analytical asymptotic solutions. Integral and differential cross sections, resonances’ properties, etc. can be obtained from these. In our calculations we have used the UKRmol suite 20 and, in order to calculate elastic differential cross sections (DCS), POLYDCS. 21 For polar targets, the fixed-nuclei cross sections diverge due to the long-range nature of the electron-dipole interaction. POLYDCS, implemented by Sanna and Gianturco, uses a frame-transformation method (based on the multipole-extracted adiabatic-nuclei procedure proposed by Norcross and Padial 22 ) to overcome this problem (see Gianturco and Jain 23 for more details). This is one of several first Born approximation based ’top up’ procedures: those partial waves not included in the ab initio calculation are incorporated in an approximate way. (Another way of seeing this approach is by interpreting the ab initio contribution as a correction to the Born cross section.)

The IAM-SCAR method The IAM-SCAR method is based on an independent atom representation (IAM) complemented with a screening corrected additivity rule (SCAR) and it has already been extensively employed to calculate electron scattering cross sections for a wide variety of molecular targets 24–26 over a broad energy range (0.1-1000 eV). We only briefly reiterate the key points of this approach here. Initially, this approximation does not consider the molecule as a single target but as an aggregate of atoms which scatter independently assuming that molecular binding does not affect their electronic distribution. The first subjects of these calculations

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Page 6 of 26

are therefore the constituent atoms, namely H, C and N. Each atomic target is represented 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

by an interacting complex (optical) potential, Vopt (r), whose real part accounts for the elastic scattering of the incident electrons while the imaginary part represents the inelastic processes that are considered as “absorption” from the incident beam. This optical potential can be expressed as: Vopt (r) = Vs (r) + Vex (r) + Vpol (r) + iVabs (r)

(2)

where Vs (r) is the static term derived from the Hartree-Fock calculation of the atomic charge density. 27 Vex (r) is the exchange term which accounts for the indistinguishability between the incident and target electrons. The expression chosen for this term is the semi-classical energy-dependent formula derived by Riley and Truhlar. 28 Vpol (r) is the polarization term which describes the long-range interactions and depends on the target dipole polarizibility, in the form given by Zhang et al. 29 Finally, the absorption potential Vabs (r), which accounts for the inelastic processes, is based on Staszewska’s quasi-free model. 30 Initially some divergences were found when results were compared to the available atomic scattering data. After including some improvements such as many-body and relativistic corrections, screening effects inside the atom, local velocity correction and in the description of the electrons’ indistinguishability, the model proved to provide a good approximation for electron-atom scattering 31,32 over a broad energy range. An excellent example of this is elastic electron– atomic iodine (I) scattering, 33 where the optical potential results compare very favourably with those from a sophisticated Dirac-B-spline R-matrix computation. Within this model we numerically integrate the radial scattering equation and obtain the complex partial wave phase shifts δl . 34 Using these phase shifts, in combination with the optical theorem, we can generate the atomic scattering amplitudes: lmax 1 X (2l + 1)(e2iδl + 1)Pl (cos θ) f (θ) = 2ik l=0

(3)

where θ is the scattering angle and k the momentum of the incident electron. From these we

6

ACS Paragon Plus Environment

Page 7 of 26

The Journal of Physical Chemistry

can determine differential (dσel /dΩ) and integral (σel ) elastic cross sections, as well as the 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

total scattering cross sections using standard scattering expressions. In order to obtain molecular cross sections, the IAM has been followed by applying a coherent addition procedure, commonly known as the additivity rule (AR). In this approach, the molecular scattering amplitude, F (θ), is derived from the sum of the above atomic molec amplitudes which lead to the differential elastic cross section for the molecule (dσel /dΩ),

according to

F (θ) =

X

fi (θ)eiq·ri

atoms

molec X dσel sin qrij = fi (θ)fj∗ (θ) dΩ qrij i,j

(4)

where q is the momentum transferred in the scattering process and rij is the distance between atoms i and j. molec Integral elastic cross sections for the molecule can be determined by integrating (dσel /dΩ).

Alternatively, elastic cross sections can be derived from the atomic scattering amplitudes in conjunction with the optical theorem 31 giving

molec σel =

X

atoms σel

(5)

atoms

Unfortunately, in its original form, we found an inherent contradiction between the integral cross section derived from those two approaches, which suggested that the optical theorem was being violated. 29 The main limitation of the AR is that no molecular structure is considered, thus it is really only applicable when the incident electrons are fast enough to effectively “see” the target molecule as a sum of the individual atoms (typically above ∼100 eV). To reduce this limitation we developed the SCAR method 32,35 which considers the geometry of the corresponding molecule (atomic positions and bond lengths) by introducing some screening coefficients which modify both differential and integral cross sections, especially for lower energies. With this correction the range of validity of the IAM-SCAR method might be

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Page 8 of 26

extended down to about 10 eV, depending on the target. For intermediate and high energies 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(20-5000 eV), this method has been proved to be a powerful tool to calculate electron scattering cross sections from a high variety of molecules of very different sizes, from diatomic to complex biomolecules. 36 From the above description of the IAM-SCAR procedure it is obvious that vibrational and rotational excitations are not considered in this calculation. However, for polar molecules such as pyridine, additional dipole-induced excitation cross sections can be calculated following the procedure suggested by Jain: 37 differential and integral rotational excitation cross sections for a free electric dipole are calculated in the framework of the first Born approximation (FBA). These are incorporated to our IAM-SCAR calculation in an incoherent way, adding the results as an independent channel. Although rotational excitation energies are, in general, very low (typically a few meV) in comparison with the incident electron energies, in order to validate the Born approximation the latter energies should be higher than about 10 eV. Under these circumstances, rotational excitation cross sections J → J’ are calculated by weighting the population for the Jth rotational quantum number at 300 K and estimating the average excitation energy from the corresponding rotational constants. We call the whole procedure, that it has been successfully used for other polar molecules such as H2 O, HCN, and pyrimidine, 38–40 the IAM-SCAR+Rotations method. Additionally, when the permanent dipole moment of the molecule is relatively large, as is the case of pyridine, the FBA also fails for medium and large scattering angles. In order to partially solve this situation, we introduced a correction based on that suggested by Dickinson, 41 which brings a substantial improvement for electron scattering cross sections with strongly polar molecules. This procedure introduces a first-order corrective term  Dck  dσ to the differential cross sections for medium and large angles but maintaining the dΩ  B FBA correction dσ for lower angles: dΩ dσ B µ2 1 θ ≤ θc ≃ 2 dΩ 6Ei sin (θ/2)

dσ Dck πµ 1 θ ≥ θc ≃ 3 dΩ 64Ei sin (θ/2) 8

ACS Paragon Plus Environment

(6)

Page 9 of 26

The Journal of Physical Chemistry

Here µ is the permanent dipole moment of the molecule and Ei the energy of the projectile. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Providing that the dipole moment is bigger than µ =0.75 Debye, both curves smoothly join together at θc , the critical angle at which they cross each other.

Details of the calculations Pyridine, C5 H5 N, has a vertical ionization energy 42 of around 9.51 eV, a dipole moment of 2.19 Debye 43 and its experimentally determined spherical polarizability 42 is 64 a30 . It has 42 electrons and belongs to the C2v point group. Pyridine is an asymmetric top, but two of its rotational constants differ by less than 5%, so it can be approximated as a symmetric top. Its first electronic excitation threshold is around 4 eV. 44 The calculations presented here are based on the experience gained in the study of electron collisions with diazines, both at low 45,46 and intermediate and high energies. 13,47

R-matrix calculations We used the ground state equilibrium geometry optimized at the QCISD level 42 (other equilibrium geometries were tested with no difference in the scattering results). Following our earlier work, 46 we performed calculations with two different basis sets: a compact, ccpVDZ, and a diffuse one, 6-311+G**. Hartree-Fock SCF (HF) orbitals were generated using MOLPRO 48 and used in the scattering calculation. The calculated ground state energy and dipole moment are listed in Table 1. The agreement of the HF energy with the accurate value for the energies is, as expected, not very good. The dipole moment (very important in the scattering description) agrees much better with the experimental value. We used an R-matrix radius a=13 a0 for calculations with the compact basis set and a=18 a0 for the diffuse one and the appropriate continuum basis sets. 49,50 For the smaller R-matrix radius, we tested the effect of including an extra partial wave (l = 5) in the calculations. We performed the calculations at the Static-Exchange plus Polarization (SEP)

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 26

Table 1: Energy (in Hartree) and dipole moment (in Debye) of the ground state of pyridine calculated at the HF level, using the basis sets indicated in the table. Also listed are the experimental value of the ground state dipole moment 43 and the calculated energy at B3LYP level using the 6-311+G(3df,2p) basis set. 42 E µ

cc-pVDZ 6-311+G** Acc. value −246.71 −246.75 -248.37 2.22 2.43 2.19

level: only the ground state of the target, described at HF level, is included in expression (1). In order to describe short range correlation-polarization effects, we include in the set of L2 functions (χi in expression (1)) configurations of this type: (core)Nd (valence)N −Nd (virtual)1

(7)

with N=42 and Nd =34. Here, ’valence’ denotes all those orbitals occupied in the ground state configuration from which single excitations are allowed, ’virtual’ denotes a subset of those HF orbitals that are unoccupied in the ground state configuration. Unfortunately, it is not possible to determine a priori the optimal number of virtual orbitals to include; if too many are used the calculations can be ’overcorrelated’. 46 Our approach for the diazines and uracil, that we have followed here, is to start with a small number of the lowest energy virtual orbitals and increase this number (always adding orbitals in increasing energy order) until we achieve good agreement with the experimental positions of the π ∗ resonances; for pyridine, we have used the resonance positions of Nenner and Schulz. We have determined as optimal the use of 35 virtual orbitals in the calculations with the compact basis set and 50 when the diffuse basis set is used. These are consistent with those required to obtain the best agreement with experiment in the case of the diazines: 35 and 40 virtuals respectively were enough in that case (although published data 46 used 25 and 35 we have since found out 51 that better agreement is reached when using a few more orbitals). The larger number of virtual orbitals required in the pyridine case when the diffuse basis is used set can probably be ascribed to the fact that these orbitals (due to the absence of the second nitrogen atom) 10

ACS Paragon Plus Environment

Page 11 of 26

The Journal of Physical Chemistry

are more diffuse than those of the diazines. It should be noted that the effect on the cross 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

sections of varying the number of virtual orbitals is weak and smaller than the differences due to employing the compact or diffuse basis sets.

Results Cross sections Figure 1 shows the integral elastic cross section for energies up to 15 eV calculated with the R-matrix method using several models. Both (rotationally summed) Born corrected and uncorrected results are shown. Above the first excitation threshold, our R-matrix SEP calculations suffer from the presence of pseudo-resonances which appear as very narrow spikes in the cross section as can be seen in the Figure. Comparison of the various R-matrix cross sections allows us to conclude that: (i) both compact and diffuse basis sets produce almost identical cross sections, the most noticeable difference being the position of the visible shape resonances (see below), (ii) including an extra partial wave in the R-matrix calculation does not affect the Born-corrected results. Schwinger multichannel results from Barbosa et al. 15 are also plotted in Figure 1. As is the case for the R-matrix results, both an uncorrected cross section and one corrected to include additional partial waves are presented. We can see in the Figure that for energies below 5 eV the uncorrected Schwinger cross section is smaller than the R-matrix one and fairly different in shape. A similar difference in the shape of the cross section obtained with both methods was found for the non-dipolar pyrazine and attributed to the differences in the description of the polarization. 51 We can see that the diffuse basis set R-matrix calculation shows some small indication of a decrease in the 2-3 eV range. Above 5 eV, the size and shape of both uncorrected cross sections is similar. There are a number of differences between the R-matrix and Schwinger calculations that could explain the difference in the uncorrected cross sections: use of different basis sets, of different configurations for the description of the 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Page 12 of 26

short-range polarization, different number of partial waves included in the expansion and 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the fact that R-matrix calculations do not take into account the polarization potential in the outer region. The Schwinger cross section obtained using a Born-closure procedure (’corrected’) is smaller than the R-matrix one in the whole energy range. This is partly due to the fact that the Schwinger integral cross section is obtained integrating the corresponding DCS between 1◦ and 180◦ , whereas the R-matrix cross sections correspond to the integration over the whole energy range (but they are not obtained from a numerical integration of the DCS). Integrating our DCS in the range 1◦ and 180◦ for 6 eV, results in a cross section of around 187 a20 , a smaller value than that of the corrected Schwinger cross section (211 a20 ). The difference between these two values, however, is the same as that between the uncorrected R-matrix and Schwinger cross sections. This means that the difference in the corrected cross sections is therefore mainly due to the different integration ranges. It also indicates that it is unlikely that a different number of partial waves included in the (uncorrected) R-matrix and Schwinger is responsable for the differences between these two cross sections. Figure 2 shows the elastic integral cross section for the whole energy range studied in this work. We see that the size of IAM-SCAR and IAM+Rotations cross sections around 10 eV is reasonably consistent with the R-matrix and Schwinger multichannel uncorrected and corrected cross sections respectively. The shape of the uncorrected cross section is different, but this is not surprising: the IAM-SCAR method is not expected to provide reliable results for energies much lower than ∼30 eV. Nonetheless, the slope of the IAMSCAR cross section above 30 eV and the low energy cross sections below ∼6 eV indicate that a smooth interpolation between both energy ranges is possible. This agreement is similar to that obtained for pyrimidine. 13 The cross section for electron scattering from pyrimidine is similar in size to the pyridine one (the dipole moment of pyrimidine is fairly similar too: 2.334D 52 ). Also plotted in Figure 2 is the total electronically inelastic cross section obtained with the IAM-SCAR method. This cross section is somewhat bigger than

12

ACS Paragon Plus Environment

Page 13 of 26

The Journal of Physical Chemistry

that for pyrimidine. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

13

Note that this cross section is zero below ∼10 eV, when in fact, the

first electronic excitation thresholds is around 4 eV. This is due to the fact that in the present IAM-SCAR model, electronic excited states lying below 10 eV are ignored. It should be noted that the IAM-SCAR+Rotations cross section is obtained, as explained above, assuming a 300 K population for the initial rotational states. The Born-corrected Rmatrix cross section is determined assuming the molecule is initially in its ground rotational state and summed over all possible final states with J < 9. The comparison of these cross sections is, in principle, not ’like-with-like’. However, recent work on positron scattering from pyrimidine 53 showed that for a rigid rotor with parameters similar to those of pyridine, temperature effects on the cross section are very small below 20 eV. Figure 3 shows DCS for several energies. We observe minor difference between the Rmatrix DCS calculated with different basis sets. For the lower energies we have plotted the R-matrix results together with those of Barbosa et al. 15 The agreement between R-matrix and Schwinger results is good, and fairly independent of the energy. Note that the agreement for small angles (≤ 10◦ ) looks very good, but this is due to the scale of the figure; in fact, although the shape of the DCS is very similar, the size is not. The Schwinger cross section is smaller than the R-matrix one. For 10 eV we have also plotted the IAM-SCAR results: once again, the differences in the shape of ab initio and IAM-SCAR DCS are not surprising as this energy is beyond the range of validity of the method. At higher energies, only the IAM-SCAR/IAM-SCAR+Rotations method is able to provide results. The shape and size of the DCS is very similar to that of pyrimidine below 50 eV. Above this energy, however, the DCS for pyridine seems to be much bigger.

Low energy results: shape resonances The integral elastic cross section in Figure 1 shows clearly the presence of three resonances in both the uncorrected R-matrix and Schwinger results (the resonant peaks are less obvious in 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Table 2: Positions (and widths), in eV, of the π resonances in pyridine calculated in this work at the SEP level using the basis sets and number of virtual orbitals (Nr. V) indicated. The parameters have been determined using RESON, 54 a Breit-Wigner profile fitting program. Also listed are the positions of the resonances calculated by Barbosa et al. 15 at the SEP level considering only singlet-coupled excitations and the experimental results of Nenner and Schulz 17 and Modelli and Burrow. 16 ∗

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nr. V Basis set cc-pVDZ 25 Basis set cc-pVDZ 35 Basis set cc-pVDZ 40 Basis set 6-311+G** 50 Barbosa et al. Nenner and Schulz Modelli and Burrow

2

B1 1.15 (0.057) 0.77 (0.022) 0.48 (0.01) 0.67 (0.027) 0.90 0.62 0.72

2

2 A2 B1 1.46 (0.058) 5.8 (0.60) 1.11 (0.025) 5.51 (0.48) 0.96 (0.016) 5.32 (0.48) 1.07 (0.030) 5.33 (0.47) 1.33 5.80 1.20 4.58 1.18 4.48

the Born corrected cross sections but visible nonetheless) below 6 eV. The pseudoresonances present in the R-matrix results prevent us from providing reliable information on any physical resonances that may appear above this energy region. In the diazines, both calculations and experiments indicate the presence of a fourth π ∗ resonance. 46 We are unable to determine whether a fourth π ∗ resonance is also present in electron-pyridine scattering. The resonance positions are listed in Table 2, together with results from ETS experiments. 16,17 Also listed is the symmetry of the resonances. For the calculations with the compact basis set, we have listed results with three different numbers of virtual orbitals. As already stated, 35 virtual orbitals give the best agreement with the results of Nenner and Schulz. Increasing the number of virtual orbitals included in the calculations with the compact basis set to 40 causes the two lowest resonances to become overcorrelated. The R-matrix resonances appear somewhat lower in energy than the Schwinger ones (addition of triplet-coupled excitations has been shown to lower the position of the 2 B1 resonances in the Schwinger calculations). The diffuse basis set seems to give positions in better agreement with both experiments. The third resonance appears too high in energy in all the calculations quoted: this resonances is likely to have a mixed shape and core-excited character (as does the third resonance in diazines 46,55 ) and therefore, the inclusion of configurations 14

ACS Paragon Plus Environment

Page 14 of 26

Page 15 of 26

The Journal of Physical Chemistry

describing electronically excited states of the target are needed to describe it accurately. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Conclusions We report elastic cross sections for electron scattering from pyridine in the energy range 1 eV to 1 keV using the R-matrix (low energy) and IAM-SCAR (higher energy) methods, together with electronically inelastic cross sections calculated with the IAM-SCAR method. These cross sections agree reasonably well. Earlier, Schwinger calculations, produce cross sections that are significantly smaller than our Born-corrected cross sections. We show that this discrepancy is mainly due to the choice of angular range ( 1◦ and 180◦ ) for which the Schwinger DCS is integrated to obtain the integral cross section. However, the uncorrected cross sections also show some quantitative and qualitative differences that we believe are mainly due to differences in the description of the polarization. Three shape resonances (of π ∗ character) are identified and characterized: their positions agree well with experimental results and prior calculations. The elastic pyridine cross sections are fairly similar in size to those of pyrimidine for energies below 50 eV. There are, however, however small differences in the shape of the DCS for angles > 30 − 40◦ . Pyrimidine has an additional nitrogen in the 6-atom ring and therefore an additional lone pair. This may increase slightly the exchange effects and lead to the changes in the DCS. Above 50 eV, the DCS seem to be larger for pyridine. Similarly the inelastic cross section obtained with the IAM-SCAR method is somewhat bigger than that for pyrimidine. On the whole, however, the substitution of a carbon atom for a nitrogen one in going from pyridine to pyrimidine seems to have a small effect on the scattering data presented here.

15

ACS Paragon Plus Environment

2

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Cross sections (a0 )

The Journal of Physical Chemistry

Page 16 of 26

compact compact (l=5) diffuse Ab intio compact (l=5) Ab intio diffuse Schwinger corrected Schwinger

1000

100 0

1

2

3

4

5

6

7 8 E (eV)

9

10

11

12

13

14

15

Figure 1: Low energy integral elastic cross section for electron scattering from pyridine calculated using the R-matrix method using the basis sets indicated in the figure. The cross sections labelled ’Ab initio’ and ’Schwinger’ do not incorporate a Born-based correction; for those labelled l=5 the ab initio calculation contains the extra partial wave. Also plotted are the multichannel Schwinger results of Barbosa et al..

16

ACS Paragon Plus Environment

Page 17 of 26

1000

2

Cross sections (a0 )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

100

10

1 0.1

diffuse Ab initio diffuse Schwinger Schwinger corrected IAM-SCAR IAM-SCAR+Rotations Inelastic (IAM-SCAR) 1

10 E (eV)

100

1000

Figure 2: Integral elastic cross section for electron scattering from pyridine for the whole energy range studied in this work for the methods and models indicated in the figure. Also plotted are the Schwinger results of Barbosa et al..The cross sections labelled ’Ab initio’ and ’Schwinger’ do not incorporate a Born-based correction.

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

3 eV

6 eV

100

100

10

10

1

1 R-matrix diffuse R-matrix compact IAM-SCAR+Rot. Schwinger

10 eV 100 100

-16

2

-1

DCS (10 cm sr )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 26

10

30 eV

10 1 1

50 eV

100 eV 200 eV 300 eV 1000 eV

100

100

10 10 1 1

0.1

0.1 0.01 0

30

60

90

120 150 180

0

30

E (eV)

60

90

120 150 180

E (eV)

Figure 3: Elastic differential cross section for electron scattering from pyridine for the energies indicated in the panels calculated using the R-matrix and IAM-SCAR+Rotations methods. The R-matrix results are Born-corrected. Also plotted are the multichannel Schwinger results of Barbosa et al.. The bottom right-hand panel shows IAM-SCAR+Rotations cross sections for several energies.

18

ACS Paragon Plus Environment

Page 19 of 26

The Journal of Physical Chemistry

Figures 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Acknowledgement This work was partially supported by the Spanish Ministerio de Econom´ıa y Competitividad (Project FIS2012-31230) and the European Commission COST Action MP1002 (NanoIBCT) and Action CM1301 (CELINA). The research carried out at The Open University was funded by EPSRC. We thank Prof. Marcio Bettega for providing us with their calculated data in order to produce the plots.

References ˇ ˇ (1) Carsky, P., Curik, R., Eds. Low-Energy Electron Scattering from Molecules, Biomolecules and Surfaces; CRC Press, 2011. (2) Baccarelli, I.; Bald, I.; Gianturco, F. A.; Illenberger, E.; Kopyra, J. Electron-induced damage of DNA and its components: Experiments and theoretical models. Phys. Rep. 2011, 508, 1–44. (3) Bouda¨ıffa, B.; Cloutier, P.; Hunting, D.; Huels, M. A.; Sanche, L. Resonanant formation of DNA strand breaks by low-energy (3 to 20 eV) electrons. Science 2000, 287, 1658– 1660. (4) Alizadeh, E.; Sanche, L. Precursors of Solvated Electrons in Radiobiological Physics and Chemistry. Chemical Reviews 2012, 112, 5578–5602. (5) Gallup, G. A.; Fabrikant, I. I. Vibrational Feshbach resonances in dissociative electron attachment to uracil. Phys. Rev. A 2011, 83, 012706. (6) Agostinelli, S. et al. Geant4–a simulation toolkit. Nucl.Instrum.Methods Phys. Res.,Sect.A 2003, 506, 250 – 303.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Page 20 of 26

(7) Pimblott, S. M.; LaVerne, J. A. Production of low-energy electrons by ionizing radia1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

tion. Radiation Physics and Chemistry 2007, 76, 1244 – 1247, Proceedings of the 11th Tihany Symposium on Radiation Chemistry. (8) Bouchiha, D.; Gorfinkiel, J. D.; Caron, L. G.; Sanche, L. Low-energy electron collisions with tetrahydrofuran. J. Phys. B: At. Mol. Opt. Phys. 2006, 39, 975–986. (9) Trevisan, C. S.; Orel, A. E.; Rescigno, T. N. Elastic scattering of low-energy electrons by tetrahydrofuran. J. Phys. B: At. Mol. Opt. Phys. 2006, 39, L255–L260. (10) Tonzani, S.; Greene, C. Radiation damage to DNA: Electron scattering from the backbone subunits. Journal of Chemical Physics 2006, 125 . (11) Fuss, M.; Munoz, A.; Oller, J. C.; Blanco, F.; Almeida, D.; Limao-Vieira, P.; Do, T. P. D.; Brunger, M. J.; Garc´ıa, G. Electron-scattering cross sections for collisions with tetrahydrofuran from 50 to 5000 eV. Phys. Rev. A (Atomic, Molecular, and Optical Physics) 2009, 80, 052709–6. (12) Gauf, A.; Hargreaves, L. R.; Jo, A.; Tanner, J.; Khakoo, M. A.; Walls, T.; Winstead, C.; McKoy, V. Low-energy electron scattering by tetrahydrofuran. Phys. Rev. A 2012, 85, 052717. (13) Sanz, A.; Fuss, M.; Blanco, F.; Maˇs´ın, Z.; Gorfinkiel, J. D.; Carelli, F.; Sebastianelli, F.; Gianturco, F.; Garc´ıa, G. Electron scattering cross section calculations for polar molecules over a broad energy range. Applied Radiation and Isotopes 2014, 83, Part B, 57 – 67. (14) Maˇs´ın, Z.; Gorfinkiel, J. D. Resonance formation in low energy electron scattering from uracil. The European Physical Journal D 2014, 68, 1–8. (15) Barbosa, A. S.; Pastega, D. F.; Bettega, M. H. F. Shape resonances in the elastic scattering of slow electrons by pyridine. Phys. Rev. A 2013, 88, 022705. 20

ACS Paragon Plus Environment

Page 21 of 26

The Journal of Physical Chemistry

(16) Modelli, A.; Burrow, P. Electron-transmission study of the temporary anion states of 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

substituted pyridines. Journal of Electron Spectroscopy and Related Phenomena 1983, 32, 263 – 276. (17) Nenner, I.; Schulz, G. J. Temporary negative ions and electron affinities of benzene and N-heterocyclic molecules: pyridine, pyridazine, pyrimidine, pyrazine, and s-triazine. J. Chem. Phys. 1975, 62, 1747. (18) Burke, P. G. R-Matrix Theory of Atomic Collisions: Application to Atomic, Molecular and Optical Processes; Springer, 2011. (19) Tennyson, J. Electron - molecule collision calculations using the R-matrix method. Phys. Rep. 2010, 491, 29–76. (20) Carr, J. M.; Galiatsatos, P. G.; Gorfinkiel, J. D.; Harvey, A. G.; Lysaght, M. A.; Madden, D.; Maˇs´ın, Z.; Plummer, M.; Tennyson, J.; Varambhia, H. N. UKRmol: a lowenergy electron- and positron-molecule scattering suite. The European Physical Journal D 2012, 66 . (21) Sanna, N.; Gianturco, F. A. Differential cross sections for electron/positron scattering for polyatomic molecules. Computer Phys. Comm. 1998, 114, 142–167. (22) Norcross, D. W.; Padial, N. T. The Multipole-Extracted Adiabatic-Nuclei Approximation for Electron-Molecule collisions. Phys. Rev. A 1982, 25, 226–338. (23) Gianturco, F. A.; Jain, A. The theory of electron scattering from polyatomic molecules. Phys. Rep. 1986, 143, 347–425. (24) Murai, H.; Ishijima, Y.; Mitsumura, T.; Sakamoto, Y.; Kato, H.; Hoshino, M.; Blanco, F.; Garc´ıa, G.; ao Vieira, P. L.; Brunger, M. J.; Buckman, S. J.; Tanaka, H. A comprehensive and comparative study of elastic electron scattering from OCS and CS2 in the energy region from 1.2 to 200 eV. J. Chem. Phys. 2013, 138, 054302. 21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Page 22 of 26

(25) Brunton, J. R.; Hargreaves, L. R.; Buckman, S. J.; Garc´ıa, G.; Blanco, F.; Zat1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

sarinny, O.; Bartschat, K.; Brunger, M. J. Anomalously large low-energy elastic cross sections for electron scattering from the CF3 radical. Chem. Phys. Lett. 2013, 568-569, 55. (26) Palihawadana, P.; Sullivan, J. P.; Buckman, S. J.; Maˇs´ın, Z.; Gorfinkiel, J. D.; Blanco, F.; Garc´ıa, G.; Brunger, M. J. A joint theoretical and experimental study for elastic electron scattering from 1,4-dioxane. J. Chem. Phys. 2013, 139, 014308. (27) Cowan, R. The Theory of Atomic Structure and Spectra; University of California Press: London, 1981. (28) Riley, M.; Truhlar, D. Approximations for the exchange potential in electron scattering. J. Chem. Phys. 1975, 63, 2182. (29) Zhang, X. Z.; Sun, J. F.; Liu, Y. F. A new approach to the correlation polarization potential-low-energy electron elastic scattering by He atoms. J. Phys. B 1992, 25, 1893. (30) Staszewska, G.; Schwenke, D.; Thirumalai, D.; Truhlar, D. Quasifree-scattering model for the imaginary part of the optical potential for electron scattering. Phys. Rev. A 1983, 28, 2740. (31) Blanco, F.; Garc´ıa, G. Improvements on the quasifree absorption model for electron scattering. Phys. Rev. A 2003, 67, 022701. (32) Blanco, F.; Garc´ıa, G. Screening corrections for calculation of electron scattering from polyatomic molecules. Phys. Lett. A 2003, 317, 458. (33) Zatsarinny, O.; Bartschat, K.; Garc´ıa, G.; Blanco, F.; Hargreaves, L. R.; Jones, D. B.; Murrie, R.; Brunton, J. R.; Brunger, M. J.; Hoshino, M.; Buckman, S. J. Electroncollision cross sections for iodine. Phys. Rev. A 2011, 83, 042702.

22

ACS Paragon Plus Environment

Page 23 of 26

The Journal of Physical Chemistry

(34) de Oliveira, E. M.; da Costa, R. F.; Sanchez, S. d.; Natalense, A. P. P.; Bettega, M. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

H. F.; Lima, M. A. P.; Varella, M. T. d. N. Low-energy electron scattering by cellulose and hemicellulose components. Phys. Chem. Chem. Phys. 2013, 15, 1682–1689. (35) Blanco, F.; Garc´ıa, G. Screening corrections for calculation of electron scattering differential cross sections from polyatomic molecules. Phys. Lett. A 2004, 330, 230. (36) Blanco, F.; Garc´ıa, G. Calculated cross sections for electron elastic and inelastic scattering from DNA and RNA bases. Phys. Lett. A 2007, 360, 707. (37) Jain, A. Theoretical study of the total (elastic+inelastic) cross sections for electron −H2 O (NH3 ) scattering at 10-3000 eV. J. Phys. B 1988, 21, 905. (38) A. Mu˜ noz, J. C. O.; Blanco, F.; Gorfinkiel, J. D.; ao Vieira, P. L.; Garc´ıa, G. Electron scattering cross sections and stopping powers in H2 O. Phys. Rev. A 2007, 76, 052707. (39) Sanz, A. G.; Fuss, M. C.; Blanco, F.; Sebastianelli, F.; Gianturco, F. A.; Garc´ıa, G. Electron scattering cross sections from HCN over a broad energy range (0.110000 eV): Influence of the permanent dipole moment on the scattering process. J. Chem. Phys. 2012, 137, 124103. (40) Fuss, M. C.; Sanz, A. G.; Blanco, F.; Oller, J. C.; Limao-Vieira, P. P.; Brunger, M. J.; Garc´ıa, G. Total electron-scattering cross sections from pyrimidine as measured using a magnetically confined experimental system. Phys. Rev. A 2013, 88, 042702. (41) Dickinson, A. S. Differential cross sections for electron scattering by strongly polar molecules. J. Phys. B 1977, 10, 967. (42) D. Johnson III Editor, R. NIST Standard Reference Database Number 101. 2011; http: //cccbdb.nist.gov/. (43) R. D. Nelson Jr., A. A. M., D. R. Lide Selected Values of electric dipole moments for molecules in the gas phase; 1967. 23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Page 24 of 26

(44) Schreiber, M.; Silva-Junior, M. R.; Sauer, S. P. A.; Thiel, W. Benchmarks for elec1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

tronically excited states: CASPT2, CC2, CCSD, and CC3. The Journal of Chemical Physics 2008, 128, 134110–25. (45) Maˇs´ın, Z.; Gorfinkiel, J. D.; Jones, D. B.; Bellm, S. M.; Brunger, M. J. Elastic and inelastic cross sections for low-energy electron collisions with pyrimidine. The Journal of Chemical Physics 2012, 136, 144310. (46) Maˇs´ın, Z.; Gorfinkiel, J. D. Shape and core excited resonances in electron collisions with diazines. The Journal of Chemical Physics 2012, 137, 204312. (47) Sanz, A. G.; Fuss, M. C.; Blanco, F.; Gorfinkiel, J. D.; Almeida, D.; da Silva, F. F.; Limao-Vieira, P.; Brunger, M. J.; Garc´ıa, G. An investigation into electron scattering from pyrazine at intermediate and high energies. The Journal of Chemical Physics 2013, 139, 184310. (48) Werner, H.-J.; Knowles, P. J.; Lindh, R.; Manby, F. R.; Sch¨ utz, M.; et al., MOLPRO, version 2009.1, a package of ab initio programs. 2009. (49) Faure, A.; Gorfinkiel, J. D.; Morgan, L. A.; Tennyson, J. GTOBAS: fitting continuum functions with Gaussian-type orbitals. Comp. Phys. Commun. 2002, 144, 224–241. (50) Tarana, M.; Tennyson, J. Polarization effects in electron collisions with Li2 : application of the molecular R-matrix method with pseudostates. J. Phys. B: At., Mol. Opt. Phys. 2008, 41, 205204. (51) Maˇs´ın, Z. Resonance formation in electron collisions with pyrimidine-like targets. Ph.D. thesis, The Open University, 2012. (52) Blackman, G.; Brown, R.; Burden, F. Microwave spectrum, dipole moment, and nuclear quadruple coupling constants of pyrimidine. Journal of Molecular Spectroscopy 1970, 35, 444. 24

ACS Paragon Plus Environment

Page 25 of 26

The Journal of Physical Chemistry

(53) Franz, J.; Gianturco, F. A. Low-energy positron scattering from gas-phase pyrimidine: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

A quantum treatment of the dynamics and a comparison with experiments. Phys. Rev. A 2013, 88, 042711. (54) Tennyson, J.; Noble, C. J. RESON–A program for the detection and fitting of BreitWigner resonances. Comput. Phys. Commun. 1984, 33, 421–424. (55) Winstead, C.; McKoy, V. Resonant Channel Coupling in Electron Scattering by Pyrazine. Phys. Rev. Lett. 2007, 98, 113201.

25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Graphical TOC Entry 1000

2

Cross sections (a0 )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

100

10

1 0.1

diffuse Ab initio diffuse Schwinger Schwinger corrected IAM-SCAR IAM-SCAR+Rotations Inelastic (IAM-SCAR)

1

10 E (eV)

100

1000

Integral elastic cross section for

electron scattering from pyridine.

26

ACS Paragon Plus Environment

Page 26 of 26

Electron scattering from pyridine.

We have calculated cross sections for elastic and inelastic electron scattering from pyridine in the energy range 1 eV to 1 keV. The R-matrix and IAM-...
326KB Sizes 2 Downloads 4 Views