## Simulating Vibrationally-Resolved X-ray Photoelectron Spectra of Flexible Molecules: Linear Alkanes C H n 2 n +2 ( n =1-8)

Xiao Cheng, 1, 2, 3 Minrui Wei, 1 Guangjun Tian, 4 and Weijie Hua 1, ∗

1 MIIT Key Laboratory of Semiconductor Microstructure and Quantum Sensing, Department of Applied Physics, School of Physics, Nanjing University of Science and Technology, 210094 Nanjing, China 2 Hefei National Laboratory for Physical Science at the Microscale, University of Science and Technology of China, 230026 Hefei, China 3 Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, Royal Institute of Technology, SE-106 91, Stockholm, Sweden 4 Key Laboratory for Microstructural Material Physics of Hebei Province, School of Science, Yanshan University, Qinhuangdao 066004, China

We integrated full core-hole density functional theory with Franck-Condon calculations, considering Duschinsky rotation, to simulate vibrationally-resolved C1s X-ray photoelectron spectra (XPS) of eight linear alkanes, from methane to octane (C n H 2 n +2 , n =1-8). Results align excellently with experimental absolute binding energies and profiles. The spectrum of ethane serves as a 'spectral seed', with each longer alkane's atom-specific spectrum displaying similar characteristics, albeit with shifts and slight intensity adjustments. Detailed assignments clarify the distinct spectra in short alkanes ( n =1-4) and their stabilization in long alkanes ( n =5-8). Carbons are classified as central or distal (C 1 and C ), with central carbons contributing nearly identically to the lowest-energy feature 2 A, while distal carbons contribute to the second lowest-energy feature B (both are 0-0 transitions). Increasing molecule size adds more central carbons, enhancing feature A and weakening feature B. Our analysis identifies two dominant Franck-Condon-active vibrations: C -H stretching ( ∗ ∼ 34003500 cm -1 ) and bending ( ∼ 1400-1500 cm -1 ) modes. Structural analysis shows that core ionization minimally affects alkane geometry, less than in ring compounds from previous studies. This work extends our protocol from rigid ring compounds to flexible molecules, contributing to building a high-resolution theoretical XPS library and enhancing the understanding of vibronic coupling.

## I. INTRODUCTION

X-ray photoelectron spectroscopy (XPS) is a powerful analytical technique that provides detailed information about the chemical composition and electronic structure of molecules and materials.[1] High-resolution XPS spectra allow for investigation of fine structures, which predominantly originate from the vibronic coupling effect. The sub-eV fine structures carry crucial insights into the interplay between nuclear and electronic motions, unveiling the fundamental physics underpinning core hole ionization dynamics and alterations in potential energy surfaces (PESs) between ground and coreionized states. Historically, vibronic fine structures in XPS spectra were observed as early as 1974 by Siegbahn and coworkers[2] while studying the spectra of CH 4 . The continuous advancements in synchrotron radiation technologies have propelled XPS to higher resolutions, unlocking more comprehensive chemical information about the chemical systems.[3-11]

XPS databases have been helpful in interpreting the structure-spectroscopy relation. However, existing experimental XPS databases[12-17] primarily focus on binding energy (BE) values, often omitting intensity data, limiting the depth of analysis. Even for BEs, different experiments sometimes report values that differ over 1 eV. An ultimate goal is to create a comprehensive, high-resolution theoretical library for XPS spectra with both accurate BEs and profiles. Preliminary studies were conducted to investigate a family of molecules (azines[18], indoles[19], polycyclic aromatic hydrocarbons (PAHs)[20], and some other fiveand sixmembered ring molecules[21, 22], biomolecules[23]) with similar structures, by combing density functional theory (DFT) with the full core-hole (FCH) approximation and Franck-Condon (FC) simulations which incorporate the Duschinsky rotation (DR) effect[24] (denoted as the FCH-DR method for simplicity).[19, 23, 25] Grounded on the harmonic oscillation (HO) approximation, the calculations yielded good agreement with experiments.[18-21] The comparative analysis on theoretical and experimental high-resolution spectra offers a more rigorous evaluation, extending beyond solely assessing binding energies, for the accuracy of the chosen theoretical methodology, validating prediction of spectra of related molecules with no experimental data or no high-resolution experimental data. The studies provide data on equal footing, serving as a reliable stand point for analyzing vibronic properties within a specific family and edge. Extensive investigations were performed including influences on functionals, basis sets, temperature effects, and tautomeric effects, as well as analyses on core-ionization induced structural changes, active vibronal modes, etc. Additionally, the

∗ wjhua@njust.edu.cn

experiment-theory agreement also validated the HO approximation used which is less sensitive for polyatomic molecules than diatomic systems.[25]

However, all previous calculations on polyatomic systems have predominantly focused on ring compounds which exhibit relatively rigid structures including PAHs. It would be natural to investigate the performance on flexible molecules. In this study, we extend our research to a series of straight-chain alkanes (C n H 2 n +2 , n =1-8; see Fig. 1): methane (CH ), 4 ethane (C H ), 2 6 propane (C H ), 3 8 n -butane (C H 4 10 ), n -pentane (C H 5 12 ), n -hexane (C H 6 14 ), n -heptane (C H 7 16 ), and n -octane (C H 8 18 ). They represent the simplest organic molecule family, composed solely of sp 3 -hybridized carbon atoms and hydrogen atoms. Gradually increasing the -CH 2 groups provide an ideal set for investigating size effects. There are access to high-resolution experimental spectra[26] (except heptane), facilitating the performance evaluation of the theoretical method. Technically, as linear molecules generally possess greater flexibility than their ring counterparts, the geometry optimization of core-ionized states becomes a more intricate task than rigid molecules.

