R-matrix calculations for elastic electron and positron scattering from pyrazine: effect of the polarization description

We present R-matrix calculations of electron and positron low energy scattering from the highly polarizable pyrazine molecule. We compare integral and differential elastic cross sections with experimental results and assess the quality of the models used for describing collisions of either projectile. Static-exchange-plus-polarization models give a good description of electron scattering (including that of shape resonances), whereas both the integral and small-angle differential cross sections are underestimated for positron collisions for the same models. We discuss whether the absence of a permanent dipole moment improves the comparison with experiment for this molecule, as well as future calculations that may improve the description of polarization effects and thus positron scattering.


Introduction
Low energy electron and positron collisions with molecules take place in a number of natural and manmade environments with the former being, due to the much higher abundance of electrons, significantly more prevalent and therefore the focus of more computational investigations [1]. These works, together with their experimental counterparts, are aimed at providing both quantitative data and fundamental insight into the collisions to help model and control the environments in which they are involved (technological plasmas, radiation used for medical treatment and imaging, etc.). Below the positronium formation threshold, the same approach can be used to model both types of collisions: e.g. the R-matrix [2], Kohn variational [3], Schwinger Multichannel (SMC) [4] methods, single centre approach [5] and convergent close-coupling approach (for H 2 ) [6]. Their computational implementations can be applied to the study of the low energy interaction of both projectiles with small and medium size molecules. In general, however, these approaches seem to fare better when describing electron scattering, with positrons collision (even below the positronium formation threshold) not always well described.
The lack of exchange interaction between projectile and target electrons in positron scattering could be seen as making the computational modelling of these collisions easier than for electron scattering. This is, however, not the case [2,7,8]. On one hand, since polarization is an attractive interaction for both projectiles, this a e-mail: J.Gorfinkiel@open.ac.uk (corresponding author) interaction has the same sign as the (short-range) static interaction in the case of an electron, but opposite sign for a positron; this makes the observables for the latter more sensitive to the polarization description. On the other, the short-range correlation in positron scattering seems to be harder to model than the (anti)correlation of electrons.
As a result, the accurate description of correlation/polarization effects is much more critical in positron calculations. Whereas on the whole, this description no longer seems to pose serious limitations to the modelling of electron scattering, this is not the case for positrons, and this holds true for several of the methods mentioned above. Attempts have been made to improve and modify the approaches: for example, the use of atomic basis functions placed on additional 'dummy' centres in the SMC method [9] or use of a semi-empirical scaling factor to increase the electronpositron attraction integrals in R-matrix calculations [8] that improved their accuracy for a range of molecules [10,11]. Use of pseudostates has also been tried in Rmatrix calculations to improve polarization description.
In this paper, we investigate the description of correlation/polarization effects in calculations for positron and electron scattering from the non-polar pyrazine (see Fig. 1) using the R-matrix method. As with all diazines, pyrazine has a polarizability of around 60 a 3 0 . With its high polarizability and the dipolar interaction absent, pyrazine is an excellent system to study the effect of polarization description: the low energy behaviour of the cross section will be strongly affected by it and the presence of a Ramsauer-Townsend (RT) minimum and/or a virtual state is possible. Pyrazine: blue indicates nitrogen atoms, dark grey carbon and white hydrogen. The coordinate system used in our calculations is indicated on the right A couple of comprehensive theoretical studies have been performed into electron scattering from pyrazine. Mašín and Gorfinkiel [12] used the R-matrix method to determine elastic and inelastic integral cross sections, elastic differential cross sections and resonance positions. Winstead and McKoy [13,14] also determined elastic cross sections as well as shape resonance energies and widths using the SMC method. Experimental elastic absolute differential and integral cross sections were measured by Palihawadana et al. [15] using a crossed electron-molecular beam spectrometer technique. Resonances in pyrazine were investigated by Nenner and Schulz [16]. The only available publications on positron scattering from pyrazine are a computational study by Moreira and Bettega [17], who used the SMC method to determine the elastic integral and differential cross sections and a recent joint theoretical/experimental work [18] that presents measured elastic, electronic excitation and positronium formation cross sections for scattering energies of 1 to 79 eV as well as cross sections calculated using the R-matrix method (elastic, low energy) and the IAM-SCAR approach (elastic and inelastic, higher energies).
The R-matrix method is described in Sect. 2, emphasizing polarization description. Electron scattering results, based on the work of Mašín and Gorfinkiel [12] are presented in Sect. 4.1 and their agreement with SMC calculations and experiments is discussed. Sect. 4.2 presents our positron calculations and, again, a critical comparison with prior data. Section 4.3 discusses the use of pseudostates to improve the polarization description. Finally, further discussion and conclusions are presented in Sect. 5.