Alkanes are naturally occurring hydrocarbons found in crude oil, serving as high-value petrochemical feedstock materials for the petroleum industry[27]. They also play a prominent role as a major component of non-methane volatile organic compounds in vehicle emissions and urban atmospheric environments[28]. Methane, propane, and butane are clean energy sources for domestic and industrial applications, with methane dominating liquefied natural gas and propane/butane used in liquefied petroleum gas. Ethane, a minor natural gas component, holds significant value as a petrochemical feedstock[29, 30], primarily for ethylene production. Pentane is used as a polystyrene foam blowing agent and nonpolar solvent, while hexane is a common affordable, safe, and rapidly evaporable non-polar organic solvent. Heptane serves in various industrial applications, including paints, coatings, and pharmaceuticals, and is a significant gasoline component. Octane, another vital gasoline component, is often blended with the highly branched isomer iso-octane to establish a standard reference for measuring the fuel's resistance to engine knock. Methane exhibits T d symmetry, while ethane possesses D 3d symmetry. Longer straight-chain alkanes have either C 2v (propane, pentane, and heptane) or C 2h (butane, hexane, and octane) symmetry, as labeled in Fig. 1 together with all non-equivalent carbon atoms.

It is necessary to acknowledge early theoretical and experimental studies on vibrationally-resolved XPS of alkanes, particularly the work by Svensson and coworkers.[26, 31] They employed equivalent-core hole (ECH) DFT, replacing the core-ionized carbon (C ∗ ) with a nitrogen atom,[32-34] and Franck-Condon simulations to compute spectra for methane, ethane, and other small hydrocarbons, achieving good agreement with experiments.[31] Subsequently, they transitioned to another electronic structure method, using a nitrogen basis set for C ∗ with scaled exponents and an effective core potential (ECP) with adjusted parameters to account for one electron in the C1s shell.[26] This approach was applied to simulate larger linear alkanes ( n = 1 -6 8 , ) and extrapolated to alkanes with n = 1000 , providing valuable insights into the vibronic coupling properties of alkanes in XPS.

However, these calculations have several limitations. The ECH-DFT approach[31] does not yield accurate absolute binding energies. While theoretical spectra can be adjusted to align with experimental values, this method is unsuitable for predictions when no corresponding experiments exist (e.g., heptane). Additionally, our earlier calculations[21] indicated that ECH produces less accurate vibronic profiles compared to FCH, particularly regarding the intensity ratios of the 0-0 and 0-1 peaks, F 01 / F 00 . In the alternative approach[26] using a nitrogen basis set and ECP, scaling issues arise with the basis sets, ECPs, and vibrational frequencies (with different factors for different modes), leading to empirical adjustments that may not be valid for other systems. Furthermore, the linear coupling method (product of Poisson distributions) does not account for mode mixing effects (i.e., Duschinsky rotation). Lastly, while analyses have been conducted on core hole-induced structural changes and potential energy curves, there is a lack of comprehensive examination of active vibrational modes.

The goal of this study is twofold. First, we aim to evaluate the performance of a more accurate FCH-DR approach[18, 20, 21] for alkane molecules. This computational protocol has been validated[18-22] to accurately predict both binding energies (within 1 eV) and vibronic fine structures for XPS of cyclic compounds, including PAHs. The computational procedure employs only standard basis sets and model core potentials (MCPs) without any empirical parameter scaling. Second, with accurate theoretical spectra at hand, we will investigate general vibronic coupling rules in gaseous linear alkanes. This comprehensive analysis of aliphatic hydrocarbons, alongside our study of PAHs,[20] will enhance our understanding of hydrocarbon systems and contribute to our ongoing goal of constructing a high-resolution XPS library.

## II. COMPUTATIONAL METHODS

The FCH-DR method[18, 20, 21] was employed for the XPS computations in this work. Details of the methodology have been outlined previously. Briefly, DFT-based electronic structure calculations were performed using the B3LYP functional[35, 36] and the GAMESS software package[37]. A double basis set approach[21] was utilized for geometry optimization and frequency calculations, involving distinct basis sets for the ground state ( χ GS ) and core-ionized state ( χ FCH ). Specifically, the cc-pVTZ basis set was used for χ GS , while for χ FCH , the IGLO-III basis set[38] was used for core-ionized carbon atoms (C ∗ )

and the MCP/NOSeC-V-TZP basis set[39-41] was used for non-excited carbon atoms, complemented by the corresponding MCPs. This choice of basis sets has been validated for its accuracy in spectral simulations of various cyclic molecules[18, 21]. Fine integration grid and a stringent gradient convergence tolerance of 1 0 . × 10 -7 Hartree/Bohr were used. The optimized structures and vibrational modes for both the ground and FCH states were used as input to simulate FCFs using our in-house software package, DynaVib[42], modified[21] to extend its capabilities from the optical to the X-ray regimes. The DR method[24] was employed in the computation of FCFs, which characterizes the differences between the normal coordinates of the initial and final states. The 〈 0 | n 〉 amplitudes were computed through an iterative procedure[43, 44] initiated with 〈 0 0 | 〉 . Lorentzian broadening with a half width at half maximum (HWHM) γ =0.05 eV (i.e., the C1s core hole lifetime[45]) was applied to all C1s ionizations, and a small uniform shift δ rel =0.2 eV[46] was added to account for the scalar relativistic effect. The total XPS profile of each alkane is a weighted combination of the XPS spectrum from its nonequivalent carbon atoms. To better compare with the experimental data,[26] theoretical spectrum of each alkane molecules was shifted by δ = -0.27 eV, except that a slightly different value of δ = -0.33 eV was applied for methane.

## III. RESULTS

## A. Overview and notations for peak labels

Figure 2 presents the simulated vibrationally-resolved C1s XPS spectra for all alkanes under study, compared to available experimental data[26] (note that no experimental data is available for heptane). As shown in Figs. 3-6, more detailed analyses-including atom-specific spectra, interpretations of vibronic transitions, and identifications of active modes-are provided for each pair of molecules with the same number of non-equivalent carbons, C H n 2 n +2 and C n +1 H 2 n +4 (where n = 1 3 5 7 , , , ). Figure 7 consolidates the atomic-specific contributions from the edge carbons C 1 and C 2 to highlight patterns in spectral evolution.

For methane and ethane, each containing only one nonequivalent carbon, we use Roman numerals i-iii and i-v to label their major peaks, respectively. The same numeral indicates the same vibronic assignment and is transferable between the two molecules. For other alkanes with multiple non-equivalent carbons, we use Latin letters A, B, C, etc., to label features in ascending order of energy in the total spectrum. These labels are molecule-specific and not necessarily transferable among all molecules, although they are indeed transferable for long alkanes with n =5-8.

## B. Methane and ethane

The T d symmetry molecule methane is the simplest alkane with a single carbon atom. Its experimental spectrum [Fig. 2(a)] exhibits three prominent peaks (i-iii) that are well-reproduced by theory [Fig. 2(b)]: peaks i (experimental, 290.70 eV; theory, 290.71 eV ), ii (experimental, 291.10 eV; theory, 291.12 eV), and iii (experimental, 291.51 eV; theory, 291.55 eV). One can see that the experiment-theory deviation in peak energies is less than 0.1 eV. These peaks correspond to the 0-0, 0-1, and 0-2 transitions of the totally-symmetric (A ) C 1 ∗ -H stretching mode ν (6) (3383 cm -1 ) at the final-state geometry min FCH , which is the only Franck-Condon-active mode [Fig. 3(a,c)]. For clarity, we label peaks ii and iii as 0-1 6 and 0-2 6 , respectively [Fig. 3(a)], with the mode index indicated as a subscript.[20] Beyond these peaks, no additional features are clearly identifiable in the experimental spectrum [Fig. 2(a), bottom], while the 0-3 6 peak is predicted at approximately 292.00 eV [Fig. 2(b), bottom; Fig. 3(a)].

The D 3d symmetry molecule ethane also contains one single non-equivalent carbon atom. Its theoretical spectrum [Fig. 2(b)] closely matches experiment [Fig. 2(a)]. The main peak i (0-0 peak) exhibits a red shift of ca. 0.15 eV compared to methane. Peaks ii-iii are again observed, assigned to the C ∗ -H stretching mode ν (16) (3428 cm -1 ) [Figs. 3(d), right]. This mode is essentially similar to the only active mode of methane, and the vibrational frequency is slightly (45 cm -1 ) larger than that of methane. Besides, two additional peaks, iv (experimental, 290.73 eV; theory, 290.76 eV) and v (experimental, 291.36 eV; theory, 291.19 eV), are resolved and primarily attributed to the C ∗ -H bending (scissoring) mode ν (7) (1473 cm -1 ), which also includes contributions from the bending vibrations of the other -CH 3 group of the nonexcited carbon [Fig. 3(d), left]. Specifically, peaks iv and v are assigned as the 0 → 1 7 and 0 → 1 ,1 7 16 transitions, respectively.

## C. Propane and butane

Propane with C 2v symmetry contains two nonequivalent carbons, distal carbon C 1 and central carbon C , in a 2:1 ratio. 2 In the total theoretical spectrum [Fig. 2(b)], it now contains rich features (labeled by A-H), since C 1 and C 2 curves are similar in shape but with different 'phase' [see also Fig. 4(a-b)]. Feature A at 290.39 eV and feature B at 290.50 eV are 0-0 transitions by C 1 and C , respectively. 2 Corresponding experimental values are 290.40 and 290.50 eV, respectively, almost equivalent to our predictions. Features C and D come from 0-1 12 and 0-1 14 transitions for C ∗ -H bending modes ν (12) [1412 cm -1 , Fig. 4(e)] for C 1 and ν (14) [1488 cm -1 , Fig. 4(f)] for C 2 , respectively. Meanwhile, peak features E and F come from 0-1 25 and 0-1 26 transitions for C ∗ -H stretching modes ν (25) [3430 cm -1 , Fig. 4(e)] for C 1 and

ν (26) [3465 cm -1 , Fig. 4(f)] for C , respectively. 2 Features G and H are both assigned as 0-2 transitions with contributions of both carbons. Feature G is assigned as mixed 0-1 transitions for both the bending and stretching modes; while feature H is assigned as pure 0-2 transitions of the stretching mode.