The R-matrix method
The R-matrix method has been extensively and successfully applied to the study of electron-molecule scattering below the ionization threshold. A number of publications describe it in detail [2,19]. The calculations in this work have been performed with the UKRmol suite [20] for positrons and the UKRmol+ suite [21] for electrons as the latter was not yet fully tested for positrons. Below, we provide a brief summary of the method and refer the reader to earlier publications for more details. We describe and apply the approach within the fixednuclei approximation.
The R-matrix method separates the scattering problem into an inner and an outer region by means of a sphere of radius a. In the inner region, all particles are considered explicitly and exchange and correlation must be taken into account. By contrast, in the outer region the scattering particle becomes distinguishable (in the case of the electron) so exchange can be neglected and the projectile-target interaction is described using a single centre multipole expansion potential. This makes the outer region computationally simpler and cheaper for most calculations. At the boundary of these two regions, the R-matrix is constructed using results from the inner region.
The inner region basis functions Ψ Γ k (x 1 , ..., x N +1 ) can be expanded using the Close-Coupling approximation: Here, x i denotes the four spin and space coordinates of the electron i. The operator A ensures antisymmetrization of the wavefunctions and is required when the projectile is an electron. The target wavefunctions Φ i and the L 2 functions χ Γ i are built as antisymmetrized products of bound molecular orbitals. The L 2 functions describe the short-range correlation-polarization. The single particle functions γ ij are known as continuum orbitals and describe the projectile. The coefficients a ijk and b ik are determined through the diagonalization of the Hamiltonian matrix in the inner region, made Hermitian by the inclusion of the Bloch operator [2]. n b is the number of target electronic states included in the calculation and the superscript Γ indicates the irreducible representation of the functions: Φ i (x 1 , ..., x N ) and γ ij must be chosen such that their product has the correct space and spin symmetry.
In this work, the continuum orbitals γ ij are constructed from Gaussian type orbitals (GTOs) centred on the centre of mass of the system. The individual continuum GTOs are first orthogonalized to the target orbitals and then to each other. A deletion threshold is required in the latter orthogonalization step. If this deletion threshold is set too high, too many continuum orbitals are deleted leading to a poor continuum description. Alternatively, if set too low, too many orbitals are kept, which can cause linear dependency issues.
The R-matrix, evaluated at the R-matrix boundary, is used to solve a simpler set of differential questions in the outer region. This is done using a propagation approach [2]. Matching the solutions obtained with their known asymptotic form enables the determination of the scattering energy-dependent K-matrix, that contains all the scattering information produced by the calculations. From the K-matrices, it is trivial to deter-mine T-and S-matrices and from them integral and differential cross sections. The latter are calculated using POLYDCS [22].
The type and number of L 2 functions, along with the number of target states included, defines the scattering model used. Well defined models regularly employed for electron and positrons scattering are the static (S)/static-exchange (SE) (for positrons and electron, respectively) and Static plus polarization (SP)/staticexchange plus polarization (SEP) models. Both of these include a single target wavefunction in the closecoupling expansion (Eq. 1). When more target state wavefunctions are used, the model is called closecoupling (CC). We briefly explain how polarization is described in these models.