Similarly, butane with C 2h symmetry also has two nonequivalent carbons, C 1 and C , in a 2:2 ratio. 2 Compared to propane, the fine structure looks simper because in the C 1 and C 2 contributions, all peaks almost coincide with each other [see also Fig. 4(c-d)]. The active mode analyses for butane [Fig. 4(g-h)] mirror the discussions for propane, identifying a C -H stretching and a C -H ∗ ∗ bending modes for each carbon. Peak A is contributed by 0-0 transitions; peaks B and C are 0-1 transitions for the bending and stretching modes, respectively; peaks D and E are 0-2 transitions, either mixed 0-1 transitions for both the bending and stretching modes (feature D) or pure 0-2 transitions of the stretching mode.

## D. Pentane and hextane

Pentane (C 2v ) and hexane (C 2h ) both have three nonequivalent carbon atoms, C -C , with weight ratios of 1 3 2:2:1 and 2:2:2, respectively. Both molecules exhibit similar spectra with six major features A-F. Peak A is assigned as 0-0 transition by the central carbon C 3 . Peak B is from 0-0 transitions of C 1 and C 2 as well as 01 transition of C 3 for a C -H bending mode (Fig. ∗ 5). Peaks C and D are assigned as 0-1 transitions, and Peaks E and F as 0-2 transitions. Active modes of the two molecules are analyzed in Fig. 5 (g)-(i) and Fig. 5 (j)-(l). Each atom-specific spectrum is identified with a mediumfrequency (1400-1500 cm -1 ) C -H bending mode and a ∗ high-frequency (3400-3500 cm -1 ) C -H stretching mode. ∗

## E. Heptane and octane

Heptane (C 2v ) and octane (C 2h ) have four nonequivalent carbon atoms, with weight ratios of 2:2:2:1 and 2:2:2:2, respectively. Although one additional nonequivalent carbon is present as compared to pentane and hextane, the total spectra do not change much, also exhibiting six features A-F. In other words, when n changes from 5 to 8, the vibronic profile almost converges into a consistent pattern. Concerning binding energies, the spectral series exhibit a slight red shift with increasing n . Despite the lack of experimental spectrum for heptane, our predicted theoretical spectrum matches the overall evolutionary trend.

The underlying reason is that in heptane and octane with relatively long length, spectra contributed by the two central carbons C 3 and C 4 almost coincide with each other, indicating their similar local chemical environment. Their 0-0 transitions both contribute to feature A in the total spectrum. 0-0 transitions of distal carbons

C 1 and C 2 make contributions to feature B some 0.2 eV beyond. The deviation of 0-0 transition energies for C 1 and C 2 makes peak B broader than peak A. In a even longer linear alkane, it is reasonable to yield that the spectral contributions for two neighboring carbons near the terminal deviates more than two neighboring carbons in the central region. Based on each atom-specific spectrum of both molecules [Fig. 6 (a)-(h)], each was identified with two active modes, including one mediumfrequency (1400-1500 cm -1 ) C -H bending mode and a ∗ high-frequency (3400-3500 cm -1 ) C -H stretching mode ∗ [Fig. 6 (i)-(p)].

## IV. DISCUSSION

## A. Red shift of the lowest peak A

In the total theoretical spectra [Fig. 2(b)], the lowestenergy peak, peak A, exhibits a red shift, with energies of 290.70 eV (methane), 290.58 eV (ethane), 290.38 eV (propane), 290.32 eV (butane), 290.12 eV (pentane), 290.09 eV (hexane), 290.07 eV (heptane), and 290.05 eV (octane), as indicated by the numbers at the top of this panel. The magnitude of the red shift is relatively large for alkanes with n =1 to 4 but becomes significantly smaller starting from n =5. As discussed earlier, for n =5 to 8, peak A for each molecule shares the same assignment (from central carbons). The energy differences between neighboring alkanes are only 0.02 to 0.03 eV (approximately 1 milli-Hartree), indicating that the values are approaching convergence.

## B. Ratio of distal and central carbons controls intensity ratios of peaks A and B

Above analyses have shown stable spectral profiles of six features A-F for linear alkanes at n = 5 -8 . Peak A and B are 0-0 peaks of central (C 3 when n = 5 -6 ; C 3 , C 4 when n = 7 -8 ) and distal (C 1 , C 2 ) carbons, respectively. As n increases, the intensity ratio of peak A to peak B increases due to the increase of more central carbons. This can be explained by the locality of core excitations and the alkane structure (without delocalized π structure), which make the atom-specific spectra dominantly determined by the nearest neighboring units. In a long alkane molecule, C 1 directly links to a -CH 2 group, C 2 directly bonds to a -CH 3 and -CH 2 groups on each side, and all other non-equivalent carbons all bond with two -CH 2 groups, which are classified as central carbons. It is worth noting that different central carbons contribute closely similar spectra, as verified by the coincidence of C 3 and C 4 spectra in heptane and octane. Therefore, in an even longer linear alkane, peak A will be extremely high compared to peak B. The total spectra can be considered as contributed from three types of carbons: C , 1 C , and central carbons. This explains the stable features 2

in longer linear alkanes, with their relative peak intensities adjusted by the ratios of different types of carbons and their relative BE differences.

## C. Evolution of distal carbons C 1 and C 2