Description of polarization
Polarization effects are modelled in two ways in Rmatrix calculations; (i) by inclusion of the appropriate L 2 functions (at S(E)P and CC levels); (ii) by inclusion of excited states in the CC expansion (CC level). The L 2 configurations used in S(E)P calculations for an N electron target can be schematically written as follows; Core' indicates a frozen core, containing N c electrons in doubly occupied orbitals, 'Val' a 'valence' space from which excitations are carried out and VO a set of virtual orbitals, i.e. bound molecular orbitals not included in the target ground state configuration. The VO are occupied by the projectile in the first set of configuration and by the projectile and a target electron in the second set. In addition, when the projectile is a positron, the positron is also allowed to occupy the Core and Val spaces.
In the case of the CC model, it is customary to expand the target wave functions in terms of a complete active space (CAS): Core Nc CAS N −Nc . The L 2 configurations correspond to allowing the projectile to occupy all the orbitals in the active space (again, the positron can also occupy a Core orbital): or one of a set of selected VO: Note that the Core spaces and N c in the CC and SEP/SP calculations need not be the same for a given target. Also, for few electron systems, usually N c =0 and no Core is used. When the polarizability of the molecule is not large, the second set of configurations is not required (this gives a more balanced model [21]).

Pseudostates and pseudocontinuum orbitals
An approach that has been employed to improve (in some cases, converge [23]) the polarization description for small molecules is the inclusion of pseudostates in the calculation. The R-matrix with pseudostates method (RMPS) was first developed to model electron-impact near-threshold ionization [24] and then employed to improve the modelling of positron scattering observables (cross sections, Z-effective) [25,26] as well as electron scattering below the ionization threshold (e.g. [27]).
Pseudostates are eigenfunctions of the molecular Hamiltonian that do not describe true states of the system. In the molecular RMPS approach [21,28] they are built from single excitations from the ground state configuration into a set of pseudocontinuum orbitals (PCOs). These PCOs are, in turn, generated from even tempered GTOs centred on the centre of mass of the system. Their exponents follow a geometrical progression: with 1 < β < 2. The exponents are the same for different angular momenta l but the number of GTOs used usually decreases as l increases (as is the case for the continuum basis set) up to a specified l PCO max . The configurations used to generate the target states and pseudostates can be written as: and the L 2 functions are thus: where, once more, the positron can also occupy Core and CAS orbitals (configurations of the type Core Nc CAS N −Nc−1 PCO 2 are not usually employed as they make the calculation too big). As mentioned above, earlier studies seem to show reasonable convergence of the polarizability when pseudostates are included. However, all molecules studied were relatively compact: diatomics, triatomics, or linear molecules like C 2 H 2 ; none had polarizabilities bigger than 25 a 3 0 , except Li 2 that has a polarizability over 200 a 3 0 .

Evaluation of the polarizability
It isn't possible, in R-matrix studies, to quantify the amount of polarization being described in S(E)P calculations. For CC calculations, one can determine the (implicit) value of the molecular polarizability being modelled in the inner region by the (pseudo)states by calculating it as a sum over said states; for the ground state the expression is: where q represents Cartesian components X, Y and Z, n and 0 label an excited and the target ground state, respectively and μ q is the component of the dipole operator along q. The summation runs over all the states included in the close-coupling expansion. Note that as the energy difference between the excited and ground states increases, the contribution of that state to the polarizability will decrease, leading to convergence. One can use this calculation to estimate how much polarization is being described in a CC calculation.

Details of the calculation
Pyrazine belongs to the D 2h point group and thus does not possess a permanent dipole moment. However, the molecule does have a large polarizability: an accurately computed (using the DFT method, the B3LYP functional and the 6-311+G(3df,2p) basis set) value for the spherical polarizability [29], along with its diagonal components are listed at the bottom of Table 3 (no experimental data is available). The geometry used in these calculations was taken from an experiment by Innes et al. [30]. The ground state configuration of pyrazine is: . Calculations were performed using a number of models. For the S(E)P calculations, HF orbitals were used to describe the ground state of pyrazine. The CC calculations used CASSCF orbitals determined with MOL-PRO: details of the active space and state-averaging can be found in [12]. Five basis sets were tested: two of them, cc-pVDZ and 6-311+G** have been used extensively in previous R-matrix studies of this and similar molecules, so that deletion thresholds for orthogonalization and the required R-matrix radius were known. For the remaining three basis sets, DZP, DZP+diffuse and 6-311++G(3df,3pd) these parameters had to be determined. The most appropriate number of VOs for each calculation was also determined and will be discussed below.
The appropriate R-matrix radius was determined from the orbital amplitudes computed as a function of the distance from the centre of mass using a functionality available in the UKRmol+ suite. The R-matrix radius is chosen to be the smallest possible that ensures the amplitudes of the target orbitals included in the calculation are less than approximately 10 −4 at the R-matrix boundary. As only GTOs were used for the continuum description, the largest R-matrix radius for which calculations could be sensibly run was 18 a 0 . This meant that for the 6-311++G(3df,3pd) basis set some orbitals had to be removed from the calculations as their amplitudes were larger than 10 −4 even for a=18 a 0 . Table 1 summarizes the R-matrix radius and deletion thresholds used for each basis set for which results are presented. The highest continuum partial waves for all calculations was l = 4.