Figure 7 clearly demonstrates the evolution of vibrationally-resolved XPS profiles for the two distal carbon atoms C 1 (orange) and C 2 (lightblue). The concept of distal carbons is evident for long alkanes with n =58, while all carbons are considered as distal carbons for short alkanes with n =1-4, providing a comprehensive and smooth picture. Methane shows three features iiii, attributed to the C -H stretching vibration. ∗ Since ethane, an additional C ∗ -H bending vibration is involved, contributing to peaks vi (pure 0-2 bending mode excitation) and v (mixed 0-1 excitations with both modes). All other atom-specific spectral profiles show great similarity to the ethane spectrum, with relative peak intensities of the five peaks preserved almost the same. Figure 7 highlights the evolution of peaks i-iii. These results indicate that the spectrum of ethane can be considered as a 'spectral seed', as different shifted spectra of ethane approximately generate the total spectra of a longer linear alkane. Additionally, a switch of BE order is found for C 1 and C . Carbon C 2 2 starts since n =3, exhibiting a higher binding energy of ca. 0.1 eV than C . 1 This energy order is quickly reversed since n =4, and at n =8, the binding energy of C 2 is ca. 0.1 eV lower than C . 1

## D. Structural changes induced by core ionizations

Structurally, the linear alkanes are generated by gradually adding methylene groups (-CH -) to methane. Each 2 molecule is constructed by linking C-C and C-H single bonds, and studying the structural change induced by the C1s core hole is crucial. In alkanes, the typical C-C bond length hovers around 1.529 Å, with a tiny difference between distal (1.528 Å) and central carbons (1.529-1.530 Å), as shown in Fig. 1. The structural changes near the excited carbon are summarized in Table I, revealing discernible patterns. Specifically, the bond lengths C ∗ -C r and C -C ∗ l in the FCH state exhibit minimal changes (not exceeding 0.02 Å) compared to the ground state structure, indicating marginal influence due to the presence of a C1s core hole. Here, C l (C ) is the left (right) carr bon atom bonded to C ∗ (see relative atomic positions in Fig. 1). The angle ∠ C -C -C , typically close to 114 , r ∗ l ◦ experiences an increment of ca. 1 ◦ from the ground- to the final-state geometry. This angle at C ∗ experiences slightly larger bending for central carbons compared to distal carbons. All C -H bonds exhibit virtually identi-∗ cal characteristics in both the ground and the final states, with minor shortening of 0.05-0.07 Å.

## V. SUMMARY AND CONCLUSIONS

To summarize, we simulated high-resolution, vibrationally-resolved C1s X-ray photoelectron spectra (XPS) for a series of linear alkanes, ranging from methane to octane (C H n 2 n +2 , with n =1-8), by combining full core-hole density functional theory with Franck-Condon calculations that account for Duschinsky rotation. Our theoretical spectra align excellently with all available experimental results, and we predicted the spectra of heptane, showing reasonable spectral evolution in this series. This approach has facilitated reliable peak assignments and enabled an investigation into the spectral evolution with increasing chain length, as well as the underlying vibronic coupling rules governing this molecular family.

This work extends our FCH-DR approach to flexible molecular systems, enhancing our understanding of vibronic coupling phenomena. The good performance of our protocol in these flexible systems should benefit from the minimal structural changes induced by core ionization. It also proves that the harmonic oscillator approximation works effectively for alkane systems. The insights gained from this analysis of simple alkane systems will be valuable for studying more complex hydrocarbons and other flexible molecules. Ultimately, this study contributes to our ongoing research goal of constructing a comprehensive, high-resolution theoretical XPS library.

The spectrum of ethane serves as a 'spectral seed', with the vibronic profile of each longer alkane exhibiting similar characteristics, albeit with shifts in peak positions and slight intensity adjustments. The carbon atoms in these longer alkanes are categorized as central or distal (C 1 and C ), with central carbons contribut2 ing nearly identically to the lowest-energy feature A, while distal carbons contribute to the second-lowest feature B. Detailed assignments elucidate the distinct spectra observed in short alkanes ( n =1-4) and their stabilization in longer alkanes ( n =5-8). Through comprehensive assignments, we tracked the evolution of specific spectral features across the molecular series. We investigated core ionization-induced structural changes and identified Franck-Condon-active vibrational modes, highlighting two key active modes present in each atomspecific spectrum: the C -H symmetric stretching (ca. ∗ 3400-3500 cm -1 ) and the C ∗ -H bending (ca. 1400-1500 cm -1 ) modes.

## ACKNOWLEDGMENTS

We thank Prof. Yi Luo for helpful discussions. Financial support from the National Natural Science Foundation of China (Grant No. 12274229) is greatly acknowledged. M.W. thanks to Fund for Fostering Talented Doctoral Students of Nanjing University of Science and Technology.

TABLE I: Local structural changes of each alkane near the core ionization center (C ∗ ) at the optimized FCH state structure ( min FCH ) as compared to its ground state structure ( min GS ). Bond lengths are in Å and angles in ◦ .