Electron scattering
Mašín and Gorfinkiel [12] investigated a number of models for electron pyrazine scattering (different basis sets, number of virtual orbitals, active spaces for the CC calculation, etc.) focusing primarily on the description of the resonances. Their best SEP and CC calculations provided a fairly good description of the shape and mixed shape-core excited resonances present: the resonance positions were in reasonable agreement with experiment [16]. Table 2 summarizes the results for the shape resonances. Figure 2 shows the total elastic cross section calculated at SEP and CC levels using the two basis sets chosen by Mašín and Gorfinkiel for their detailed investigations: cc-pVDZ and 6-311+G**, compared to the SMC [13] and experimental results [15]. It is clear that the CC calculations provide very similar cross sections with small shifts in the π * resonance positions. The SEP results with the diffuse basis agree reasonably well with Winstead and McKoy's, whereas those using the compact basis set fail to describe the RT minimum (this is more clearly visible in the 2 A g contributions to the integral elastic cross section, presented in Figs. 3 and 5 Table 2 Energy (in eV) of the 3 lowest energy π * resonances in pyrazine. R-matrix data calculated at SEP level for the basis sets indicated [12]. Winstead [13] are presented as a purple dotted line. R-matrix calculations computed in this work using the models of Mašín and Gorfinkiel [12] and are shown as dashed and dot-dashed lines with the basis set and method used indicated in the figure of [12]). The position of this minimum in the total elastic cross section determined with the best SEP model is a little higher than in SMC calculations [31].
The comparison with experiment in Fig. 2 would lead to the conclusion that it's the CC calculations (and the SEP one with the compact basis set) that provide better agreement with experiment at 3 eV, potentially casting doubt on the presence of the RT minimum, and the conclusion, derived from the resonance positions, that the SEP model with the diffuse basis set provides the best agreement with experiment. However, the quantity actually measured by the experiment is the differential cross section (DCS); below 15 eV the measurements were performed for angles between 30 • and 129 • , so the integral cross sections in Fig. 2 are the result of extrapolation of these DCS in order to integrate over all angles. A better comparison of experiments and calculation is therefore given by the DCS. Figure 3 shows excitation functions (DCS for a specific scattering angle as a function of energy). Note that the jagged look of the SEP results is due to the presence of non-physical resonances above ∼6 eV that also lead to the overestimation of the cross section in that energy range. This, together with the fact that inclusion of excited states in the CC calculations usually leads to a somewhat lower elastic cross section than that produced by SEP calculations (and the appearance of physical core-excited resonances), accounts for the noticeable differences between calculated excitation functions; the shift in physical resonances due to different polarization descriptions is clearly visible for 90 • and 120 • . Nonetheless, we see that the excitation functions calculated at CC and SEP level with the diffuse basis set show similar agreement with those measured  [15] are shown as brown data points with angular distribution measurements in blue. The solid and dashed lines are R-matrix calculations computed in this work using the models of Mašín and Gorfinkiel [12] by Palihawadana et al., with the SEP results being slightly closer to experiment. Therefore, the difference in the integral cross section at 3 eV must come from the contribution of angles below 60 • and/or above 120 • .
To investigate this in more detail, Fig. 4 shows the DCS for 3 eV, both for positron and electron scattering. We see clearly here than in the case of electron scattering (the positron one will be discussed in the next section), the experimental DCS at 3 eV agrees quite well with both the SMC results and the R-matrix ones determined at SEP level with the diffuse basis set (the other R-matrix models, not shown, produce DCS that are significantly higher in the angular range 30 • -60 • ). This confirms that, at 3 eV, the higher value of the experimental integral cross section does not come from the measured values of the DCS; the extrapolation of the experimental values below 30 • and/or above 130 • must overestimate the DCS. This comparison also confirms that the RT minimum is probably physical. Elastic differential cross section for electron and positron scattering from pyrazine at 3 eV. Experimental data: blue circles, positron data by Edwards et al. [18] (the measurement at 10 • is outside the plotted range); green squares, electron data by Palihawadana et al. [32]. SMC SEP/SP results: purple dot-dashed curve, folded positron data by Moreira and Bettega [17]; cyan dot-dash-dashed curve, electron data by Winstead and McKoy [13]. R-matrix SEP/SP results: green dash-dot-dotted curve, electron scattering; black solid (unfolded) and red dashed (folded) curves, positron data. The basis sets used are indicated in the figure The best model identified for the description of elastic scattering (SEP using a diffuse basis set) describes electron scattering well: the resonance positions, the presence of an RT minimum and the size and shape of the DCS elastic cross sections are well described. It is therefore this model that we expect to provide the best description of positron scattering from pyrazine.

Positron scattering
Initial positron scattering calculations were based on the SEP models selected by Mašín and Gorfinkiel to study electron scattering from diazines. The SP elastic integral cross sections produced using both the compact and diffuse basis set, as well as the DZP basis, are shown in Fig. 5. One can see that the most diffuse basis set produces the largest cross section for the whole energy range presented and that increasing the number of VOs from that used for electron scattering (40) to 80 has a noticeable effect below 3-4 eV, whereas when increasing them to 120 the effect is significantly smaller.
We also see that the shape of the positron integral cross section is very different to that for electron scattering. The former increases significantly as the energy decreases, so that it is larger than the electron scattering one below ∼3 eV. This increase originates in the 2 A g symmetry contribution to the cross section and can be attributed to the presence of a virtual state. This behaviour is also observed in the experimental cross section [18] and that calculated by Moreira and Bettega  [17] and the cyan dots are experimental results from Edwards et al [18] using the SMC method [17], both shown in Fig. 6. Moreira and Bettega determined the scattering length (SL) associated with their results and concluded that it may actually be indicating a bound state, not a virtual state. The quality of the description of this state was dependent on how well the polarizability was described. Our estimates of the scattering length, determined from the A g eigenphase sum using the same equations provided by Morrison [33] employed by Moreira and Bettega, also point at a bound state, though more weakly bound than that of Moreira and Bettega: we obtain a scattering length of 190-200 a 0 (depending on the energy range fitted) for the calculation using the 6-311+G** basis set and 80 VOs and a value of 68-70 a 0 for the calculation using the 6-311+G** basis set and 120 VOs.
These values correspond to estimated binding energies (−1/(2*SL 2 )) of 0.34-0.38 meV and 2.8 meV, respectively. Given the approximate nature of these very small binding energies, they should not be taken as definite confirmation of the presence of a bound state. Figure 5 shows that the choice of target basis set has a significant effect, as was the case for electron scattering: the more compact DPZ and cc-pVDZ calculations produce cross sections that are much smaller than those determined with more diffuse basis sets. This effect is not too different in size to that of increasing VOs from 40 to 80: both influence the quality of the polarization description. This result led us to test the use of an even bigger basis set, 6-311++G(3df,3pd); this result is included in Fig. 6. The cross section shows some nonphysical structure at 3 and 4 eV, probably due to linear dependence or small leaking outside the sphere: some of the orbitals included had larger amplitudes at the boundary (7.7×10 −5 ) than for all other calculations. This prevented us from attempting to include more VOs in the calculation. Use of this basis set (and 120 VOs) does not bring the computed value of the cross section much closer to experiment, indicating that something else is required to ensure polarization and correlation effects are sufficiently well described.
We see in Fig. 6 that the R-matrix cross section calculated with the diffuse basis sets agree reasonably well with the SMC result, but that both methods significantly underestimate the experimental cross section. Interestingly, this disagreement is larger at higher energies, rather than at 1 eV where the dependence on the polarization description would be expected to be strongest. The difference for 1 eV is around 20% and close to 50% at 3 eV: for this latter energy the electron scattering cross section was underestimated by the SEP models by less than 30%, and this difference was attributed to the extrapolation used to determine the experimental integral cross sections. For 5 eV, the experimental positron cross section is about twice the calculated one, whereas for electrons, the differences at 6 eV are no bigger than 10%.
Unlike electron scattering, the positron experiment measures DCS and the integral cross section separately. The integral cross section is measured down to an energy-dependent angle: 13.7 • at 1 eV decreasing to 4.3 • at 10 eV. An energy-dependent 'correction' that corresponds to a percentage of the integrated cross section is employed to account for the missing angles (see [18] for more details). As a result, in the case of positrons it is not possible to make a direct connection between the angular range for which the DCS is measured and the experimental value of the integral cross section. Nonetheless, it is still instructive to compare the measured and calculated DCS. A look at the DCS for 3 eV in Fig. 4 shows that, as is the case for electrons, the SMC and R-matrix SP calculations with the diffuse basis agree reasonably for the whole angular range. The experimental DCS is significantly bigger than the calculated results between the first angle measured (10 • ) and ∼30 • (by almost a factor of 3 at 10 • ), whereas the calculated DCS are bigger than the exper-iment between ∼80 • and 90 • . Note that positron DCS measurements are 'folded' around 90 • : the folded DCS for a specific θ correspond to DCS(θ) + DCS(180 • -θ).
Nonetheless, Fig. 4 shows for 3 eV that the calculated DCS below 30 • is around an order of magnitude bigger than above 150 • , so that the differences between folded and unfolded DCS is less than 20% for small angles. Around 90 • , the DCS is similar in size below and above this angle, so that the folded DCS is around twice the unfolded one.
If one integrates the 3 eV folded experimental and R-matrix DCS in the angular range for which the DCS is measured, the experimental cross section is around 12% bigger: a significantly smaller difference than the factor of ∼2 observed in the integral cross section. Interestingly, if one performs the same partial integration for 1 eV, the experimental result is close to 12% smaller than the R-matrix one. In this reduced angular range, there is a compensation for positrons between the underestimation of the experimental DCS at small angles (where the sin(θ) is smaller) and its overestimation around 90 • . Once again, agreement between experimental and calculated DCS is better than between the integrated cross section: the large differences seen in the integral elastic cross section seem to arise either from the applied correction significantly overestimating the contribution of the missing angular range in the integral cross section measurement or from calculations underestimating the forward-peaked nature of the collisions. Perhaps both. The smaller disagreement at 1 eV may be due to the compensation mentioned earlier. Figure 7 provides a more detailed comparison of the different models, as it shows the contributions of the different D 2h irreducible representations to the integral cross section. The largest contribution below 1 eV is that of the 2 A g symmetry and it is clearly bigger for the more diffuse basis sets. Above 1 eV this contribution is smaller than 20Å 2 and the largest contribution is that of the 2 B 1u followed by the 2 B 3u symmetry for the diffuse bases but the 2 B 2u one in the case of the compact basis set. All other contributions are significantly smaller.
The agreement of the symmetry contributions for the diffuse basis sets with those calculated with the SMC approach is, overall, fairly good, although whether the 6-311+G** or the 6-311++G(3df,3pd) one produces better agreement is symmetry dependent. The more noticeable differences can be seen at very low energies, where the energy dependence of the cross section contributions is quite different for symmetries 2 B 1g and 2 B 2u (the spike at around 3 eV in this latter symmetry is a non-physical pseudoresonance as is the structure around 4 eV in 2 A g ). The symmetryby-symmetry agreement shows that the fact that both methods produce integral elastic cross sections of very similar size is not fortuitous. If the discrepancy with experiment is due to limitations in the calculation, these are clearly very similar for the models used in both approaches. Symmetry contributions to the integral elastic cross section for positron scattering from pyrazine using SP models. Red dot-dashed curve, cc-pVDZ basis set with 80 VOs; blue dotted curve, 6-311+G** basis set with 120 VOs and green solid curve, 6-311++G(3df,3pd) basis set using 120 VOs. Purple dashed curve, SP3 SMC model by Moreira and Bettega [17]