| Molecule   | C ∗   |   C ∗ -H | C ∗ -H     | C ∗ -C r a   | C ∗ -C r a   | C ∗ -C l a   | C ∗ -C l a   | ∠ C r -C ∗ -C l   | ∠ C r -C ∗ -C l   |
|------------|-------|----------|------------|--------------|--------------|--------------|--------------|-------------------|-------------------|
| Methane    | C 1   |    1.037 | (-0.051) b |              |              |              |              |                   |                   |
| Ethane     |       |    1.035 | (-0.056)   | 1.526        | (-0.002)     |              |              |                   |                   |
| Propane    | C 1   |    1.034 | (-0.058)   | 1.543        | (+0.015)     |              |              |                   |                   |
| Propane    | C 2   |    1.033 | (-0.061)   | 1.515        | (-0.013)     | 1.515        | (-0.013)     | 113.884           | (+0.878)          |
| Butane     | C 1   |    1.034 | (-0.056)   | 1.543        | (+0.015)     |              |              |                   |                   |
| Butane     | C 2   |    1.033 | (-0.062)   | 1.530        | (+0.001)     | 1.512        | (-0.016)     | 114.302           | (+0.933)          |
| Pentane    | C 1   |    1.034 | (-0.058)   | 1.544        | (+0.015)     |              |              |                   |                   |
| Pentane    | C 2   |    1.033 | (-0.062)   | 1.530        | (+0.000)     | 1.511        | (-0.017)     | 114.447           | (+1.109)          |
| Pentane    | C 3   |    1.033 | (-0.063)   | 1.526        | (-0.003)     | 1.526        | (-0.003)     | 114.780           | (+1.019)          |
| Hexane     | C 1   |    1.034 | (-0.058)   | 1.544        | (+0.016)     |              |              |                   |                   |
| Hexane     | C 2   |    1.033 | (-0.062)   | 1.530        | (+0.001)     | 1.511        | (-0.017)     | 114.483           | (+1.160)          |
| Hexane     | C 3   |    1.033 | (-0.063)   | 1.526        | (-0.004)     | 1.526        | (-0.004)     | 114.925           | (+1.192)          |
| Heptane    | C 1   |    1.034 | (-0.058)   | 1.544        | (+0.016)     |              |              |                   |                   |
| Heptane    | C 2   |    1.033 | (-0.062)   | 1.531        | (+0.001)     | 1.511        | (-0.018)     | 114.494           | (+1.168)          |
| Heptane    | C 3   |    1.033 | (-0.063)   | 1.526        | (-0.003)     | 1.525        | (-0.005)     | 114.957           | (+1.236)          |
| Heptane    | C 4   |    1.033 | (-0.063)   | 1.525        | (-0.004)     | 1.525        | (-0.004)     | 115.059           | (+1.356)          |
| Octane     | C 1   |    1.034 | (-0.058)   | 1.545        | (+0.016)     |              |              |                   |                   |
| Octane     | C 2   |    1.033 | (-0.062)   | 1.531        | (+0.001)     | 1.511        | (-0.018)     | 114.511           | (+1.181)          |
| Octane     | C 3   |    1.033 | (-0.063)   | 1.527        | (-0.003)     | 1.525        | (-0.005)     | 114.966           | (+1.241)          |
| Octane     | C 4   |    1.033 | (-0.063)   | 1.526        | (-0.004)     | 1.524        | (-0.005)     | 115.088           | (+1.397)          |

- a C (C ) is the left (right) carbon atom bonded to C l r ∗ (see relative atomic positions in Fig. 1).
- b Numbers in parentheses denote the relative values compared with those at the ground state.

FIG. 1: Optimized ground-state structures (i.e., min GS ) of the eight linear alkanes under study: (a) methane, (b) ethane, (c) propane, (d) butane, (e) pentane, (f) hextane, (g) heptane, and (h) octane. The chemical formula, symmetry, all non-equivalent C sites, and all C-C bond lengths (in Å) are specified for each molecule.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00589v1-figure-05f6f060-1dd5-45a1-8b21-6e3ef3724b4d.jpg)

FIG. 2: Vibrationally-resolved C1s XPS spectra of linear alkane molecules. (a) Experimental spectra recaptured from Karlsen et al.[26] (note that heptane has no available experiment). (b) Our spectral results simulated by using the FCH-DR method.[18, 20, 21] For methane and ethane with one non-equivalent carbon, major peaks are labeled by Roman numerals and interpreted by active modes (visualized in Fig. 3). For propane to octane with multiple non-equivalent carbons, atom-specific contributions from carbons C n (defined in Fig. 1) are shown with dotted curves, and major peaks in each total spectrum are labeled with Latin letters. Vertical dashed lines (gray) indicate positions of peak A, with values (in eV) given explicitly as a table above panel b.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00589v1-figure-66c3e73a-94a4-4de0-bf83-f8b53162f1e1.jpg)

## Methane 0\_n transitions

FIG. 3: Vibrational analysis for methane and ethane. (a,b) Simulated C1s XPS spectra of methane and ethane with vibronic transitions interpreted for major peaks. For instance, in the notation 0 → n k , n stands for the vibrational quantum number and k denotes the mode index. Roman numerals i-v recapture peak labels defined in Fig. 2(b). (c,d) Franck-Condon-active vibrational mode(s) interpreted at the optimized fine-state structure ( min FCH ) of methane and ethane. Main vibrational characteristics near C ∗ (the core ionization center) is specified under each frequency.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00589v1-figure-03ca24fa-0dc2-4a5d-9cf7-051227bfaabd.jpg)

Propane 0\_n transitions

FIG. 4: Vibrational analysis for propane and butane. (a-d) Atom-specific contributions of simulated C1s XPS spectra of propane and butane. Underlying vibronic transitions are interpreted for major peaks. Vertical lines (green) marked with Latin letters A,B, C, · · · denote positions of peaks in the total theoretical spectra [Fig. 2(b)]. (e-h) Active vibrational modes interpreted at the optimized fine-state structures ( min FCH ) of propane and butane. Main characteristics near C ∗ is specified under each frequency. Major peaks for C 1 (labeled as i-v) are interpreted.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00589v1-figure-7832b825-dd18-45ea-b7ca-ffef9e24e567.jpg)