Testing PCO use
Inclusion of pseudostates in the calculation can improve the modelling of polarization effects so we investigated the effect of using of PCOs on the target polarization description: different l PCO max and number of pseudostates were used for the calculation of the polarizability using Eq. 3. The parameters to generate the PCO exponents (Eq. 2) and deletion thresholds (used in the orthogonalization of these orbitals) were chosen, following previous studies [25,34], to be α 0 = 0.17, β = 1.4 and 2 × 10 −4 , respectively. HF target orbitals were used for these tests. Table 3 summarizes the results. Even the largest calculations including PCOs and pseudostates (models 1 and 8) produce a spherical polarizability of 36-37 a 3 0 , less than 65% of the calculated accurate value of around 59 a 3 0 . One can also see that the effect of using a different target basis set is not large: about 1% for l PCO max =5 and 400 pseudostates, and around 10% for l PCO max =4 and 200 pseudostates The results also seem to show that increasing l PCO max has a small effect if the number of pseudostates included in the calculation remains the same (compare, for example, models 2 and 3). However, increasing the number of states for a specific l PCO max , does have a significant effect on the value of the polarizability (see, in this case, models 1 and 2). Model 9 shows the effect of using a different number of PCOs for different partial waves: here 8 PCOs are used for l ≤ 3 and 4 PCOs for l > 3. The spherical polarizability is almost identical to the one obtained when 8 PCOs are used for all l (model 8). This may point at a way of reducing the size of the integral calculation without a noticeable reduction in size of the polarization even for bigger models.
We note that from the practical point of view, increasing l PCO max implies a bigger integral calculation, while increasing the number of pseudostates increases the size of Hamiltonian to diagonalize and can also lead to the inclusion of more channels in the outer region. Nonetheless, it is clear from Table 3 that using more pseudostates leads to a better description of the polarization but that the best values obtained are still around 30% smaller than accurate results Given the unsatisfactory description of the polarizability of pyrazine when PCOs are used, we performed an additional type of test. This involves generating the same type of configurations as those used in the above calculations, but using VOs instead of PCOs. The configurations generated are no different to those used in standard SEP/SP calculations; the difference is that in this case these configurations are used to describe pseudostates in addition to being included as L 2 ones. We call this model Close-coupling Hartree-Fock (it is related to the Polarization Consistent Coupled Hartree-Fock, PCCHF, used in R-matrix photoionization calculations [35]), to indicate that the (pseudo)states included in the CC expansion are all described in terms of a single configuration/Slater determinant. It is clear from the results at the bottom of Table 3, generated using 80 VOs, that this produces significantly better values for the polarizabilities than the use of PCOs: use of 400 states produces a polarizability almost identical to the accurate value. Even use of 200 states gives a polarizability of around 44 a 3 0 , better than that obtained for the same target basis set and 400 states when PCOs are used (model 1). This model leads to a smaller calculation, both because fewer integrals need to be computed as no pseudocontinuum basis is included, and because fewer target states provide a better description of the polarizability. Electron scattering calculations with this model and 40 VOs showed a small low- For the models identified with a number, the total number of PCOs is obtained as 8 l PCO max (2l PCO max + 1) (except for model 9 that used 8 PCOs for l ≤ 3 and 4 PCOs for l > 3). The models labelled CC-HF did not use PCOs (see text for details). The number of CSFs is averaged over all irreducible representations. The accurate value is from the CCCBDB and was calculated using the DFT method B3LYP and 6-311+G(3df,2p) basis set [29]. The spherical polarizability and its components are given in a 2 0 . The X-, Y-and Z-axes correspond to those indicated in Fig. 1 ering of the π* resonances (around 0.25 eV) that turned the 2 B 3u into a bound state; it also leads to changes at very low energies in the A g contribution that might indicate the presence of a virtual state. We think it likely that both these states (the bound and virtual) are not physical and that this model is overcorrelated for electron scattering. Unfortunately, attempts to perform positron calculations using these CC-HF models lead to significant linear dependence. Use of a Bspline basis set for the continuum [21] may alleviate this problem.

Discussion and conclusions
We have performed R-matrix calculations for electron and positrons scattering from pyrazine and compared our integral and differential elastic cross sections with experiment and SMC calculations in order to assess how well polarization and correlation effects are described.
In the case of electron scattering, we have confirmed that discrepancies at low energies (3 eV) in the integral cross section are mainly due to the extrapolations required to generate the experimental integral value. Both experiment and theory produce a significantly more forward peaked DCS for positron than electrons. However, the calculated positron DCS does not increase as much as the angle decreases. It is unlikely the calculations are underestimating the backward scattering DCS by a factor close to 10, so it seems to be the modelling of the small angle (forward) scattering that, although qualitatively correct, it is quantitatively inconsistent with experiment for positrons. It is perhaps this underestimation of the forward-angle contribution to the calculated integral cross section (particularity for the very small angles, for which there are no measurements) that is a significant contributor to the disagreement of the positron integral cross sections.
Both for positron and electron scattering the SMC and R-matrix SEP/SP calculations produce integral elastic cross sections in very good agreement. We conclude that the models used in both calculations are of a similar quality: this is sufficient to capture the main effects influencing the collision for electron scattering, but this is not the case for positron as a projectile. Whether the issue is related to the description of polarization effects or correlation/anticorrelation is not possible to say as these can't be separated in our calculations.
When choosing pyrazine for these studies, an assumption was made that, due to the lack of dipolar interaction, the comparison with experiment would be more straightforward and therefore better than for a similar, but polar, target. The reason is that, for polar molecules, use of a Born-based top-up is necessary to fully account for the projectile-permanent dipole interaction [2] and model the behaviours of the DCS, that is very forward peaked, and the integral cross section, that grows rapidly as the energy tends to 0 eV. This top-up is understood to be less accurate than the ab initio calculations. In this respect, a comparison with results for pyrimidine, an isomer of pyrazine with a significant dipole moment, can shed some light.
The R-matrix method has been applied to the study of electron and positron scattering from pyrimidine. For electron scattering, elastic DCS agree very well with experiment for the measured angular range, except for very low energies (i.e. 1 eV) [36]; the Born-corrected integral cross section overestimated experiment, but the agreement was very good if the integral cross section was determined integrating the uncorrected DCS with their behaviour between 0 • and 2 • modified (see details in [36]). The corrected and uncorrected DCS are virtually identical for angles above a certain value (30 • for 3 eV decreasing to 10 • for 10 eV). So the better agreement when uncorrected DCS were used (with some modifications) seemed to indicate that the Born correction overestimates the forward scattering cross sections.
For pyrazine, it is the experiment that overestimates the integral cross section at low energies. For this target, the inaccuracy of the (unnecessary) Born correction cannot be the source of the discrepancy. The agreement between R-matrix results and experiment seems, overall, of a similar quality for both targets.
Positron scattering from pyrimidine has also been studied with the R-matrix method [37]: comparison with experiment [38] was analysed by Barbosa et al. [39] who also presented SMC results. The comparison at angles around 90 • is similar to that for pyrazine: the calculations overestimate the experimental DCS. At smaller angles, however, agreement is better: the experimental cross section is larger than the calculated ones, but the forward-peaked nature of the Born-corrected DCS means they are closer to the experimental values. As a result, the calculated integral elastic cross section is larger than the experimental one, the R-matrix one (calculated as SP level with a diffuse basis set) more so than the SMC and with the disagreement increasing as the scattering energy decreases. In this case, therefore, it seems that R-matrix method (and the SMC) describe the collisions with the dipolar target better than for pyrazine.
It is interesting to note that the models used for the electron-pyrimidine calculations mentioned above produce both accurate electronic excitation cross sections (to bands of states) and core-excited resonances: both agreed very well with high-quality experiments [40]. The position of the latter however, was shifted to higher energies in the calculations: this was interpreted to be the result of an insufficient description of the polarization effects.
Calculations of the polarizability of pyrazine using the sum-over-states expression and including a number of pseudostates and PCOs show that, at least for the PCOs and number of pseudostates tested, the approach does not significantly improve the value of the polarizability. A better approach for this seems to be to include in the CC expansion more target states built as single excitations from the ground state configuration to virtual orbitals (CC-HF model). This may prove a route to improving the polarization description and, perhaps agreement with experiments. Future work will focus on running these type of calculations, as well as standard CC ones, with the UKRmol+ suite, where use of B-splines for the continuum description should eliminate linear dependence problems; the suite is also computationally much more efficient making it easier to run bigger calculations with many more target states.