## Pentane 0\_n transitions

FIG. 5: Vibrational analysis for pentane and hextane. (a-f) Atom-specific contributions of simulated C1s XPS spectra of pentane and hextane. Underlying vibronic transitions are interpreted for major peaks. Vertical lines (green) marked with Latin letters A,B, C, · · · denote positions of peaks in the total theoretical spectra [Fig. 2(b)]. (g-l) Active vibrational modes interpreted at the optimized fine-state structures of pentane and hextane. Main characteristics near C ∗ is specified under each frequency.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00589v1-figure-420f74fe-ba63-4ffc-84ab-4e7ecc9dca7d.jpg)

## Heptane 0\_~n transitions

Binding Energy (eV)

FIG. 6: Vibrational analysis for heptane and octane. (a-h) Atom-specific contributions of simulated C1s XPS spectra of heptane and octane. Underlying vibronic transitions are interpreted for major peaks. Vertical lines (green) marked with Latin letters A,B, C, · · · denote positions of peaks in the total theoretical spectra [Fig. 2(b)]. (i-p) Active vibrational modes interpreted at the optimized fine-state structures of heptane and octane. Main vibrational characteristics near C ∗ is specified.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00589v1-figure-2be0f80e-f262-4fc8-9bbc-131ab16e53f5.jpg)

FIG. 7: Simulated vibrationally-resolved C1s XPS spectra of distal carbons C 1 (orange) and C 2 (sky blue) in linear alkanes, showing spectral evolution with molecule length. Peaks i-iii are labeled for methane C 1 , and peaks i-v for other molecules. Yellow lines highlight selected C 1 peak trends. C 2 spectra are similar, so no labeling is shown. Molecule structures are schematically depicted in the top left of each panel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00589v1-figure-4c93d16b-6e9f-4d2f-91a9-4d817256b2f7.jpg)

- [1] P. van der Heide, X-Ray Photoelectron Spectroscopy: An Introduction to Principles and Practices (John Wiley &amp; Sons, Inc., Hoboken, NJ, USA, 2011). [26] T. Karlsen, K. J. Børve, L. J. Sæthre, K. Wiesner, M. Bässler, and S. Svensson,
- J. Am. Chem. Soc. 124 , 7866 (2002).
- [3] V. Carravetta, R. C. Couto, and H. Ågren, J. Phys.: Condens. Matter 34 , 363002 (2022).
- [2] U. Gelius, S. Svensson, H. Siegbahn, E. Basilier, Å. Faxälv, and K. Siegbahn, Chem. Phys. Lett. 28 , 1 (1974).
- [27] E. Stauffer, J. A. Dolan, and R. Newman, in Fire Debris Analysis (Elsevier, 2008) pp. 49-83.
- [4] U. Hergenhahn, J. Phys. B: At. Mol. Opt. Phys. 37 , R89 (2004). ]
- [28] R. Atkinson, J. Arey, and S. M. Aschmann, Atmos. Environ. 42 , 5859 (2008).
- [5] S. Svensson, J. Phys. B: At. Mol. Opt. Phys. 38 , S821 (2005).
- [29 Y. Xiang, H. Wang, J. Cheng, and J. Matsubu, Catal. Sci. Technol. 8 , 1500 (2018).
- [6] K. Siegbahn, ESCA applied to free molecules (NorthHolland Publishing, 1969).
- [30] A. M. Gaffney and O. M. Mason, Catal. Today 285 , 159 (2017).
- [8] M. Mendolicchio, A. Baiardi, G. Fronzoni, M. Stener, C. Grazioli, M. de Simone, and V. Barone, J. Chem. Phys. 151 , 124105 (2019).
- [10] J.-E. Rubensson, A. Pietzsch, and F. Hennies,
- [7] K. Ueda, J. Phys. B: At. Mol. Opt. Phys. 36 , R1 (2003).
- [31] T. D. Thomas, L. J. Saethre, S. L. Sorensen, and S. Svensson, J. Chem. Phys. 109 , 1041 (1998).
- [9] F. Montorsi, F. Segatta, A. Nenov, S. Mukamel, and M. Garavelli, J. Chem. Theory Comput. 18 , 1003 (2022).
- J. Electron Spectros. Relat. Phenomena 185 , 294 (2012).
- [11] C. Miron and P. Morin, in Handbook of High-resolution Spectroscopy , edited by M. Quack and F. Merkt (John Wiley &amp; Sons, Ltd, Chichester, UK, 2011) p. hrs066.
- [12] J. R. Rumble Jr., D. M. Bickham, and C. J. Powell, Surf. Interface Anal. 19 , 241 (1992).
- [13] C. D. Wagner, Handbook of x-ray photoelectron spectroscopy: a reference book of standard data for use in x-ray photoelectron spectroscopy (Perkin-Elmer, Shelton, CT, 1979).
- [14] W. Jolly, K. Bomben, and C. Eyermann, Atom. Data Nucl. Data 31 , 433 (1984).
- [15] Https://xpsdatabase.com/, Accessed on 2024-7-20.
- [17] Http://www.lasurface.com/xps/index.php/, Accessed on 2024-7-20.
- [16] Http://www.sasj.jp/, Accessed on 2024-7-20.
- [18] M. Wei, X. Cheng, L. Zhang, J.-R. Zhang, S.-Y. Wang, G. Ge, G. Tian, and W. Hua, Phys. Rev. A 106 , 022811 (2022).
- [20] X. Cheng, M. Wei, G. Tian, Y. Luo, and W. Hua, J. Phys. Chem. A 126 , 5582 (2022).
- [19] M. Wei, L. Zhang, G. Tian, and W. Hua, Phys. Rev. A 108 , 022816 (2023).
- [21] W. Hua, G. Tian, and Y. Luo, Phys. Chem. Chem. Phys. 22 , 20014 (2020).
- [23] M. Wei, J. Zuo, G. Tian, and W. Hua, arXiv preprint arXiv:2403.09109 (2024).
- [22] M. Wei, J. Zuo, G. Tian, and W. Hua, Phys. Rev. A 109 , 022820 (2024).
- [24] F. Duschinsky, Acta Physicochim. URSS 7 , 551 (1937), and its English translation by C. W. Müler available from http://www.chem.purdue.edu/zwier/pubs/Duschinsky.pdf.
- [25] L. Zhang, M. Wei, G. Ge, and W. Hua, Phys. Rev. A 109 , 032815 (2024).
- [33] D. W. Davis and D. A. Shirley, Chem. Phys. Lett. 15 , 185 (1972).
- [32] W. L. Jolly and D. N. Hendrickson, J. Am. Chem. Soc. 92 , 1863 (1970).
- [34] O. Plashkevych, T. Privalov, H. Ågren, V. Carravetta, and K. Ruud, Chem. Phys. 260 , 11 (2000).
- [35] A. D. Becke, J. Chem. Phys. 98 , 5648 (1993).
- [36] P. J. Stephens, F. J. Devlin, C. F. Chabalowski, and M. J. Frisch, J. Phys. Chem. 98 , 11623 (1994).
- [38] W. Kutzelnigg, U. Fleischer, and M. Schindler, in Deuterium and shift calculation (Springer, 1990) pp. 165-262.
- [37] G. M. J. Barca, C. Bertoni, L. Carrington, D. Datta, N. De Silva, J. E. Deustua, D. G. Fedorov, J. R. Gour, A. O. Gunina, E. Guidez, T. Harville, S. Irle, J. Ivanic, K. Kowalski, S. S. Leang, H. Li, W. Li, J. J. Lutz, I. Magoulas, J. Mato, V. Mironov, H. Nakata, B. Q. Pham, P. Piecuch, D. Poole, S. R. Pruitt, A. P. Rendell, L. B. Roskop, K. Ruedenberg, T. Sattasathuchana, M. W. Schmidt, J. Shen, L. Slipchenko, M. Sosonkina, V. Sundriyal, A. Tiwari, J. L. Galvez Vallejo, B. Westheimer, M. Wloch, P. Xu, F. Zahariev, and M. S. Gordon, J. Chem. Phys. 152 , 154102 (2020).
- [39] Y. Sakai, E. Miyoshi, M. Klobukowski, and S. Huzinaga, J. Chem. Phys. 106 , 8084 (1997).
- [41] Http://sapporo.center.ims.ac.jp/sapporo/.
- [40] T. Noro, M. Sekiya, and T. Koga, Theor. Chem. Acc. 98 , 25 (1997).
- [42] G. Tian, S. Duan, W. Hua, and Y. Luo, DynaVib, version 1.0, Royal Institute of Technology: Sweden, 2012.
- [43] P. T. Ruhoff, Chem. Phys. 186 , 355 (1994).
- [44] P. T. Ruhoff and M. A. Ratner, Int. J. Quant. Chem. 77 , 383 (2000).
- [46] L. Triguero, O. Plashkevych, L. Pettersson, and H. Ågren,
- [45] C. Nicolas and C. Miron, J. Electron. Spectrosc. Relat. Phenom. 185 , 267 (2012).
- J. Electron Spectros. Relat. Phenomena 104 , 195 (1999).
Xiao Cheng Minrui Wei Guangjun Tian Weijie Hua
Simulating Vibrationally-Resolved X-ray Photoelectron Spectra of Flexible Molecules: Linear Alkanes C$_{n}$H$_{2n+2}$ ($n$=1-8)
We integrated full core-hole density functional theory with Franck-Condon
calculations, considering Duschinsky rotation, to simulate
vibrationally-resolved C1s X-ray photoelectron spectra (XPS) of eight linear
alkanes, from methane to octane (C$_{n}$H$_{2n+2}$, $n$=1--8). Results align
excellently with experimental absolute binding energies and profiles. The
spectrum of ethane serves as a ``spectral seed'', with each longer alkane's
atom-specific spectrum displaying similar characteristics, albeit with shifts
and slight intensity adjustments. Detailed assignments clarify the distinct
spectra in short alkanes ($n$=1--4) and their stabilization in long alkanes
($n$=5--8). Carbons are classified as central or distal (C$_1$ and C$_2$), with
central carbons contributing nearly identically to the lowest-energy feature A,
while distal carbons contribute to the second lowest-energy feature B (both are
0-0 transitions). Increasing molecule size adds more central carbons, enhancing
feature A and weakening feature B. Our analysis identifies two dominant
Franck-Condon-active vibrations: C$^*$--H stretching ($\sim$3400-3500
cm$^{-1}$) and bending ($\sim$1400-1500 cm$^{-1}$) modes. Structural analysis
shows that core ionization minimally affects alkane geometry, less than in ring
compounds from previous studies. This work extends our protocol from rigid ring
compounds to flexible molecules, contributing to building a high-resolution
theoretical XPS library and enhancing the understanding of vibronic coupling.